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Abstract 

Magnetic fields appear everywhere in the universe. Prom stars and galaxies, all the 
way to galaxy clusters and remote protogalactic clouds, magnetic fields of consid- 
erable strength and size have been repeatedly observed. Despite their widespread 
presence, however, the origin of cosmic magnetic fields is still a mystery. The galac- 
tic dynamo is believed capable of amplifying weak magnetic seeds to strengths like 
those measured in ours and other galaxies. But the question is where do these seed 
fields come from? Are they a product of late, post-recombination, physics or are they 
truly cosmological in origin? The idea of primordial magnetism is attractive because 
it makes the large-scale magnetic fields, especially those found in early protogalactic 
systems, easier to explain. As a result, a host of different scenarios have appeared 
in the literature. Nevertheless, early magnetogenesis is not problem free, with a 
number of issues remaining open and matter of debate. We review the question of 
the origin of primordial magnetic fields and consider the limits set on their strength 
by the current observational data. The various mechanisms of pre-recombination 
magnetogenesis are presented and their advantages and shortcomings are debated. 
We consider both classical and quantum scenarios, that operate within as well as 
outside the standard model, and also discuss how future observations could be used 
to decide whether the large-scale magnetic fields we see in the universe today are 
truly primordial or not. 

Key words: Magnetogenesis, Cosmological Magnetic Fields. 



1 Introduction 



Observations have well established the widespread presence of magnetic fields in the uni- 
verse [1,2,3,4,5,6,7]. In fact, as the technology and the detection methods improve, it seems 
that magnetic fields are everywhere. The Milky Way, for example, possesses a coherent B- 
field of /iG strength over the plain of its disc. These fields are a very important component 
of the interstellar medium, since they govern the gas-cloud dynamics, determine the energy 
of cosmic rays and affect star formation. Similar magnetic fields have also been detected in 
other spiral and barred galaxies. Cosmic magnetism is not confided to galaxies only how- 
ever. Observations have repeatedly verified 5-fields of /iG-order strength in galaxy clusters 
and also in high redshift protogalactic structures. Recently, in particular, Kronberg et al and 
Bernet et al reported organized, strong 5-fields in galaxies with redshifts close to 1.3 [8,9]. 
Also, Wolfe et al have detected a coherent magnetic field of approximately 100 /iG in a galaxy 
at z ~ 0.7 [10]. All these seem to suggest that magnetic fields similar to that of the Milky 
Way are common in remote, high-redshift galaxies. This could imply that the time needed by 
the galactic dynamo to build up a coherent 5-field is considerably less than what is usually 
anticipated. On the other hand, the widespread presence of magnetic fields at high redshifts 
may simply mean that they are cosmological (pre-recombination) in origin. Although it is still 
too early to reach a conclusion, the idea of primordial magnetism gains ground, as more fields 
of micro-Gauss strength are detected in remote protogalaxies. Further support may also come 
from very recent reports indicating the presence of coherent magnetic fields in the low density 
intergalactic space, where typical dynamo mechanisms cannot operate, with strengths close to 
10 -15 G [11,12,13,14,15]. The measurements of [12], in particular, are based on halos detected 
around Active Galactic Nuclei (AGN) observed by the Fermi Gamma-Ray Space Telescope. 
Complementary studies seem to limit the strength of these 5-fields between ~ 10~ 17 and 
~ 10~ 14 Gauss [15]. Analogous lower limits were also reported by [11,14] and [13], after mea- 
suring radiation in the GeV band (7-rays) produced by the interaction of TeV photons from 
distant blazars with those of the Cosmic Microwave Background (CMB). If supported by fu- 
ture surveys, these measurements will render considerable credence to the idea of primordial 
magnetism. It is possible, however, that the matter will not settle unless magnetic imprints 
are found in the CMB spectrum [16]. 

Among the attractive aspects of cosmological magnetic fields, is that they can in principle ex- 
plain all the large-scale fields seen in the universe today [17,18,19,20,21,22]. Nevertheless, early 
magnetogenesis is not problem free. The galactic dynamo needs an initial magnetic field in or- 
der to operate. Such seed fields must satisfy two basic requirements related to their coherence 
scale and strength [23,24,25,26,27,28,29]. The former should not drop below 10 Kpc, otherwise 
it will destabilise the dynamo. The latter typically varies between 10~ 12 and 10 -22 Gauss. 
It is conceivable, however, that in open, or dark-energy dominated, Friedmann- Robert son- 
Walker (FRW) cosmologies the minimum required magnetic strength could be pushed down 
to ~ 10~ 30 G [30]. Producing magnetic seeds that comply with the above mentioned specifi- 
cations, however, has so far proved a rather difficult theoretical exercise. There are problems 
with both the scale and the strength of the initial field. Roughly speaking, magnetic seeds 
generated between inflation and recombination have too small coherence lengths. The reason 



is causality, which confines the scale of the 5-field within the size of the horizon at the time 
of magnetogenesis. This is typically well below the dynamo requirements. If we generate the 
seed at the electroweak phase transition, for example, the size of the horizon is close to that of 
the Solar System. Assuming that some degree of turbulence existed in the pre-recombination 
plasma, one can increase the coherence scale of the initial field by appealing to a mechanism 
that in hydrodynamics is known as 'inverse cascade'. In magnetohydrodynamics (MHD), the 
process results from the conservation of magnetic helicity and effectively transfers magnetic 
energy from small to successively larger scales. The drawback is that inverse cascade seems to 
require rather large amounts of magnetic helicity in order to operate efficiently [31,32,33,34]. 
Inflation can solve the scale problem, since it naturally creates superhorizon-sized correlations. 
There, however, we have a serious strength issue. Magnetic fields that were generated during 
a period of typical de Sitter-type inflation are thought to be too weak to seed the galactic 
dynamo [35] . The solution to the strength problem is usually sought outside the realm of clas- 
sical electromagnetic theory, or of conventional FRW cosmology. There is a plethora of articles 
that do exactly that, although mechanisms operating within standard electromagnetism and 
the Friedmann models have also been reported in the literature. 

The aim of this review is to present the various mechanisms of early magnetogenesis, outline 
their basic features and discuss their advantages and weaknesses. The next section starts with 
a brief overview of the observation techniques that have established the ubiquitous presence of 
large-scale magnetic fields in the universe. We then provide the limits on cosmological -B-fields, 
imposed by primordial nucleosynthesis and the isotropy of the CMB. Section three sets the 
mathematical framework for the study of large-scale magnetic fields in relativistic cosmological 
models. There, for completeness, we also outline the typical magnetic effects on structure 
formation and how the latter could have backreacted on the .B-field itself. In section four we 
discuss cosmic magnetogenesis within the realm of standard electrodynamics and within the 
limits of conventional Friedmannian cosmology. After a brief review of the FRW dynamics, 
we explain why it is theoretically difficult to generate and sustain astrophysically relevant In- 
fields in these models. More specifically, why inflationary magnetic fields in FRW cosmologies 
are generally expected to have residual strengths less than 10 -50 G (far below the current 
galactic dynamo requirements) at the epoch of galaxy formation. At the same time, it is also 
pointed out that the standard picture can change when certain general relativistic aspects of 
the magnetic evolution are accounted for. More specifically, it is shown that curvature effects 
can in principle slow down the standard 'adiabatic' magnetic decay and thus lead to -B-fields 
with residual strengths much stronger than previously anticipated. 

In section five we describe the generation of magnetic fields by nonlinear and out-of-equilibrium 
processes, which are believed to have taken place in the early universe. We begin by analyzing 
several mechanisms of magnetogenesis that could have operated during the reheating epoch 
of the universe, namely parametric resonance, generation of stochastic electric currents and 
the breaking of the conformal invariance of the electromagnetic field by cosmological per- 
turbations. Then, we address the generation of magnetic fields during cosmological phase 
transitions. It is believed that at least two of such phase transitions have occurred in the early 
universe: the EW (Electroweak) and the QCD (Quantum Chromodynamical). In general, the 
problem with the (post-inflationary) early-universe magnetogenesis is that the generated B- 



fields have high intensity but very short coherent scale (in contrast to what happens during 
inflation), which amounts to performing certain line averages to obtain the desired large-scale 
intensities. This procedure generally results in weak magnetic fields. To a certain extent, the 
uncertainty in the obtained residual magnetic values reflects our limited knowledge of the 
dissipative processes operating at those times. Thus, better understanding of the reheating 
physics is required, if we are to make more precise predictions. Phase transitions, on the other 
hand, are better understood. Note, however, that despite the fact that the EW phase tran- 
sition in the standard model is a second order process, extensions to other particle physics 
models treat it as first order. The QCD phase transition, on the other hand, was recently 
established to be a smooth crossover [36]. To the best of our knowledge, however, no work on 
primordial magnetogenesis in this scenario has been reported in the literature. 

Section six provides an overview of magnetic generation mechanisms operating outside the 
standard model. In all scenarios the magnetic fields are created in the very early universe, 
during inflation. Then, subhorizon quantum fluctuations in the electromagnetic field become 
classical superhorizon perturbations, manifesting themselves as current supported magnetic 
fields during the subsequent epochs of standard cosmology. In order to overcome the problem 
of not creating strong enough magnetic seeds, which is known to plague standard electro- 
dynamics, different theories are explored. There are two basic classes of models, depending 
on whether electrodynamics is linear or nonlinear. In the first case, magnetic fields of astro- 
physically relevant strengths are usually achieved after breaking the conformal invariance of 
electromagnetism. This can be achieved by coupling the electromagnetic with a scalar field (as 
it naturally happens with the dilaton in string cosmology), by introducing dynamical extra 
dimensions, through quantum corrections leading to coupling with the curvature tensor, by 
inducing symmetry breaking, or by means of trace anomaly. When dealing with nonlinear 
electrodynamics, on the other hand, the conformal invariance of Maxwell's equations is nat- 
urally broken in four dimensions. Recall that the concept of nonlinear electrodynamics was 
first introduced by Born, in his search for a classical, singularity-free theory of the electron. 
Another example is provided by the description of virtual electron pair creation, which induces 
a self-coupling of the electromagnetic field. Linear and nonlinear models of electrodynamics 
are discussed in detail and the parameter space, for which strong enough magnetic seeds are 
generated, is determined. 

Finally, in section seven, we briefly summarise the current state of research on primordial 
magnetogenesis and take look at the future expectations. 



2 Magnetic fields in the universe 



Magnetic fields have long established their ubiquitous presence in the universe. They are a 
major component of the interstellar medium, contributing to the total pressure, affecting the 
gas dynamics, the distribution of cosmic rays and star formation. It also seems very likely 
that large-scale magnetic fields have played a fundamental role during the formation of galaxy 
clusters. Despite our increasing knowledge, however, many key questions related to the origin 
and the role of these fields remain as yet unanswered. 



2. 1 Large-scale magnetic fields in the universe 



Most galaxies, including the Milky Way, carry coherent large-scale magnetic fields of /iG-order 
strength. Analogous fields have also been detected in galaxy clusters and in young, high-redsift 
protogalactic structures. In short, the deeper we look for magnetic fields in the universe, the 
more widespread we find them to be. 



2.1.1 Detection and measuring methods 

The key to magnetic detection is polarized emission at the optical, the infrared, the submil- 
limeter and the radio wavelengths. Optical polarization is due to extinction along the line 
of sight, caused by elongated dust grains aligned by the interstellar magnetic field. The net 
result is that the electromagnetic signal has a polarisation direction parallel to the intervening 
B- field. This physical mechanism is sometimes referred to as the Davis-Greenstein effect [37]. 
Although optical polarization is of limited value, it has unveiled the magnetic structure in the 
spiral arms of the Milky Way and in other nearby galaxies [38,39,40,41]. 

Most of our knowledge about galactic and intergalactic magnetic fields comes from radio- 
wave signals. The intensity of synchrotron emission is a measure of the strength of the total 
magnetic field component in the sky plane. Note that polarized emission is due to ordered 
.B-fields and unpolarized comes from turbulent ones. The Zeeman splitting of radio spectral 
lines is the best method to directly measure the field strength in gas clouds of our galaxy [42] , 
OH masers in starburst galaxies [43] and in dense HI clouds in distant galaxies on the line of 
sight towards bright quasars [10]. The drawback is that the Zeeman effect is very weak and 
can mainly be used for detecting intersellar magnetic fields. This is due to the small line shift, 
which given by 

^ = 1^, (2.1.D 

v v 

and is extremely difficult to observe at large distances. Note that in the above the 5-field is 
measured in /zG and the line frequency in Hz. Also, the parameter g represents the Lande 
factor that relates the angular momentum of an atom to its magnetic moment. 

When polarized electromagnetic radiation crosses a magnetized plasma its orientation is 



changed by Faraday rotation. The latter is caused by the left and right circular polarisation 
states traveling with different phase velocities. For linearly polarised radiation, the rotation 
measure (RM) associated with a source at redshift z s is (cf. , e.g., [44]) 



RM {z a ) ~ 8 x fO 5 / .; "\; dL{z) (rad/m 2 ) , (2.1.2) 



o (1 + z) 2 

where n e is the electron density of the intervening plasma (in cm -3 ), B» is the magnetic 
intensity along the line of sight (in /iG) and dL is the distance traveled by the radio signal. 
The latter is given by 

dL (z) = 10- 6 l lJ_^ dz Mpc , (2.1.3) 

H y/l + \l Q Z 

with H$ and Qq representing the present values of the Hubble constant of the density param- 
eter respect ively|_j As the rotation angle is sensitive to the sign of the field direction, only 
ordered 5-fields can give rise to Faraday rotation. Multi-wavelength observations determine 
the strength and the direction of the line-of-sight magnetic component. Then, the total inten- 
sity and the polarization vectors yields the three-dimensional picture of the field and allow us 
to distinguish between its regular, anisotropic and random components. 

Some novel detection methods try to exploit the effects that an intervening magnetic field 
can have upon the highly energetic photons emitted by distant active sources [45,46,47,48]. 
Using such techniques, together with data from state-of-the-art instruments (like the Fermi 
Gamma- Ray Space Telescope for example), three independent groups have recently reported 
the detection of intergalactic magnetic fields with strengths close to 1CT 15 G (see § 2.1.2 next). 



2.1.2 Galactic and extragalactic magnetic fields 

The strength of the total magnetic field in galaxies can be determined from the intensity of 
the total synchrotron emission, assuming equipartition between the magnetic energy density 
and that of the cosmic rays 2 1. This seems to hold on large scales (both in space and time), 
though deviations occur locally. Typical equipartition strengths in spiral galaxies are around 
10/iG. Radio-faint galaxies, like the M31 and M33, have weaker total fields (with B ~ 5/xG), 
while gas-rich galaxies with high star-formation rates, such as the M51, M83 and NGC6946, 
have magnetic strengths of approximately 15//G. The strongest fields, with values between 50 
/iG and 100 /iG, are found in starburst and merging galaxies, like the M82 and NGC4038/39 
respectively [51]. 



1 Conventionally, positive RM values indicate magnetic fields directed towards the observer and 
negative ones correspond to those pointing away. 

2 Determining the magnetic strength from the synchrotron intensity requires information about 
the number density of the cosmic-ray electrons. The latter can be obtained via X-ray emission, by 
inverse-Compton scattering, or through 7-ray bremsstrahlung. When such data is unavailable, an 
assumption must be made about the relation between cosmic-ray electrons and magnetic fields. This 
is usually the aforementioned principle of energy equipartition [49,50] 



Spiral galaxies observed in total radio emission appear very similar to those seen in the far- 
infrared. The equipartition magnetic strength in the arms can be up to 30/iG and shows a 
low degree of polarization. The latter indicates that the fields are randomly oriented there. 
On the other hand, synchrotron radio-emission from the inter-arm regions has a higher degree 
of polarization. This is due to stronger (10/xG - 15/xG) and more regular 5-fields, oriented 
parallel to the adjacent optical arm. The ordered fields form spiral patterns in almost every 
galaxy, even in ringed and flocculent galaxies. Therefore, the magnetic lines do not generally 
follow the gas flow (which is typically almost circular) and dynamo action is needed to explain 
the observed radial magnetic component. In galaxies with massive bars, however, the field lines 
appear to follow the gas flow. As the gas rotates faster than the bar pattern of the galaxy, 
a shock occurs in the cold gas. At the same time, the warm gas is only slightly compressed. 
Given that the observed magnetic compression in the spiral arms and the bars is also small, it 
seems that the ordered field is coupled to the warm diffuse gas and is strong enough to affect 
its flow [7]. 

Spiral dynamo modes can be identified from the pattern of polarization angles and Faraday 
rotation measures from multi- wavelength radio observations of galaxy disks [52], or from RM 
data of polarized background sources [53] . The disks of some spiral galaxies show large-scale 
RM patterns, but many galaxy disks posses no clear patterns of Faraday rotation. Faraday 
rotation in the direction of QSOs helps to determine the field pattern along the line of sight 
of an intervening galaxy [53,54]. Recently, high resolution spectra have unambiguously associ- 
ated quasars with strong Mgll absorption lines to large Faraday rotation measures. As Mgll 
absorption occurs in the haloes of normal galaxies lying along the line of sight to the quasars, 
this implies that organized strong B- fields are also present in high-redshift galaxies [8,9,10]. 

Magnetic fields have also been detected within clusters of galaxies, where X-ray observations 
have revealed the presence of hot gas [5]. There are several indications that favour the existence 
of cluster magnetic fields. In particular, galaxy clusters are known to have radio halos that 
trace the spatial distribution of the intra-cluster gas found in the X-ray observations. The 
radio signals are due to synchrotron emission from relativistic electrons spiralling along the 
field lines. In addition, there have been reports of Faraday rotation measurements of linearly 
polarized emissions crossing the intracluster medium. The first detection of a cluster magnetic 
field was made in the Coma cluster [55]. The Very Large Array (VLA) was used to compare 
Faraday rotation measures of radio sources within and directly behind the Coma cluster with 
radiation not crossing the cluster. Since then, there have been more analogous detections. It 
turns out that the observed cluster-field strengths vary slightly with the type of cluster. In 
particular, the magnetic field strength depends on whether we are dealing with cooling flow or 
non-cooling flow clusters. Faraday observations indicate turbulent field strengths of /iG-order 
in non-cooling flow clusters, such as the Coma. For cooling flow clusters, like the Hydra for 
example, the i?-fields are of the order of a few 10 /iG [56]. In fact, the cool core region of the 
Hydra A cluster is associated with a magnetic field of 7/xG with correlation length of 3 Kpc. 
Non-cooling flow clusters like the Coma, on the other hand, have weaker fields of the order 
of 3/iG but with larger correlation lengths (between 10 Kpc and 30 Kpc) [57]. In general, the 
magnetic structure is not homogeneous but is rather patchy on small scales (5 Kpc - 20 Kpc), 
indicating the presence of tangled magnetic fields [56]. An alternative way of determining the 



strength of cluster 5-flelds is to compare the radio synchrotron emission with inverse Compton 
X-ray emission [5]. The former comes from spiraling electrons along the cluster magnetic field. 
The latter is mainly due to CMB photons being upwardly scattered by the relativistic electrons 
of the intracluster gas. 



In view of the accumulating observational evidence for magnetic presence on all scales up to 
that of a galaxy cluster, the idea of a truly cosmological origin for cosmic magnetism gains 
ground. The potential detection of such primordial 5-fields in the intergalactic medium may 
also change our understanding of the way structure formation has progressed. Note that an 
intergalactic magnetic field ordered on very large scales would pick out a preferred direction, 
which should then manifest itself in Faraday rotation measurements from distant radio sources. 
This puts an upper limit on any cosmological intergalactic magnetic field of -Bigm i$ 1CT 11 [58]. 
Assuming that such a field has a characteristic scale, galaxy rotation measures suggest a size of 
1 Mpc and an upper limit of the order of 1 nG [44]. Indications of intergalactic magnetic fields 
have come from observations of radio-galaxy groupings near the Coma cluster, suggesting the 
presence of 5-fields with strengths between 0.2 /iG and 0.4 /iG and a coherence scale close to 
4 Mpc [59]. There is also evidence for an intergalactic magnetic field around 0.3 /zG on scales 
of the order of 500 Kpc, from excess rotation measures towards the Hercules and the Perseus- 
Pisces superclusters [60,61]. In addition, intergalactic 5-fields close to 30 nG and spanning 
scales of approximately lh^ 1 Mpc were recently suggested after cross-correlating the galaxy 
density field, obtained from the 6th Data Release of the Sloan Digital Sky Survey, with a large 
sample of Faraday rotation measures supplied by the NRAO-VLA Sky Survey [62] . 



Additional reports of intergalactic magnetic fields have appeared within the last year, using 
techniques that exploit the magnetic effects on the highly energetic photons emitted by distant 
sources (e.g. see [45,46,47,48]). More specifically, TeV-energy photons from a distant AGNs 
interact with the low frequency photons of the extragalactic background and lead to electron- 
positron pair creation. These produce (GeV-level) 7-rays through the inverse Compton scat- 
tering of the CMB photons. Observation wise, the key point is that a magnetic presence, 
even a very weak one, can affect the profile of the the resulting 7-ray spectra. For instance, 
the .B-field can cause the formation of an extended halo around the 7-ray images of distant 
AGNs. Such halos were first reported by Ando and Kusenko, using combined data from the 
Atmospheric Cherenkov Telescopes and Fermi Gamma- Ray Space Telescope [12]. Subsequent, 
complementary analysis indicated the presence of an intergalactic magnetic field with strength 
between 10~ 17 and 10~ 14 Gauss [15]. In addition to halo formation, the -B-field can also reduce 
the observed flux of the secondary GeV-photons by deflecting them into larger solid angles. 
Using observations of the Fermi/Large Area Telescope and assuming that the original TeV- 
photons were strongly beamed, a lower limit of ~ 10~ 15 G was imposed on the intergalactic 
magnetic field [11,14]. A similar lower limit of ~ 10~ 16 G was also obtained assuming that the 
blazar source radiated isotropically [13]. 
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2.2 Limits on primordial magnetic fields 



Any primordial magnetic field must comply with a number of astrophysical constraints, the 
most significant of which come from Big Bang Nucleosynthesis (BBN) and the isotropy of the 
CMB. The latter probes B- fields with coherence scales larger than the particle horizon during 
nucleosynthesis, while the BBN limits apply in principle to all scales. 



2.2.1 Nucleosynthesis limits 

The main effects of a magnetic presence on the output of primordial nucleosynthesis are related 
with: (a) the proton-to-neutron conversion ratio; (b) the expansion and cooling of the universe; 
and (c) the electron thermodynamics. Here we will only provide a very brief summary of these 
effects. For a detailed review, the reader is referred to [20]. 

(a) In the early universe, the weak interaction is responsible for maintaining chemical equilib- 
rium between protons and neutrons. The main effect of a strong magnetic presence at the time 
of nucleosynthesis is to enhance the conversion rate of neutrons into protons. As a result, the 
neutron-to-proton ratio would freeze-out at a lower temperature. This in turn would lead to 
a less efficient production of 4 He and of heavier elements [63,64]. In fact, the magnetic effect 
would be catastrophic if B ^> rn^/e ~ 10 17 G at the time of nucleosynthesis. 

(b) The temperature at which the proton-to-neutron ratio freezes-out is determined by the bal- 
ance between the timescale of the weak interaction and the expansion rate of the universe [65] . 
Equilibrium is attained when r n _» p ~ H, where T n ^ p is the cross-section of the interaction 
and H is the Hubble constant at the time. The latter is proportional to the total energy den- 
sity of the universe, where the -B-field contributes as well. Thus, a strong magnetic presence 
will increase the value of the Hubble parameter. This would cause an earlier freeze-out of the 
proton-to- neutron ratio and result into larger residual amounts of 4 He [66,67]. 

(c) The magnetic presence will also change the phase-space volume of electrons and positrons, 
since their momentum component normal to the B- field will become discrete (Landau levels). 
Therefore, the energy density, the number density and the pressure of the electron gas augment, 
relative to their magnetic- free values [64]. The rise happens at the expense of the background 
photons, which transfer energy to the lowest Landau level. This delays the electron-positron 
annihilation, which in turn increases the photon-to-baryon ratio and finally leads to lower 3 He 
and D abundances [68]. 

All of the above need to be accounted for when calculating the BBN limits on primordial 
magnetic fields. This is done by means of numerical methods, which seem to conclude that 
the main magnetic effect on the light-element abundances comes from the field's contribution 
to the expansion rate of the universe (i.e. case (b)). The overall constraint on the magnetic 
strength is B < 10 11 G at the time of nucleosynthesis, which (roughly) translates to B < 
7 x 1CT 7 G at the time of galaxy formation 
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2.2.2 Cosmic microwave limits 

Observations of the CMB temperature anisotropies and polarization provide valuable tools to 
constrain cosmological models. As such, they also play an important role in the diagnostic of 
early universe magnetic fields. 

In comparison to the data of the angular power spectra of polarization, Cf E , and temperature- 
polarization, Cj E , the temperature angular power spectrum, Cf T , is known at higher preci- 
sion. For example, the 7- year WMAP ( Wilkinsion Microwave Anisotropy Probe) power spec- 
trum is limited only by cosmic variance up to £ f« 548 [69]. Moreover, on smaller-scales, 
observations from CBI (Cosmic Background Imager) [70] and VSA (Very Small Array) [71], 
ACBAR (Arcminute Cosmology Bolometer Array Receiver) [72] and the forthcoming SPT 
(South Pole Telescope) [73] missions will determine the Cf T to even higher accuracies. The 
PLANCK satellite is expected to extend the region limited only by cosmic variance to £ ~ 1500. 

At the moment, the high isotropy of the Cosmic Microwave Background appears to exclude 
homogeneous cosmological magnetic fields much stronger than ~ 10~ 9 G [74]. A similar limit is 
found for stochastic magnetic fields as well [75,76,77,78]. It has been shown that the tempera- 
ture angular power spectrum, Cf T , from magnetically-induced vector and scalar perturbations 
increases slightly across all angles. The extra pressure that the -B-field adds into the system 
can change the position and the magnitude of the acoustic peaks, thus producing a poten- 
tially observable effect [79,80]. The presence of small-scale magnetic fields appears to leave 
undamped features on small angular scales and may also lead to distinctive polarisation struc- 
tures [81,82]. In addition, large-scale primordial fields could be related to the low-quadrupole 
moment problem [83]. Nevertheless, the magnetic signal remains subdominant to that from 
standard scalar perturbations until around £ ~ 2000, depending on the field strength and 
spectral index [84,85,86,87,88]. Magnetic fields also source tensor modes, which however are of 
relatively low-amplitude. The signal is similar to that of inflationary gravitational wave, but 
probably weaker in strength [89]. Therefore, the direct magnetic impact onto the CMB (TT) 
correlation does not generally provide an ideal probe of primordial magnetism. However, the 
CBI mission observed a weak increase of power on small scales, as compared to the concor- 
dance model [71]. Provided this is real and not a statistical or systematic artifact, it could 
be partly explained by the presence of a cosmological .B-field [84,90]. Nucleosynthesis bounds, 
however, imply that a primordial field is unlikely to account for all the increase. 

Besides contributing to the CMB temperature fluctuations, a primordial magnetic field also 
produces E-mode polarisation that can significantly change the angular power spectrum of 
the standard ACDM model [85]. However, the polarisation limits are not as strong as those 
coming from the temperature anisotropy. Due to the presence of both vector and tensor 
perturbations, the magnetic field also leads to 5-mode polarisation. Moreover, 5-modes are 
also induced in the scalar sector by Faraday rotation if a magnetic field is present at de- 
coupling [91,92,93,94,95,96]. Taking into account that in the standard picture 5-modes are 
produced only by lensed £'-modes and by inflationary gravitational waves, in principle, the 
observation of a distinct S-mode power spectrum would be the clearest indication of a pri- 
mordial magnetic field. However, the CMB polarisation maps are poorly known, compared 
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to the temperature ones. While we currently possess a power spectrum Cj E , this is by no 
means cosmic-variance limited on any scale. The observations of the 5-modes yield bounds 
consistent with zero [97,98]. These are on relatively small scales, directly observing the region 
at which magnetic effects may come to dominate. Nevertheless, we are far from the required 
accuracy, particularly for the 5-modes. Given the limitations of the power spectra, the non- 
Gaussianity of the temperature map is a reasonable place to look for further constraints on 
primordial magnetic fields [99,100,101,102,103,104,105]. Although up to now the observations 
are entirely consistent with Gaussian initial conditions, there are non-Gaussian features in the 
WMAP maps [106]. Also, the number of non-Gaussian features could well increase with the 
next generation of CMB experiments. 



2.2.3 Limits from gravitational waves 

A strong limit on stochastic magnetic fields produced before nucleosynthesis has been derived 
in [107]. The anisotropic stress of the magnetic field acts as a source term in the evolution 
equation of gravity waves. This causes the conversion of magnetic field energy into gravity 
waves above a certain critical value of the magnetic field strength. In particular, the field 
strength smoothed over a scale A of magnetic fields generated during inflation must be smaller 
than B\ ~ 10~ 20 G for spectral indices n-Q > —2, where ne = —3 corresponds to a scale invari- 
ant magnetic field energy spectrum. If the magnetic field is produced by a causal mechanism, 
for example during the electroweak phase transition, n^ > 2, its strength has to be below 
10~ 27 G in order not to loose all its energy density to gravitational waves. 

The magnetic-strength limits asserted in [107] are the strongest reported in the literature, 
far more restrictive than those coming from nucleosynthesis or the CMB. However, analogous 
studies of magnetically produced gravity waves have reached different conclusions. It has been 
claimed, in particular, that the limits on cosmological magnetic fields set by the latest LIGO S5 
data lie close to those obtained by BBN and the CMB [108]. 

Finally, we should also note the possibility of constraining primordial -B-fields using the ioni- 
sation history of the post-recombination universe and, in particular, the observed re-ionisation 
depth. Thus, based on the 5-year WMAP data, upper limits of nGauss order have been re- 
ported in the literature [109]. 



13 



3 Magnetic fields in cosmology 

3.1 Relativistic magnetised cosmologies 



Although the study of large-scale magnetic fields goes a long way back into the past, the first 
systematic attempts to incorporated magnetism in cosmology appeared in the late 60s and 
the early 70s [110,111,112,113]. Next, we will provide the basic background for the relativistic 
study of cosmological 5-fields. For the details and a recent review the reader is referred 
to [114,115,116,117,118,119]. 



3.1.1 The gravitational field 

In the geometrical framework of general relativity, gravity is a manifestation of the non- 
Euclidean geometry of the spacetime. The gravitational field is therefore described by the 
Riemann curvature tensor (R a bcd), which satisfies the Ricci identities 

2V [a V b] v c = R abcd v d , (3.1.1) 

applied here to an arbitrary vector v a (with V a representing the familiar covariant derivative 
operator). The Riemann tensor also assumes the invariant decomposition 

1 1 

Rabcd = 77 [gacRbd + QbdRac ~ QbcRad + g a dRbc) ~ 77 R (g aC gbd ~ 9ad9bc) + C a bcd , (3.1.2) 

2 6 

and obeys the symmetries R a bcd = Rcdab and R a bcd = R\ab]{cd}- Note that g a b is the spacetime 
metric, R ab = R c ac b is the Ricci tensor, R = R a a is the associated Ricci scalar and C abc( i 
is the Weyl (or conformal curvature) tensorjjj The Ricci component of the Riemann tensor 
determines the local gravitational field through the Einstein field equations 

Rab - - Rg a b = KT ab - Ag ab , (3.1.3) 

where T ab is the energy-momentum tensor of all the matter fields involved, k = 8ttG and A the 
cosmological constant o The Weyl tensor, on the other hand, has to do with the long-range 
component of the gravitational field (i.e. tidal forces and gravity waves), shares the same 
symmetries with R abc d and it is also trace- free (i.e. C c acb = 0). 



3 Unless stated otherwise, we consider a general 4-dimensional (pseudo) Riemannian spacetime with 
a Lorentzian metric of signature (— , +, +, +). Also, throughout this review, Latin indices take values 
between and 3, while their Greek counterparts vary from 1 to 3. 

4 We use the Heaviside-Lorentz units for the electromagnetic field in this section. Natural units, 
with c = 1 = ks = h and energy as the fundamental dimension, are used throughout this review. 
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3.1.2 Kinematics 

We introduce a family of observers with worldlines tangent to the timelike 4-velocity field u a 
(i.e. u a u a = —1). These are the fundamental observers that define the direction of time. Then, 
the tensor h ab = g a b + UaU b projects orthogonal to u a and into the observers' instantaneous 
3-dimensional rest-spacejjj Together, u a and h ab introduce an 1+3 'threading' of the spacetime 
into time and space, which decomposes physical quantities, operators and equations into their 
irreducible timelike and spacelike parts (see [120,121] for further details). 

For example, splitting the covariant derivative of the observers' 4-velocity, leads to the irre- 
ducible kinematic variables of the motion. In particular, we arrive at 

V b u a = -Qh ab + a ab + u ab - A a u b , (3.1.4) 

where = V a u a = D a u a is the volume scalar, a ab = D( b u a ) is the shear tensor, u ab = D[ b u a ] 
is the vorticity tensors and A a = ii a = u b V b u a is the 4-acceleration vectorial The volume 
scalar describes changes in the average separation between neighbouring observes. When G 
is positive this separation increases, implying that the associated fluid element expands. In 
the opposite case we have contraction. The volume scalar also defines a representative length 
scale (a) by means of a/a = 6/3. In cosmological studies, a is commonly referred to as the 
'scale factor'. We use the shear to monitor changes in the shape of the moving fluid under 
constant volume, while the vorticity traces its rotational behaviour. Note that we can replace 
the vorticity tensor with the vorticity vector u a = e abc co bc /2, where e abc represents the 3-D 
Levi-Civita tensor. Finally, the 4-acceleration reflects the presence of non-gravitational forces 
and vanishes when the observers worldlines are timelike geodesies. 

The time evolution of the volume scalar, the vorticity vector and the shear tensor is determined 
by a set of three propagation equations, supplemented by an equal number of constraints. Both 
sets are onbtained after applying the Ricci identities (see (3.1.1) in § 3.1.1) to the fundamental 
4-velocity field [120,121]. 



3.1.3 Matter fields 

Analogous decompositions apply to the rest of the kinematical and dynamical variables. Thus, 
relative to the w a -frame, the energy-momentum tensor of a general (imperfect) fluid splits as 

T a ( ™ } = pu a u b + ph ab + 2q^ a u b) + -n ab . (3.1.5) 



By construction, h ab is a symmetric spacelike tensor, with h a a = 3 and h ab h c = h ac . The projector 
coincides with the metric of the observers' 3-dimensional space in non-rotating spacetimes. 
6 Overdots indicate (proper) time derivatives along the u a -field, while the gradient D a = h a b V b de- 
fines the 3-dimensional covariant derivative operator. Round brackets denote symmetrisation, square 
antisymmetrisation and angled ones indicate the symmetric and traceless part of projected tensors 
and vectors. For example, a ab = D^-u^ = D( b u a ) — (D c u c /3)h ab and Ei a \ = h a b E b - see Eq. (3.1.8). 
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Here, p = T ab u a u b represents the energy density, p = T a bh ab /3 the isotropic pressure, q a = 
—h a b T bc u c the total energy flux and 7f ab = h^ a c h b ) d T cd the anisotropic pressure of the matter, 
as measured by the fundamental observers [120,121]. When dealing with a perfect fluid, both 
q a and n ab vanish and the above reduces to 

Ta^ = P U a U b + Phab • (3.1.6) 

The remaining degrees of freedom are determined by the equation of state, which for a 
barotropic medium takes the simple p = p(p) form. 



3.1.4 Electromagnetic fields 

Magnetic and electromagnetic fields introduce new features to any cosmological model through 
their energy density and pressure contributions and due to their generically anisotropic nature. 
The Maxwell field is invariantly described by the antisymmetric Faraday tensor. Relative to 
the fundamental observers introduced in § 3.1.2, the latter decomposes as 

F a b = 2u[ a E b ] + e abc B c , (3.1.7) 

where E a = F ab u b and B a = e a b c F bc /2 are respectively the electric and magnetic components. 

The inherit anisotropy of the electromagnetic field is reflected in the form of its energy- 
momentum tensor. The latter has the invariant form T a e h = —F ac F c b + (F ca -F cd /4:)g ab , which 
relative to the M a -frame recasts to 

T a ( b em) = \ (B 2 + E 2 ) u a u b + i (B 2 + E 2 ) h ab + 2V {a u b) + U ab . (3.1.8) 

with E 2 = E a E a , B 2 = B a B a , V a = e abc E b B c and 

1 



n a & = —E( a E h ) — -B(a-Bfe) = - [E + B J h a b — E a E b — B a B b , (3.1.9) 

Comparing the above to (3.1.5) in § 3.1.3, we conclude that the Maxwell field corresponds to 
an imperfect fluid with energy density (E 2 + B 2 )/2, isotropic pressure (E 2 + B 2 )/6, energy 
flux given by the Poynting vector V a and anisotropic stresses represented by the symmetric 
and trace-free Il a fe-tensor. This fluid-like description of the Maxwell field has been proved 
particularly helpful in many applications [114,115,116,117,118,119]. 



3.1.5 Conservation laws 

In the case of charged matter, the total energy-momentum tensor is T ab = T a ™ + T a l , , 
with the individual components given by (3.1.5) and (3.1.8) respectively. The total stress- 
tensor satisfies the familiar conservation law V b T ab = 0, while the electromagnetic stress- 
energy tensor obeys the additional constraint V fe T a ^ m ^ = —F ab J b , with the right-hand side 



16 



representing the Lorentz 4-force. Combining the two, we obtain the conservation laws of the 
total energy and momentum. These are given by the continuity equation 

p = _0(p +p )- B a q a - 2A a q a - a ab n ab + E a J a , (3.1.10) 

and by the Navier-Stokes equation 

(p + p)A a = -D a p - q {a) - - Qq a - (a ab + u ab )q b - B b n ab - n ab A b 

+fiE a -e abc B b J c , (3.1.11) 

respectively [119]. Note that // = —J a u a is the electric charge density and J a = h a b J b is the 
associated 3-current, so that J a = jiu a + J a - 

An additional conservation law is that of the 4-current, which satisfies the invariant constraint 
V a J a = 0. The latter translates into the conservation law for the charge density, given by [119] 

p = -Opt - D a J a - A a J a . (3.1.12) 

3.2 Evolution of the electromagnetic field 



The vector nature of the electromagnetic components and the geometrical approach to gravity 
that general relativity introduces, mean that the Maxwell field is the only known energy source 
that couples directly to the spacetime curvature through the Ricci identities as well as the 
Einstein Field Equations. Both sets are therefore necessary for the full relativistic treatment 
of electromagnetic fields. 



3.2.1 Maxwell's equations 

We monitor the evolution of the electromagnetic field using Maxwell's formulae. In their 
invariant form these read 

V [c F ab] =0 and V b F ab = J a , (3.2.1) 

with the first manifesting the existence of a 4-potential. Relative to the -u a -frame, the above 
set splits into two pairs of propagation and constraint equations. The former consists of 



2 
3 



E {a) = --BE a + (a ab + u ab ) E b + e abc A b B c + curl^ - J a , (3.2.2) 



and 



B {a) = -- QB a + (a ab + u ab ) B b - e abc A b E c - cm\E a , (3.2.3) 
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which may be seen as the 1+3 covariant analogues of the Ampere and the Faraday laws 
respectively. The constrains, on the other hand, read 

D a E a = fi - 2u a B a and D a B a = 2co a E a . (3.2.4) 

These provide the 1+3 forms of Coulomb's and Gauss' laws respectively. Note that Eqs. (3.2.2)- 
(3.2.4) contain relative motion effects, in addition to the standard 'curl' and 'divergence' terms 
of their more traditional versions. This is an essentially built-in property of the 1+3 formalism, 
which should be always kept in mind when applying expressions like the above. 



3.2.2 The electromagnetic wave equations 

Maxwell's equations also provide the wave equation of the electromagnetic tensor. This can 
be obtained by applying the Ricci identities to the Faraday tensor and takes the invariant 
form (e.g. see [118,122]) 

V 2 F ab = -2R abcd F cd - 2R [a c F b]c - 2V [a J b] , (3.2.5) 

where V 2 = V a V a is the d'Alembertian and R abc di Rab are the Riemann and Ricci tensors 
respectively (see § 3.1.1). The above results from the vector nature of the electromagnetic com- 
ponents and from the geometrical interpretation of gravity that general relativity advocates. 
The two guarantee that the Maxwell field is the only known source of energy that couples 
directly to gravity through both the Einstein equations and the Ricci identities. 

Expression (3.2.5), which clearly shows how spacetime curvature drives electromagnetic dis- 
turbances, can also provide the wave-equations of the individual components of the Maxwell 
field. Alternatively, one may obtain these relations using the set (3.2.2)-(3.2.4) together with 
the Ricci identities. Assuming matter in the form of a single perfect fluid, we obtain [118] 



E {a) - B 2 E a =(a ab - e abc u c - - Qh ab ) E b + - k P (\ + 3w)E a 



+ - (o-ab + £ab c u c ~ - ®h ab j QE a - a {a c o- b)c E b + e abc E b a cd u d 

+^ (V - jj c/) E a + i u {a u b) E b + A b A b E a - | e abc A b cm\E c + D {a E b) A b 

+ jj e abc B b D c e + s abc B d D b a cd + B {a cu b) B b + ^ e abc B b cur\tu c + 2D {a B b) cu b 

-2e abc a b d D^ c B d > + e abc A b B c + 7 - A b uo b B a + ^ B b uj b A a - 3A b B b co a 

+3e abc A <y c B d — 1Z ab E — E ab E + H ab B 

+-fiA a -D af i-j a -QJ a , (3.2.6) 



for the electric field, and 
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/ 5 \ • 1 

B {a) - D 2 B a = [a ab - e abc co c - - Bh ab j B b + - «p(l + 3w)B a 

+ 3 \ ab + £abcUjC ~ 3 ® hab ) QBa " a ^ C(7h ) cBh + £abcB b (i cd uj d 

+^ [a 2 - jj c/) B a + 1 w (a u; 6>J B 6 + A%£? a - ^ e abc A b cni\B c + D (aJ B 6> A 6 

-H £ a6c E b D c 6 - e abc E d D b a cd - D {a co h) E b - ^ e ahc E b cm\co c - 2D {a E b) co b 

+2e abc a b d rt c E d ^ - e abc A b E c - l A b u b E a - ± E b u b A a + 3A b E b u a 

— 3e abc A a c E d — TZabB — E ab B — H ab E 

-- fiu a + cmlJa + 2e abc A b J c , (3.2.7) 

for its magnetic counterpart. Here, TZ ab is the Ricci tensor of the observer's 3-D rest-space, 
while E ab and H ab are the electric and the magnetic parts of the Weyl tensor respectively 
(see [121] for details). Note the curvature terms in the right-hand side of the above, which 
show explicitly how the different parts of the spacetime geometry can affect the propagation 
of electromagnetic signals. 



3.3 Cosmological magnetohydrodynamics 



The evolution and the implications of cosmological electromagnetic fields depend on a number 
of parameters. A decisive factor is the electrical properties of the medium that fills the universe, 
which in turn are determined by the specific form of Ohm's law. 



3.3.1 Ohm's law 

The literature contains various expressions of Ohm's law, which in principle is the propagation 
equation of the electric 3-current (e.g see [123,124,125]). For a single fluid at the limit of 
resistive MHD, Ohm's law takes the simple form [126,127] 

J a = sE a , (3.3.1) 

with q representing the electric conductivity of the matter. In highly conducting environments, 
q — > oo and the electric field vanishes. This is the familiar ideal-MHD approximation, where 
the electric currents keep the magnetic field frozen-in with the charged fluid. At the opposite 
end, namely when the conductivity is very low, q — > 0. Then, the 3-currents vanish despite 
the presence of nonzero electric fields. 
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3.3.2 The ideal MHD limit 

At the ideal MHD limit Maxwell's formulae reduce to one propagation equation and three 
constraints. The former comes from (3.2.3) and is the familiar magnetic induction equation 

B {a) = --QB a + {(x ab + uj ab )B b . (3.3.2) 

The constraints, on the other hand, are obtained from Eqs. (3.2.2) and (3.2.4). In particular, 
eliminating the electric field form these relations, we arrive at 

J a = e abc A b B c + cm\B a , n = 2uu a B a and D a B a = , (3.3.3) 

respectively. Following (3.3.3b), the inner product u a B a corresponds to an effective charge 
density, triggered by the relative motion of the 5-field. 

In the absence of the electric field, the electromagnetic energy-momentum tensor simplifies as 
well. To be precise, expression (3.1.8) reduces to 

T a ( f } = I B 2 u a u h + \ B 2 h ab + U ab , (3.3.4) 

2 6 

with U ab = —B( a B b ). This means that the magnetic field can be seen as an imperfect fluid 
with energy density given by B 2 /2, isotropic pressure equal to B 2 /6, zero energy flux and 
anisotropic stresses given by 

n afe = i B 2 h ab - B a B b . (3.3.5) 

At the MHD limit, the matter energy and that of the residual magnetic field are conserved 
separately, with the induction equation (see formula (3.3.2) above) providing the conservation 
law of the magnetic energy. At the same time, the momentum conservation law reads 

2 \ 4 

p + p + -B 2 ) A a = -B a p - q {a) - - Qq a - (cr ab + u ab )q b - B b n ab - n ab A b 

-e abc B b cm\B c - U ab A b . (3.3.6) 

When dealing with a perfect fluid with zero pressure, we may set p = = q a = ir ab . Then, 
starting from (3.3.6), one can show that A a B a = and recast the latter into the form 

p + B 2 ) A a = -e abc B b cm\B c = ~ D a B 2 + B b D b B a . (3.3.7) 

where in the right-hand side we have two expressions for the Lorentz force. Note that the first 
term in the last equality is due to the field's pressure, while the second carries the effects of the 
magnetic tension. The former reflects the tendency of the field lines to push each other apart 
and the latter their elasticity and tendency to remain 'straight'. As we will explain below, the 
majority of studies analysing the magnetic effects on structure formation do not account for 
the tension contribution to the Lorentz force. 
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3.4 Magnetism and structure formation 



Despite the widespread presence of large-scale magnetic fields in the universe, their role and 
implications during structure formation are still not well understood. Here, we will briefly 
summarise the way 5-fields could have altered the linear and the mildly non-linear stages of 
galaxy formation and also how the latter might have backreacted on the magnetic evolution. 



3.4-1 The linear regime 

Scenarios of magnetised structure formation typically work within the ideal-MHD approxima- 
tion, looking at the effects of the magnetic Lorentz force on density inhomogeneities. The bulk 
of the available inhomogeneous treatments are Newtonian, with the relativistic approaches 
being a relatively recent addition to the literature. Although the role of the magnetism as a 
source of density and vorticity perturbations was established early on [128,129,130], the com- 
plicated action of the -B-field did not allow for analytic solutions (with the exception of [131] 
for the case of dust). Full solutions were provided by means of covariant techniques, which 
considerably simplified the mathematics [114,115,116,117,118,119]. 

Magnetic fields generate and affect all three types of density inhomogeneities, namely scalar, 
vector and (trace-free) tensor inhomogeneities. The former are those commonly referred to as 
density perturbations and represent overdensities or underdensities in the matter distribution. 
Vector inhomogeneities describe rotational (i.e. vortex-like) density perturbations. Finally, 
tensor-type inhomogeneities correspond to shape distortions. Following [114,115,116,117,118,119], 
the scalar 

A = — D 2 p, (3.4.1) 

P 
describes linear perturbations in the density (p) of the matter and corresponds to the more 

familiar density contrast Sp/p. Note that positive values for A indicate overdensitites and 

negative ones underdensities. In a perturbed, weakly magnetised and spatially flat Friedmann- 

Robertson- Walker (FRW) universe, the above defined scalar evolves according td_j 

A = 3wHA - (1 + w)Z + | c 2 (l + w)HB , (3.4.2) 

where Z = a 2 D 2 G and B = (a 2 /B 2 )D 2 B 2 . The first of the last two variables describes linear 
inhomogeneities in the smooth Hubble expansion and the second represents perturbations in 
the magnetic energy density. Then, to linear order, 

Z = -2HZ - - pA + - c 2 (l + w)pB ^— D 2 A - - c 2 D 2 £ (3.4.3) 

2 F 4 aV JH 1 + w 2 a v; 

and 

B = — — — -A, (3.4.4) 

3(1 + w) ' K J 



7 Basic aspects of the FRW dynamics and evolution are discussed in § 4.1.1 and § 4.1.2. 
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respectively. In the above w = p/p is the (constant) barotropic index of the matter, H = a/a 
is the background Hubble parameter, c 2 = p/p is the square of the adiabatic sound speed and 
c 2 = B 2 / p{l + w) is that of the Alfven speed. We have also assumed that B 2 <C p, given the 
relative weakness of the magnetic field. 

Expression (3.4.2) verifies that 5-fields are generic sources of linear density perturbations. 
Indeed, even when A and Z are zero initially, A will generally take nonzero values due to 
the magnetic presence. Also, Eq. (3.4.4) shows that linear perturbations in the magnetic 
energy density evolve in tune with those in the density of the matter (i.e. B oc A). This 
means that a 5-field residing in an overdense region of an Einstein-de Sitter universe will 
grow by approximately two to three orders of magnitude (see solution (3.4.7) below). Note 
that the aforementioned growth occurs during the linear regime of structure formation and is 
independent of the (nonlinear) increase in the field's strength due to the adiabatic compression 
of a protogalactic cloud (see § 3.4.2 for more details). Finally, we should emphasise that only 
the pressure part of the Lorentz force contributes to Eqs. (3.4.2) and (3.4.3). To account for the 
tension effects at the linear level, one needs to allow for FRW backgrounds with non- Euclidean 
spatial geometrypU 

The system (3.4.2)-(3.4.4) has analytical solutions in the radiation and the dust eras. Before 
equipartition, when w = 1/3 = c 2 , H = l/2t, p = 3/4t 2 and c 2 = 3B 2 /4p = constant, large- 
scale magnetised density perturbations obey the power-law solution. In particular, keeping 
only the dominant growing and decaying modes one arrives at [116,119] 



A = c 1 r 1 / 2+10c '/ 9 + c 2 t 1 - 



■4^/9 



(3.4.5) 



In the absence of the 5-field, we recover the standard growing and decaying modes of A oc t 
and A oc £ -1 / 2 respectively. So, the magnetic presence has reduced the growth rate of the 
density contrast by 4c 2 /9. Also, since B oc A - see Eq. (3.4.4), the above describes the linear 
evolution of the magnetic energy- density perturbations as well. 

Well inside the horizon we can no longer ignore the role of the pressure gradients. There, the 
fc-mode oscillates like a magneto-sonic wave with 



A 



(k) 



oc sin 




(3.4.6) 



where A& = a/k is the perturbed scale and Xh = 1/H the Hubble horizon [116,119]. Here, the 
magnetic pressure increases the effective sound speed and therefore the oscillation frequency. 
The former makes the Jeans length larger than in non-magnetised models. The latter brings 
the peaks of short-wavelength oscillations in the radiation density closer, leaving a potentially 
observable signature in the CMB spectrum [79]. 



8 The Lorentz force splits according to e a b c B cur\B c = (l/2)D a B 2 — B Dt,B a , with the former term 
carrying the effects of the magnetic pressure and the latter those of the field's tension (see Eq. (3.3.7) 
in § 3.3.2). Given that T) a B a = at the ideal-MHD limit, only the pressure part affects linear density 
perturbations, unless FRW backgrounds with nonzero spatial curvature are involved. 
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When dust dominates, w = = c 2 , H = 2/3£, p = 4/3£ 2 and c 2 = B 2 /p oc t 2 / 3 . Then, on 
superhorizon scales, the main growing and decaying modes of the density contrast are [114,119] 

A = &t ai + C 2 t a2 , (3.4.7) 



with a 1;2 = — [1 ±5yl — (32/75)(c a Xh/X^qI/Q. In the absence of the 5-field we recover again 
the standard solution with oc\ = 2/3 and a 2 = —1. Thus, as with the radiation era before, the 
magnetic presence slows down the growth rate of density perturbations. In addition, the field's 
pressure leads to a magnetically induced Jeans length, below which density perturbations 
cannot grow. As a fraction of the Hubble radius, this purely magnetic Jeans scale is 

Xj ~ c a \ H . (3.4.8) 

Setting B ~ 10~ 9 G, which is the maximum homogeneous field allowed by the CMB [74], 
we find that Xj ~ 10 Kpc. Alternative, magnetic fields close to 1CT 7 G, like those found in 
galaxies and galaxy clusters, give Xj ~ 1 Mpc. The latter lies intriguingly close to the size of 
a cluster of galaxies. 

Overall, the magnetic effect on density perturbations is rather negative. Although 5-fields 
generate this type of distortions, they do not help them to grow. Instead, the magnetic presence 
either suppresses the growth rate of density perturbations, or increases the effective Jeans 
length and therefore the domain where these inhomogeneities cannot grow. The negative role 
of the i?-field, which was also observed in the Newtonian treatment of [131], reflects the fact 
that only the pressure part of the Lorentz force has been incorporated into the equations. 
When the tension component (i.e. the elasticity of the field lines) is also accounted for, the 
overall magnetic effect can change and in some cases it could even reverse [117]. 

Magnetic fields also induce and affect rotational, vortex-like, density inhomogeneities. To linear 
order, these are described by the vector W a = — (a 2 /2p)e a6c D fe D c p. Then, on an spatially flat 
FRW background, 

W a = -4HW a - l - pW a + i c 2 D 2 W a , (3.4.9) 

after matter-radiation equality [116,119]. Defining A a = c a A# as the 'Alfven horizon', we may 
write the associated solution in the form 

W (fc )=dt ai +C 2 t a2 , (3.4.10) 



where ai i2 = — [5±y 1 — (48/9)(A a /Afc)g]/6. On scales well outside the Alfven horizon, namely 
for A a ^C Afc, the perturbed mode decays as W oc £ -2 / 3 . This rate is considerably slower than 
W oc t _1 , the decay rate associated with magnetic-free dust cosmologies. In other words, the 
-B-field has reduced the standard depletion rate of the vortex mode. An analogous effect is 
also observed on co a , namely on the vorticity proper [116,119]. Hence, magnetised cosmologies 
rotate faster than their magnetic-free counterparts. In contrast to density perturbations, the 
field seems to favour the presence of vorticity. This qualitative difference should probably be 
attributed to the fact that the tension part of the Lorentz force also contributes to Eq. (3.4.9). 
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In addition to scalar and vector perturbations, magnetic fields also generate and affect tensor- 
type inhomogeneities that describe shape-distortions in the density distribution [116]. An 
initially spherically symmetric inhomogeneity, for example, will change shape due to the mag- 
netically induced anisotropy. All these are the effects of the Lorentz force. Even when the 
latter is removed from the system, however, the 5-field remains active. Due to its energy den- 
sity and anisotropic nature, for example, magnetism affects both the local and the long-range 
gravitational field. The anisotropic magnetic pressure, in particular, leads to shear distortions 
and subsequently to gravitational-wave production. Overall, magnetic fields are a very versa- 
tile source. They are also rather unique in nature, since 5-fields are the only known vector 
source of energy. An additional unique magnetic feature, which remains relatively unexplored, 
is its tension. When we add to all these the widespread presence of magnetic fields, it makes 
sense to say that no realistic structure formation scenario should a priori exclude them. 



3.4-2 Aspects of the nonlinear regime 

The evolution of large-scale magnetic fields during the nonlinear stages of structure formation 
has been addressed primarily by means of numerical methods. The reason is the high com- 
plexity of the nonlinear MHD equations, which makes analytical studies effectively impossible, 
unless certain simplifying assumptions are imposed. 

The simplest approximation is to assume spherically symmetric compression. Realistic col- 
lapse, however, is not isotropic. In fact, when a magnetic field is present, its generically 
anisotropic nature makes the need to go beyond spherical symmetry greater. Anisotropic 
contraction can be analytically studied within the Zeldovich approximation [132,133]. The 
latter is based on a simple ansatz, which extrapolates to the nonlinear regime a well known 
linear result. The assumption is that the irrotational and acceleration- free linear motion of 
the dust component, also holds during the early nonlinear stages of galaxy formation. This 
approximation allows for the analytical treatment of the nonlinear equations, the solution 
of which describe anisotropic (one dimensional) collapse and lead to the formation of the 
well-known Zeldovich 'pancakes'. 

Suppose that a magnetic field is frozen into a highly conductive protogalactic cloud that is 
falling into the (Newtonian) potential wells formed by the Cold Dark Matter (CMB) sector] 9 ! 
Relative to the physical coordinate system {r a }, the motion of the fluid velocity is u a = 
3Hr a + v a , where H = a/ a is the Hubble parameter of the unperturbed FRW background 
and v a is the peculiar velocity of the fluid (with a = 1,2,3). Then, the magnetic induction 
equation reads [134] 

B a = -2HB a --#B a + a aP B^ , (3.4.11) 

where overdots now indicate convective derivatives (i.e. ' = dt + vrdp). Also, d = d a v n and 
o"a/3 = d(pv a ) are the peculiar volume scalar and the peculiar shear tensor respectively! 10 1 The 



9 The Newtonian theory is a very good approximation, since the scales of interest are well below the 
curvature radius of the universe and we are dealing with non-relativistic matter. 

10 When dealing with purely baryonic collapse, the Zeldovich ansatz only holds during the early 
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former takes negative values (i.e. d < 0), since we are dealing with a protogalactic cloud that 
has started to 'turn around' and collapse. Note that the first term in the right-hand side of 
(3.4.11) represents the background expansion, the second is due to the peculiar contraction 
and the last reflects the anisotropy of the collapse. Introducing the rescaled magnetic field 
B a = a 2 B a , the above expression recasts into 

& a = -l#B a + a af) B f >, (3.4.12) 

with primes indicating differentiation with respect to the scale factor. Also d = aHd and 
a a p = aHa aj 3, where ■& = d a v a and a a p = d^v a ) (with v a = ax' a and v a = aHv a ). Relative to 
the shear eigenframe, a a p = (cr n , a 2 2, 033) and expression (3.4.12) splits into 

B[ = -\ m + a ll B 1 , B' 2 = -\ m + ~a 22 B 2 (3.4.13) 

and 

B' 3 = -^B 3 + a 33 B 3 . (3.4.14) 

This system describes the second-order evolution of a magnetic field, which is frozen in with 
the highly conductive matter of a collapsing protogalaxy, within the limits of the Zeldovich 
approximation. 

In order to solve the set of Eqs. (3.4.13), (3.4.14), we recall that in the absence of rotation 
and acceleration, the peculiar volume scalar is given by 

^ = 7tV + 7tV + 7tV' ( 3A15 ) 

1 + a\\ 1 + a\ 2 1 + a\ 3 
Similarly, the shear eigenvalues are 

°n = —^--U, Zn = -^—-\# (3.4.16) 

1 + a\i 3 1 + a\ 2 3 

and 

^ = TxV"^- ( 3A17 ) 

1 + a\ 3 3 

where Ai, A2 and A3 are the eigenvalues of the initial tidal field and determine the nature of the 
collapse [135,136]. In particular, one-dimensional collapse along, say, the third eigen-direction 
is characterised by Ai = = A2 and by A3 < 0. In that case, the pancake singularity is reached 
as a — > — I/A3. Spherically symmetric collapse, on the other hand, has Ai = A2 = A3 = A < 0. 
Then, we have a point-like singularity when a — > — 1/A. 



stages of the nonlinear regime, when the effects of the fluid pressure are negligible. Assuming that 
the contraction is driven by non-baryonic CDM, means that we can (in principle) extend the domain 
of the Zeldovich approximation beyond the above mentioned mildly nonlinear stage. 
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Substituting, the above expressions into the right-hand side of Eqs. (3.4.13) and (3.4.14), we 

obtain the solutions 

^l + a A 2 )(l + a A 3 ) 



B 1 = B° 



Bo = Br, 



and 



B, 



Bl 



(l + aA 2 )(l + aA 3 ) _ 

(1 + aoAi)(l + a A 3 ) 

(l + aAi)(l + aA 3 ) 

1 + o Ai)(l + a A 2 ) 



a 

'Oo 

a 
a 



(3.4.18) 
(3.4.19) 

(3.4.20) 



(1 + oAi)(1 + oA 2 ) 

with the zero suffix corresponding to a given time during the protogalactic collapse. Note that 
the ratios in parentheses reflect the magnetic dilution due to the background expansion, while 
the terms in brackets monitor the increase in the field's strength caused by the collapse of 
the protogalactic cloud. According to the above solution, when dealing with pancake collapse 
along the third eigen-direction, the B 3 -component decays as a -2 , while the other two increase 
arbitrarily. Alternatively, during a spherically symmetric contraction the B-field evolves as 



B = B 



1 + a A \ /oq\ 



1 + aA 



a / 



(3.4.21) 



Here, all the magnetic components diverge as we approach the point singularity (i.e. for 
a — > — 1/A). Comparing the two results, we deduce that the anisotropic (pancake) collapse 
leads to a stronger increase as long as A 3 < A. The latter is always satisfied, provided that 
the initial conditions are the same for both types of collapse, given that A 3 = #/(l — ao^) and 
A 3 = #/(3 — Ooi?) - see expression (3.4.15) above. 

The above given qualitative analysis indicates that a magnetic field trapped in an anisotropi- 
cally contracting protogalactic cloud will increase beyond the limits of the idealized spherically 
symmetric scenario. Note that this type of amplification mechanism appears to be the only 
alternative left if the galactic dynamo (see § 4.2.1 below) fails to operate. Quantitatively, 
the achieved final strength depends on when exactly the backreaction of the B-field becomes 
strong enough to halt the collapse [25]. Thus, the longer the anisotropic collapse lasts, the 
stronger the residual B- field. The analytical study of [134], in particular, showed that (real- 
istically speaking) the anisotropy could add one or two orders of magnitude to the magnetic 
strengths achieved through conventional isotropic compression. These results are in very good 
agreement with numerical studies simulating shear and tidal effects on the magnetic evolution 
in galaxies and galaxy clusters [137,138,139]. 
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4 Magnetogenesis in conventional FRW models 



In order to operate successfully, the galactic dynamo needs magnetic seeds that satisfy two 
specific requirements. The first refers to the (comoving) coherence length of the initial -B-field 
and the second is related to its strength. The scale must not drop below ~ 10 Kpc. The 
strength typically varies between ~ 1CT 12 G and ~ 10~ 22 G, depending on the efficiency of 
the dynamo amplification. At first, these requirements may seem relatively straightforward 
to fulfill. Nevertheless, within classical electromagnetism and conventional FRW cosmology, 
magnetic seeds with the aforementioned desired properties are very difficult to produce. 



4-1 The Friedmann models 



Current observations, primarily the isotropy of the CMB, strongly support a universe that 
is homogeneous and isotropic on cosmological scales. In other words our universe seems to 
be described by the simplest cosmological solution of the Einstein field equations, the FRW 
models. Before proceeding to discuss the magnetic evolution on FRW backgrounds, it helps 
to summarise some basic features of these models. 



4-1.1 The FRW Dynamics 

The high symmetry of the Friedmann models means that all kinematical and dynamical vari- 
ables are functions of time only, while every quantity that represents anisotropy or inhomo- 
geneity vanishes identically. Thus, in covariant terms, an FRW model has = 3H(t) ^ 0, 
Cab — = oo a = A a , E a b = = H a b, where H = a/ a is the familiar Hubble parameter. The 
isotropy of the Friedmann models also constrains their matter content, which can only have 
the perfect-fluid form (with p = p(t) and p = pit)). In addition, due to the spatial homogene- 
ity, all orthogonally projected gradients (e.g. D a p, D a p, etc - see § 3) are by definition zero. 
This means that the only nontrivial equations left, are the FRW version of Raychaudhuri's 
formula, the equation of continuity and the Friedmann equation. These are given by [121] 

H = -H 2 - \n{p + 3p) + \ A , p = -3H(p + p) (4.1.1) 

O 3 

and 

H 2 = -k P -- + -A, 4.1.2 

6 a z 6 

respectively. Note that K = 0, ±1 is the 3-curvature index. The latter is associated to the 
Ricci scalar (1Z) of the spatial sections by means of the relation 1Z = 6K/a 2 [121]. 

In FRW spacetimes with non-Euclidean spatial geometry, the scale factor also defines the 
curvature scale (Xk = a) of the model. This marks the threshold at which the effects of spatial 
curvature start becoming important (e.g. see [140]). Lengths smaller than the curvature scale 
are termed subcurvature, while those exceeding \ K are referred to as supercurvature. The 
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former are essentially immune to the effects of spatial geometry, with the latter dominating 
on supercurvature lengths. The relation between the curvature scale and the Hubble radius is 
determined by Eq. (4.1.2). In the absence of a cosmological constant, the latter reads 

where Xh = H~ x and Q p = np/3H 2 are the Hubble radius and the density parameter re- 
spectively. Hence, hyperbolic 3-geometry (i.e. K = —1) ensures that X K > \ H always, with 
\k — > oo as Q p — y 1 and \ K —¥ \ H for Q p —¥ 0. In practice, this means that supercurvature 



scales in spatially open FRW cosmologies are always outside the Hubble radiuslHI This is not 
the case in closed models, where A^ > Xh when Q p < 2 and Aj^ < Ajy if O p > 2. Finally, we 
note that, since the curvature scale simply redshifts with the expansion, the importance of 
spatial geometry within a comoving region does not change with time. 



4-1-2 Scale-factor evolution in FRW models 

In order to close the system (4.1.1), one needs to introduce an equation of state for the 
matter. Here, we will only consider barotropic perfect fluids, mainly in the form of non- 
relativistic 'dust' or isotropic radiation (with p = and p = p/3 respectively). When w = p/p 
is the (constant) barotropic index of the medium, the continuity equation (see (4.1.1b)) gives 
p oc a~ 3( - 1+w \ Then, setting K = = A and normalising so that a(t = 0) = 0, we obtain 

ft \ 2/3(1+™) 

a = a °{-f) > ( 4 - L4 ) 

when w ^ — 1. For non-relativistic matter with w = (e.g. baryonic dust or non-baryonic cold 
dark matter), we have the Einstein-de Sitter universe with a oc t 2 ^ 3 . Alternatively, a oc t}l 2 in 
the case of relativistic species (e.g. isotropic radiation) and a oc t 1 ^ 3 for a stiff medium with 
w — 1. When w = —1/3, which corresponds to matter with zero gravitational mass, the above 
leads to 'coasting' expansion with a oc t. Solution (4.1.4) does not apply to the w — — 1 case. 
There, both p and H are constant to ensure de Sitter-type inflation with a oc e H °^' ta \ 

When the FRW spacetime has non-Euclidean spatial geometry it helps to use conformal rather 
than proper time. Then, for K = +1, A = and w ^ —1/3 Eqs. (4.1.1), (4.1.2) combine 
to [121] 

where i] is the conformal time (with f] = 1/a). Also, (1 + 3w)r]/2 G (0,7r) and we have 
normalised so that a(r] — > 0) — > 0. For non-relativistic matter w = and the above solution 



11 In Friedmann universes with hyperbolic spatial sections, the particle horizon can exceed the Hubble 
scale (e.g. see [141]). This means that wavelengths larger than the curvature radius of an open FRW 
cosmology can be causally connected, despite the fact that they always lie outside the Hubble length. 
Here, to keep things simple, we will treat the Hubble scale as our causal horizon as well. 
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reduces to a oc sin 2 (77/2) [142], while we obtain a oc sin 77 if radiation dominates. When 
w = —1/3 one can no longer use solution (4.1.5). Instead, Eq. (4.1.1a) leads immediately to 
the familiar coasting-expansion phase with a oc t. Finally, in the w — — 1 case, expression 

(4.1.1b) ensures that p = p = constant and (4.1.2) leads to a(l + J3/p H) oc e ^ {po/3)t . 



Applied to FRW cosmologies with hyperbolic spatial geometry, zero cosmological constant 
and w 7^ —1/3, the same analysis leads to 

f sinh[(l + 3 W )V2n 2/(1+3w) ,, lfi . 

a = ao isinh[(l + 3 W )V2]| ' (4L6) 

where now (1 + 3w)t]/2 > [121]. Assuming pressure-free dust and normalising as before, we 
find a oc sinh (i]/2) [142]. On the other hand, solution (4.1.6) implies a oc sinh?7 for a open 
FRW universe dominated by relativistic species. We finally note that, similarly to the K — +1 



case, the system (4.1.1), (4.1.2) ensures that a oc t when w = —1/3 and o(l + J3/p H) oc 

e y(po/3)t for w = — 1. A special spacetime with open spatial geometry is the vacuum Milne 
universe, where the absence of matter guarantees that a = t. 



4-2 Late vs early-time magnetogenesis 



The various mechanisms of magnetogenesis have been traditionally classified into those op- 
erating at late times, that is after recombination, and the ones that advocate an early (pre- 
recombination) origin for the 5-field. In either case, the aim of the proposed scenarios is to 
produce the initial magnetic fields that will successfully seed the galactic dynamo. 



4-2.1 The galactic dynamo paradigm 

The belief that some kind of nonlinear dynamo action is responsible for amplifying and sus- 
taining the galactic magnetic fields has long roots in the astrophysical community [23,24,25]. 
Dynamos provide the means of converting kinetic energy into magnetic energy and the reader 
is referred to [29] for a recent extended review. Nevertheless, one can get a quick insight of 
how the mechanism in principle works by looking at the magnetic induction equation. In the 
Newtonian limit and assuming resistive MHD, the latter reads 

B a = -- $B a + (<7 Q/3 + uj a0 ) B 13 - rWja , (4.2.1) 

where overdots indicate convective derivatives (see also § 3.4.2) and a, (5 =1,2,3. Contracting 
the above along the magnetic field vector and recalling that curl_B a = J ai leads to 

(B 2 ) = -\$B 2 + 2o- a pB a B? + 2^~ 1 d a F a - 2^ J 2 , (4.2.2) 
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with F a = Safi^B^J'^ representing the magnetic Lorentz force. The latter contributes to the 
kinetic energy of the fluid via the Navier-Stokes equation (see expression (3.1.11) in § 3.1.5). 
Following (4.2.2), the action of the Lorentz force can in principle enhance the magnetic energy 
at the expense of the fluid's kinetic energy. The amplification can happen provided that 
the dissipative effects, carried by the last term on the right-hand side of Eq. (4.2.2), are 
subdominant. The first term, on the other hand gives the magnetic increase caused by the 
adiabatic galactic collapse (typically d < in gravitationally bounded systems), while the 
second conveys the effects of the shearing stresses (see § 3.4.2). 

Dynamos are typically powered by the differential rotation of the galaxy. The latter combines 
with the small-scale turbulent motion of the ionised gas causing the exponential increase of the 
large-scale mean 5-field in the plane of the galactic disc. The growth continues until it reaches 
saturation, which typically occurs when B ~ 10~ 6 G and the backreaction of the magnetic 
stresses suppresses any further increase. The amplification factor, however is quite sensitive 
to the specific parameters of the dynamo model. This sensitivity leads to serious uncertainties 
regarding the total ammount of magnetic amplification and has been the subject of ongoing 
discussions. 

The pattern and the orientation of the galactic magnetic field, especially of those seen in spiral 
galaxies, seem to support the dynamo idea. On the other hand, the detection of strong mag- 
netic fields in high-redshift protogalactic structures has raised a number of question regarding 
the role of dynamos [8,9]. In any case, the galactic dynamo needs the presence of an initial 
magnetic seed in order to operate. These seeds must satisfy certain requirements regarding 
their coherence length and strength. The minimum required scale for the magnetic seed is 
comparable to the size of the smallest turbulent eddy. This lies close to 100 pc at the time 
the host galaxy is formed, which translates to a comoving length of approximately 10 Kpc 
before the collapse of the protogalactic cloud. The strength of the seed-field, on the other 
hand, varies, depending on the efficiency of the dynamo amplification and on the cosmological 
model it operates in [26,27,28,29]. Typically, the required values range between 10~ 12 G and 
10~ 22 G. It is conceivable, however, that the lower limit could be brought down to 10~ 30 G 
in spatially open or dark-energy dominated FRW universes [30]. Note that, in the absence of 
the dynamo, protogalactic collapse (spherically symmetric or anisotropic - see § 3.4.2) seems 
the only alternative means of magnetic amplification. Then, 5-seeds as strong as 10~ 9 G may 
be needed in order to meet the observations. So, provided that galactic dynamos work, the 
question is where do the initial magnetic seeds come from? 



4-2.2 Late-time magnetogenesis 

Post recombination mechanisms of magnetic generation appeal to astrophysical processes and 
battery-type effects. It has been proposed, in particular, that the Biermann-battery mechanism 
can produce seed S-fields, which the dynamo could subsequently amplify on galactic scales 
and to the observed strengths. The Biermann effect [143], which was originally discussed in 
the stellar context, exploits differences between the electron and the ion acceleration that are 
triggered by pressure gradients. These will first give rise to electric currents and subsequently 
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lead to magnetic fields by induction. 

The literature contains several alternative scenarios using battery-type mechanisms to generate 
magnetic seed-fields in the post recombination era. Supernovae explosions of the first stars, 
for example, could eject into the interstellar medium 5-fields that could seed the galactic 
dynamo [144,145]. Active galaxies and Active Galactic Nuclei (AGN) can also channel away 
jets of magnetised plasma [146,147]. Thermal-battery processes operating in (re)ionisation 
fronts can also lead to magnetic seeds that can sustain the dynamo [148,149]. Analogous results 
could be achieved through turbulent motions or shocks developed in collapsing protogalactic 
clouds [150,151]. 

Nevertheless, while Biermann battery effects can produce the seed fields that the dynamo 
will subsequently amplify to the observed strengths, the whole process operates on galactic 
scales. For this reason, it is less straightforward to invoke the Biermann mechanism when 
trying to explain the magnetic fields found in galaxy clusters. To a certain extent, this also 
weakens the overall position of the Biermann battery as a likely candidate for generating the 
galactic magnetic fields. Indeed, the possibility that the galactic and the cluster 5-fields have 
a different origin seems rather unlikely, in view of their similarities. 



4-2.3 Early-time magnetogenesis 

The idea that cosmic magnetism might have pre-recombination origin is attractive because 
it makes the widespread presence of magnetic fields in the universe easier to explain. Espe- 
cially the origin of the fields observed in high-redshift proto-galactic condensations. However, 
generating cosmological 5-fields that will also successfully seed the galactic dynamo is not a 
problem-free exercise. 

In the early 1970s, Harrison proposed that battery- type effects, operating during the radiation 
era, could generate B- fields with strengths capable of sustaining the galactic dynamo [152]. 
The mechanism is based on conventional physics and does not need any new postulates. The 
disadvantage in Harrison's idea is that it requires significant amounts of primordial vorticity. 
The latter is essentially absent from the standard cosmological model. Note that the possibility 
of simultaneously generating both vorticity and magnetic fields in the late radiation era and 
around recombination (when the tight-coupling between photons and baryons is relaxed) was 
recently investigated in [153,154]. 

An alternative approach is to generate the magnetic seeds during phase transitions early in the 
radiation era (see § 5 below). There are problems, however, primarily related to the coherence 
length of the initial 5-field. The difficulties arise because the size of the post-inflationary mag- 
netic seeds, namely those created between inflation and (roughly) recombination is typically 
too small and will destabilise the dynamo. The reason is causality, which confines the scale 
of the field within that of the horizon at the time of magnetogenesis. For example, -B-fields 
produced during the electroweak phase transition have coherence lengths of the order of the 
astronomical unit. The size of the magnetic field can increase if the host plasma has some 
degree of MHD turbulence. In such environments "cascade" processes are known to occur, 
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whereby certain ideally conserved quantities flow from larger towards smaller scales (direct 
cascade) or the other way around (inverse cascade). In three dimensional MHD turbulence, the 
total (kinetic plus magnetic) energy cascades toward smaller scales, where it is dissipated by 
viscosity and resistivity. However, the other important ideal invariant, the magnetic helicity, 
inverse-cascades towards larger scales. The magnetic helicity is defined by the integral [155] 

H M = ^JA a B a dV, (4.2.3) 

where A a is the electromagnetic vector potential (recall that B a = cm\A a ) and is equivalent 
to the Chern-Simon number of particle physics [32]. Besides being an ideal invariant, the 
magnetic helicity is also asymptotically conserved within the resistive MHD approximation. 
Physically, %m describes the topology of the field lines, that is their degree of withering 
and twisting [156]. As mentioned before, magnetic helicity inverse-cascades and evolves from 
smaller to larger scales [155], while its conservation has profound effects in the operation 
of MHD dynamos [157]. The aforementioned inverse-cascade effect makes primordial helicity 
very important, because it allows the magnetic energy to shift from smaller to larger scales, 
as the system tries to minimize its energy while conserving magnetic helicity. For example, 
Pouquet et al carried out a study in which nonhelical kinetic energy and maximally helical 
magnetic energy were injected into the plasma at a constant rate [158]. The outcome was a 
well defined wave of magnetic energy and helicity, propagating from smaller to larger scales. 
Similar results were also obtained in the case of steady turbulence [31,32,33,34,159,160,161] 
and for freely decaying MHD turbulence [162]. 

Although helical magnetic fields can enhance their original length, inflation seems to be the 
only effective solution to the scale problem faced by fields generated during the early radiation 
era. The reason is that inflation can naturally generate correlations on superhorizopn lengths. 
There are still problems, however, this time with the magnetic strength. In particular, 5-fields 
that has survived a period of standard de Sitter inflation are typically to weak to sustain the 
galactic dynamo. 



4-3 Typical inflation produced magnetic fields 



Inflation is known to produce long wavelength effects from microphysical processes that oper- 
ate well inside the Hubble radius. For this reason, inflation has long been seen as the best can- 
didate for producing large-scale, cosmological magnetic fields. Here, we will look at scenarios 
operating within standard electromagnetic theory and conventional FRW models. Alternative 
approaches are given in § 6. 



4-3.1 Quantum-mechanically produced magnetic seeds 

The inflationary paradigm provides the dynamical means of producing long-wavelength elec- 
tromagnetic fluctuations, by stretching subhorizon-sized quantum mechanical fluctuations to 
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superhorizon scales. Roughly speaking, quantum fluctuations in the Maxwell field are excited 
inside the horizon and cross the Hubble horizon approximately 

N = N{\) = 45 + ln(A) + \ In (^) + \ (g|) , (4.3.1) 

e-folds before the end of the de Sitter phase [163]. In the above A is the comoving scale of 
the mode (measured in Mpc and normalised to coincide with the mode's current physical 
length), M is the scale of inflation and Trh is the reheat temperature (both measured in 
GeV). Assuming that p is the energy density of the electromagnetic mode, then 

%-H*. (4.3.2) 

at the first horizon crossing. Once outside the Hubble radius, the aforementioned quantum- 
mechanically excited modes are expected to freeze-out as classical electromagnetic waves. The 
latter, which initially appear like static electric and magnetic fields, can subsequently lead to 
current-supported magnetic fields. This happens after the modes have re-entered the horizon 
in the radiation era, or later during the dust-dominated epoch. Note that, after the second 
horizon crossing, the currents of the highly conductive plasma will also eliminate the electric 
component of the Maxwell field, leaving the universe permeated by a large-scale 5-field of 
primordial origin. 

The fast expansion of the de Sitter phase means that, by the end of inflation, the initial 
electromagnetic quantum fluctuations have achieved correlation lengths much larger than the 
current size of the observable universe. Thus, inflation produced 5-fields have no scale problem 
whatsoever. Nevertheless, magnetic seeds that have survived a period of de Sitter expansion 
are generally too weak to sustain the dynamo. In particular, the typical strength of the residual 
B- field (in today's values) is less than 10~ 50 G [35]. To understand why and how this happens, 
we first need to consider the linear magnetic evolution on FRW backgrounds. 



4-3.2 The adiabatic magnetic decay 

The evolution of large-scale electromagnetic fields on FRW backgrounds depends on the elec- 
tric properties of the medium that fills the universe. Here, we will consider the two limiting 
cases of poorly and highly conductive matter. For any intermediate case, one needs a model 
for the electrical conductivity of the cosmic medium. 

In poorly conductive environments, q — > and the electric currents vanish despite the presence 
of nonzero electric fields (see Ohm's law (3.3.1) in § 3.3). Then, the wave equation (3.2.7) 
linearises to 

B a - D 2 B a = -5HB a + i k P {1 + 3w)B a - AH 2 B a - TZ ab B b , (4.3.3) 

where H = a/a is the Hubble parameter of the unperturbed model. To simplify the above we 
introduce the rescaled magnetic field B a = a 2 B a and employ conformal, rather than proper, 
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time [35]. Then, on introducing the harmonic splitting B a = Y,nB(n)QaP ~ with D a B( n ) = 0, 
expression (4.3.3) takes the compact form 12 1 

B'{ n) + n 2 B H = -2KB {n) , (4.3.4) 

with the primes denoting conformal-time derivatives and K = 0, ±1 [118]. The above describes 
the linear evolution of the rescaled magnetic field on a Friedmannian background with any 
type of spatial curvature. Note the magneto-curvature term on the right-hand side of (4.3.4), 
which results form the direct coupling between the 5-field and the geometry of the 3-space. 
The interaction is monitored by the Ricci identities and reflects the fact that we are dealing 
with an energy source of vector nature within a geometrical theory of gravity. We will discuss 
the implications of this interaction, which is largely bypassed in the literature, for the evolution 
of large-scale magnetic fields in § 4.4. 

When the FRW host has Euclidean spatial hypersurfaces, the 3-curvature index is zero 
(i.e. K — 0) and expression (4.3.4) assumes the Minkowski-like form 

B'( n) +n 2 B (n) =0. (4.3.5) 

The latter accepts the oscillatory solution 

B(n) = Ci sin(m/) + C 2 sin(n?7) , (4.3.6) 

with the integration constants depending on the initial conditions. Then, recalling that B( n ) = 
a 2 S( n ), the above given solution translates into 



5 (n) = (^) 2 |< ) sin(nr /0 ) + i 



B' (n) + 2 



4 (^) 2 {4 n) cos(m/o)-- 
a J \ n 




b'^ + 2 r-\ Bt ] 



cos(n?7o) \ sin(nry) 

sin(ri77o) > cos(n?7) . (4.3.7) 



This guarantees an adiabatic (-B( n ) oc a~ 2 ) depletion for the magnetic field, irrespective of the 
equation of state of the matter, as long as the background spacetime is a spatially flat FRW 
model and the electrical conductivity remains very poor. 

The adiabatic magnetic decay is also guaranteed in highly conductive environments, namely 
at the ideal-MHD limit. There, <^ — > oo and, according to Ohm's law - see Eq. (3.3.1) in § 3.3, 
the electric field vanishes in the frame of the fluid. As a result, when linearised around a FRW 



12 We use pure- vector harmonics that satisfy the constrains Qa = = D a Qa and the associated 
Laplace-Beltrami equation, namely D 2 Qa ' = —(n/a) 2 Qa- Following [164], the (comoving) eigen- 
values depend on the background spatial curvature according to re 2 = v 2 + 2K, where v represents 
the associated wavenumber. Also, n has a continuous spectrum, with n 2 > 0, when K = 0, —1 and 
a discrete one, with n 2 > 2, if K = +1. 
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background, Faraday's law (see Eq. (3.2.3) in § 3.1) reduces to 

B a = -2HB a . (4.3.8) 

The latter guarantees that B a oc a~ 2 on all scales, regardless of the equation of state of the 
matter and of the background 3-curvature. 

The universe is believed to have been a very good electrical conductor throughout its classical 
Big-Bang evolution, at least on subhorizon scales. During inflation, on the other hand, the 
conductivity is expected to be very low. However, as the universe leaves the inflationary phase 
and starts reheating, its conductivity grows. So, by the time we have entered the radiation era, 
the currents have eliminated the electric fields and frozen their magnetic counterparts in with 



the matter J 13 1 These arguments essentially guarantee that the set (4.3.5) and (4.3.8) monitors 
the evolution of cosmological magnetic fields throughout the universe's lifetime. This in turn 
has led to the widespread belief that the adiabatic decay-rate of cosmological magnetic fields 
is ensured at all times, unless classical electromagnetism is modified or the FRW symmetries 
are abandoned. As we will see in § 4.4, however, this is not necessarily the case. 



4-3.3 The residual magnetic field 

The immediate implication of (4.3.7) is that magnetic fields that survived a period of typical 
de Sitter-type inflation have been drastically diluted by the accelerated expansion of the 
universe. Together with (4.3.8), this means that 5-fields of primordial origin are too weak to 
be of astrophysical relevance today. To demonstrate the dramatic magnetic depletion during 
the de Sitter phase we follow [35]. As a first step, recall that the relative energy density stored 
in the n-th magnetic mode at the (first) horizon crossing is {pb/p)hc — {H/Mp{) 2 . Here, 
Pb = -B( 2 „) 5 P is the energy density of the background universe and Mpi is the Planck mass. 
During inflation the total energy density is dominated by that of the vacuum (i.e. p ~ M 4 , 
with M representing the energy scale of the adopted inflationary scenario). Consequently, the 
relative strength of the n-th magnetic mode at horizon crossing is (pb/p)hc — (M/M P i) 4 . 
Also, throughout inflation the universe is believed to be a very poor electrical conductor. 
This means that any magnetic field that may be present at the time decays adiabatically 
(see solution (4.3.7) in § 4.3.2). As a result, B 2 , = (B? n ,) HC e~ 4N by the end of inflation, 
with N = ln(ajN/aHc) representing the number of e-folds between horizon crossing and 
the end of the de Sitter era. This number depends on the scale of the mode in question 
and, in typical inflationary scenarios, is given by Eq. (4.3.1). Using the latter and recalling 
that (pb/p)rh = (pb / p) in(Thh / M) 4 ' 3 is the relative change of the magnetic energy density 
between the end of inflation proper and that of reheating, we find that [35] 

-104 \ -4 




l(T iU4 A- 4 , (4.3.9) 



13 Causality ensures that the post-inflationary currents are confined within the horizon. This means 
that outside the Hubble radius the low conductivity assumption still holds. 
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at the onset of the radiation era. Note that p 7 ~ prh — Trh represents the energy density of 
the relativistic species and A is the current (comoving) scale of the B- field. Also, the above ratio 
is independent of the energy scale of the adopted inflationary scenario and of the associated 
reheat temperature. Moreover, given that p#, p 7 oc a~ 4 throughout the subsequent evolution 
of the universe, the same ratio remains unchanged until the time of galaxy formation. 

Once the scale of the magnetic mode is given, we can use (4.3.9) to evaluate the residual 
strength of any primordial 5-field that underwent an era of (typical) de Sitter inflation. For 
example, in order to operate successfully, the galactic dynamo requires magnetic seeds with 
a minimum coherence scale of approximately 10 Kpc. Substituting this scale into Eq. (4.3.9) 
and recalling that p 7 ~ 10~ 51 GeV today, we find that the corresponding magnetic field 
has strength of ~ 10~ 53 G [35]. This value is well below the galactic dynamo requirements, 
which leads to the conclusion that magnetic fields that have survived a period of standard, de 
Sitter-type inflationary expansion are (for all practical purposes) astrophysically irrelevant. 



4-4 Magnetic amplification in conventional FRW models 



The "negative" results of the previous section have been been widely attributed to the con- 
formal invariance of Maxwell's equations and to the conformal flatness of the Friedmannian 
spacetimes. The two have been thought to guarantee an adiabatic decay-rate for all large-scale 
magnetic fields at all times. This, in turn, has led to the widespread perception that infla- 
tion produced 5-fields are astrophysically unimportant, unless standard electromagnetism is 
modified or the FRW symmetries are abandoned. 



4-4-1 Superadiabatic amplification 

Strictly speaking, the adiabatic magnetic depletion seen in solution (4.3.7) of § 4.3 has only 
been proved on Friedmannian backgrounds with Euclidean spatial sections. Although it is true 
that all three FRW universes are conformally flat, they are not the same. There are differences 
in their 3-curvature, which mean that only the spatially flat model is globally conformal to 
Minkowski space. For the rest, the conformal mappings are local [165,166]. Another way of 
putting it is that, when dealing with spatially curved Friedmann universes, the conformal 
factor is no longer the cosmological scale factor but has an additional spatial dependence. All 
these imply that the wave equation of the rescaled magnetic field (B a = o 2 B a ) takes the simple 
Minkowski-like form (4.3.5), which guarantees an adiabatic decay for the actual 5-field, only 
on FRW backgrounds with zero 3-curvature. In any other case, there is an additional curvature- 
related term (see expressions (4.3.3) and (4.3.4) in § 4.3), reflecting the non-Euclidean spatial 
geometry of the host spacetime. As a result, when linearised around an FRW background with 
nonzero spatial curvature, the magnetic wave equation reads [118] 

&( a) + (n 2 ±2)B {n) = 0, (4.4.1) 
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with the plus and the minus signs indicating the spatially closed and the spatially open model 
respectively. Recall that in the former case the eigenvalue is discrete (with n 2 > 2), while in 
the latter it is continuous (with n 2 > 0). In either case, the curvature- related effects fade away 
as we move down to successively smaller scales (i.e. for n 2 ^> 2). 

According to (4.4.1), on FRW backgrounds with spherical spatial hypersurfaces, the S-field 
still decays adiabatically. The curvature term only modifies the frequency of the magnetic 
oscillation in accord with the solution [167] 



B 



(n) 



C\ sin (yn + 2r/j + C 2 cos f \Jn + 2 



aoV 
a J 



(4.4.2) 



Overall, the adiabatic decay-rate of the -B-field remains. Also, as expected, the smaller the 
scale the less important the role of the background 3-geometry. 

The standard picture, and the adiabatic-decay law, change when the background FRW model 
has open spatial sections. In particular, the hyperbolic geometry of the 3-D hypersurfaces 
alters the nature of the magnetic wave equation on large enough scales (i.e. when < n 2 < 2). 
These wavelengths include what we may regard as the largest subcurvature modes (i.e. those 
with 1 < n 2 < 2) and the supercurvature lengths (having < n 2 < 1). Recall that eigenvalues 
with n 2 = 1 correspond to the curvature scale with physical wavelength X = \k = a (see 
§ 4.1.1). 

Following [167,168,169], we introduce the scale-parameter k 2 = 2 — n 2 , with < k 2 < 2. 
Then, k 2 = 1 indicates the curvature scale, the range < k 2 < 1 corresponds to the largest 
subcurvature modes and their supercurvature counterparts are contained within the 1 < k 2 < 
2 interval. In the new notation and with K = —1, Eq. (4.4.1) reads 



while its solution leads to 



B'L - k 2 B {n) = 



B(k) = C\ smh(\k\r]) + C 2 cosh(\k\r]) . 



(4.4.3) 



(4.4.4) 



Written with respect to the actual magnetic field, the above takes the form 



B 



(k) 




B'M + 2 ( a - 



B'rS k) + 2 




\k\(v-Vo) 



(4.4.5) 



Magnetic fields that obey the above evolution laws can experience superadiabatic amplification 
without modifying conventional electromagnetism and despite the conformal flatness of the 
FRW hostHU For instance, throughout the radiation and the dust eras (as well as during 



14 Superadiabatic amplification does not imply amplification per se, but decay at a slower than the 
adiabatic pace. The concept was originally introduced in gravitational-wave studies [170,171]. 
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reheating), the scale factor of a FRW universe with hyperbolic spatial geometry evolves as 
a oc sinh(r]) and a oc smh 2 (r]/2) respectively (see solution (4.1.6) in § 4.1.2). Focusing on the 
curvature scale, for simplicity, we may set \k\ = 1 in (4.4.4). It is then clear that, on that scale, 
the magnetic mode never decays faster than Bm oc a -1 [168]. In other words, large-scale in- 
fields are superadiabatically amplified throughout the post-inflationary evolution of an open 
Friedmann universe. 

Although in the above examples we only considered the cases of radiation and dust, the 
amplification effect is essentially independent of the type of matter that fills the universe. In 
particular, i?-fields in open FRW models containing a barotropic medium with p/p ^ —1/3 



are superadiabatically amplified on large scales [169] J 15 1 This means that the mechanism also 
operates during reheating (when p = 0) and also throughout a phase of slow-roll inflation, 
namely in spatially open FRW models with a false-vacuum equation of state (i.e. when p = 
—p). In the latter case, the background scale factor evolves as [121] 

/l- e 2 w\ 
a = a (j— - j e™ , (4.4.7) 

where ^,^ < 0. Substituting the above into Eq. (4.4.5), we find that near the curvature scale 
(i.e. for \k\ — >■ 1) the magnetic evolution is given by 



Jv\ I a \ ,r.*-v l a ° 



S(1)=C3(1 _ e ^j +C4e -,^j\ (4.4.8) 

with C 3 , C4 depending on the initial conditions [167,168]. This result also implies a superadi- 
abatic type of amplification for the 5-field, since the dominant magnetic mode never decays 
faster than B(i) oc a -1 . The adiabatic decay rate is only recovered at the end of inflation, 
as rj — > 0~ . It should be noted that the magnetogeometrical interaction triggering the above 
described effect is possible because, when applied to spatially curved FRW models, inflation 
does not lead to a globally flat de Sitter space. Although the inflationary expansion dramat- 
ically increases the curvature radius of the universe, it does not change its spatial geometry. 
Unless the universe was perfectly flat from the beginning, there is always a scale where the 
3-curvature effects are important. It is near and beyond these lengths that primordial B- fields 
are superadiabatically amplified. 



15 Mathematically, the easiest way of demonstrating the amplification effect is by adopting the Milne 
universe as our background spacetime. The latter corresponds to an empty spacetime with hyperbolic 
spatial geometry (see § 4.1.2) and can be used to describe a low density open universe. There, the 
scale factor evolves as a oc e v , which substituted into solution (4.4.5) leads to [168] 

Consequently, all magnetic modes spanning scales with < k 2 < 2 are superadiabatically amplified. 
Close to the curvature scale, that is for k 2 — > 1, the dominant magnetic mode is Bn\ oc a -1 ; a rate 
considerably slower than the adiabatic a _2 -law. The latter is only restored at the k 2 = limit, namely 
on small enough scales where the curvature effects are no longer important. Stronger amplification 
is achieved on supercurvature lengths, with B/ k \ oc a^ 2 ~ 2 at the homogeneous limit (i.e. as k 2 — > 2). 
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The magnitude of the residual magnetic field is calculated in a way analogous to that given in 
§ 4.3.3. Now, however, there is an additional parameter due to the non-Euclidean geometry of 
the FRW background. In particular, near the curvature scale - where B oc a -1 , we find that 

at the beginning of the radiation era [169]. Note that, in deriving the above, we also used the 
auxiliary relation (pb/p)rh = {pb/p)in{M/Trh)^ 3 - The latter provides the relative change 
in the energy density of the (superadiabatically amplified) magnetic field between the end 
of inflation and that of reheating. Comparing (4.4.9) to expression (4.3.9), one can see that 
the (superadiabatic) magnetic amplification is already substantial by the end of reheating, 
despite its dependence on the energy scale of the inflationary model and of the corresponding 
reheat temperature. Moreover, large-scale 5-fields are superadiabatically amplified during 
the subsequent evolution of the universe. This means that on scales close to the curvature 
radius of our background model, the ratio r = p B j p 7 is no longer constant but increases as 
r oc a 2 oc T~ 2 . Consequently, recalling that X K = \ H /\Jl — Vt is the curvature scale of a 
spatially open FRW cosmology, we obtain 

/ A/f \ 2 

B ° ~ 1(r13 (io*) (1 " fio)1/2 G - (4A10) 

for the present strength of the residual B- field [169]. Therefore, the higher the energy scale 
of inflation the stronger the superadiabatic amplification. On the other hand, the closer the 
density parameter to unity, the weaker the final field. 

Currently, the WMAP reports indicate that 1 - Q < 1(T 2 [172,173,174]. On these grounds, 
and provided that the universe is spatially open, expression (4.4.10) gives 

£ ~ 1(T 14 G, (4.4.11) 

when M ~ 10 14 GeV and 1 — Qq ~ 10~ 2 [169]. The last parameter choice implies a curvature 
radius of the order of 10 4 Mpc at present. These lengths are far larger than 10 Kpc, which is 
the minimum magnetic size required for the dynamo to work. Nevertheless, once the galaxy 
formation starts, the field lines should break up and reconnect on scales similar to that of 
the collapsing protogalaxy. According to (4.4.10), the above quoted magnetic strength will 
increase if the energy scale of inflation is greater than 10 14 GeV. On the other hand, the 
magnitude of the residual 5-field will drop if the current curvature scale of the universe is 
much larger than the Hubble horizon (i.e. for 1 — Q <C 1). Nevertheless, the ^-dependence 
in Eq. (4.4.10) is relatively weak, which means that B- fields capable of seeding the galactic 
dynamo (i.e. with Bq > 10~ 22 G) are possible even when 1 — Qq ~ 10~ 18 (or lower - when the 
scale of inflation is higher than 10 14 GeV). Overall, FRW universes with hyperbolic spatial 
geometry seem capable of sustaining astrophysically important magnetic fields under a fairly 
broad range of initial conditions! 16 1 



16 If the results of [30] are taken at face value, the 1 — Qq difference can drop down to ~ 10 (or even 
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Large-scale, primordial magnetic fields with residual magnitudes like those quoted above are 
far stronger than any other of their conventionally produced counterparts. Such strengths 
are usually achieved outside classical electromagnetism or beyond the standard model (see 
§ 5 and § 6 below) J 17 1 Moreover, although magnetic seeds with strengths < 10~ 10 G cannot 
affect primordial nucleosynthesis or the CMB spectrum, their strength lies within the galactic 
dynamo requirements (see § 4.2.1) and are therefore of astrophysical interest. Finally, we 
should also note very recent reports indicating the presence of coherent magnetic fields in 
empty intergalactic space with strengths intriguingly close to those quoted here [12,13,14,15]. 



lower) and still produce magnetic fields able to sustain the galactic dynamo (i.e. with Bq > 10~ 30 G). 
We also note that to these magnitudes one should add the magnetic amplification during the linear 
and the nonlinear regime of structure formation - see § 3.4.1 and § 3.4.2 respectively. 
17 p r i m ordial magnetic fields can be superadiabatically, or even resonantly, amplified through their 
interaction with cosmological gravitational waves (see [175] and references therein). The former type 
of amplification is typically associated with highly conductive environments, but requires rather large 
amounts of shear anisotropy to operate efficiently. Resonant amplification, on the other hand, occurs 
in media of poor electrical conductivity and can lead to substantially strong I?-fields with relatively 
minimal shear anisotropy. Both mechanisms are essentially nonlinear in nature and their detailed 
discussion lies beyond the limits of this review. 
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5 Magnetogenesis in the standard model 



In this section we review mechanisms of primordial magnetic field generation in the framework 
of the standard model of particle physics and of non-linear and out of equilibrium processes 
that may have happened in the very early universe. In the first subsection we address mag- 
netogenesis during reheating and in the second subsection we review magnetogenesis due to 
the EW (electroweak) and QCD (quantum chromodynamical) phase transitions. Although the 
EW phase transition is likely to have taken place during reheating, due to its importance and 
to the fact that it is framed in particle physics field as well as the QCD transition, we treat it 
together with the latter, in a specific subsection. 



5.1 Magnetogenesis during reheating 

Accepting the inflationary paradigm, the reheating phase of the Universe was one of the 



richest epochs in its evolutionj 18 1 It is usually treated as the intermediate phase between the 
exponential expansion and the radiation dominated expansion, during which almost all the 
matter that constitutes the Universe was produced. This period can be divided roughly into 
two or three stages: preheating, heating and thermalisation. Of these the most interesting ones 
are the first and the third. During the first stage, the dominant effect is parametric particle 
creation. The importance of this process for reheating was first realized in 1990 by Traschen 
and Brandenberger [176] and also by Dolgov and Kirilova [177] and later developed in Refs. 
[178,179,180,181]. The thermalization process is a difficult and complex one. The interested 
reader is referred to specific works, such as, e.g., [182], with references therein. 

In few words, the process of reheating the Universe is due to the profuse creation of par- 
ticles, caused by inflaton oscillations around the minimum of the effective potential. Those 
particles self-interact and ultimately reach a state of thermal equilibrium, when all (or almost 
all) the inflaton energy has been transformed into thermal energy of the created elementary 
particles, with temperature Trh, the so-called reheating temperature. Being a strong out 
of equilibrium process (and also turbulent, according to theoretical and numerical studies 
[183,184,185,186,187,188]), the reheating period is therefore a suitable scenario for primordial 
magnetogenesis. 

It is important to observe that, irrespective of the mechanism that may have generated the 
magnetic field during inflation, the quantity r = Pb/p-y can be diluted during reheating be- 
cause, during that phase, the radiation density increases by a factor of at least e AN , with N 
being the number of e-foldings. So, unless the gauge field is also amplified by the same amount, 
r is likely to decrease during reheating. 



18 The term "reheating" was coined in the first inflationary models, in which the Universe was hot 
before the onset of inflation, and was "reheated" again, after super-cooling during inflation. 
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5.1.1 Parametric resonance 

Although the conformal invariance of the gauge fields is the main drawback for their ampli- 
fication by the expansion of the universe, it also opens up the possibility of amplification by 
parametric resonance, if the conditions in the early universe are favourable. In this sense, the 
preheating stage of reheating offers a suitable scenario of magnetic amplification, through its 
parametric resonance with a scalar field [189] . From the Lagrangian density of scalar electro- 
dynamics conformally coupled to gravity, 

C = -^-F ab F ab - (d a + ieA a ) $t (d a - ieA a ) $ - V (V$) , (5.1.1) 

l07T v j 

where A a is the gauge potential, F a b = d a Ab — d\,A a and $ a complex (i.e., charged) scalar 
field, we obtain the evolution equation for the gauge field. Working in the Coulomb gauge 
(V • A = 0) and writing the homogenous part of the scalar field as $ = exp [iQ (£)] p (£) /y/2, 
the evolution equation of the fluctuations in the spatial component of the gauge potential, 
SA, reads 

V 2 

5A+H5A -6 A + Aire 2 p 2 5 A = 4n5j , (5.1.2) 

a 2 

where <5j is a source term that can be non null when statistical correlations of the electric 
currents are considered [190,191]. In Fourier space and in terms od conformal time, 77 with 
di] = alt J a (t), the homogeneous part of eq. (5.1.2) is 

SAl + u 2 (ri)SAl = 0, (5.1.3) 

with uj\ (f]) = k 2 + 4ire 2 a 2 (77) p 2 (77), and where the primes denote derivatives with respect to 
77. Expression (5.1.3) describes a harmonic oscillator with time dependent frequency, which 
during the oscillation of the complex scalar field can be rewritten as a Mathieu-like equation. 
The solutions of that equation exhibit exponential increase, i.e., they are proportional to 
e Mfc?? (/ifc is the Floquet exponent) for some frequency intervals known as resonant bands. The 
features of the parametric resonance depend strongly on the time evolution of the homogeneous 
scalar field, which in turn depends on the form of V f$^<I>). For potentials of the form V = 

. 2n 

Finelli and Gruppuso found that, for a quadratic potential, parametric resonance 
is efficient when 4-7re 2 p 2 ^> Ai and it is stochastic and broad, with the largest pk occurring 
for small k [178]. For a quartic potential, eq. (5.1.3) becomes a Lame equation [192]. The 
resonance features and pt depend on e 2 /A2 [192]. In particular modes with k 2 <^i A2P0 are 
resonant for 1/2-7T < e 2 /\2 < 3/2tt (po being the initial value of p). For a symmetry-breaking 

potential, i.e., V = m 2 && + A (^(I?) and in the case of m 2 < 0, the gauge field acquires 



A n ($t$y 



an effective mass propotional to Jm 2 /\, a fact that can completely inhibit the resonance in 
an expanding universe. For a general value of m 2 , the gauge coupling affects the resonance 
structure of the scalar field and it is not possible to determine the resonant bands for the 
imaginary part of the scalar field, as it would be the case if the charged scalar fields were not 
coupled to the gauge field. 
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Another possible coupling is the one described by the Lagrangian density 

C = --^-F ab F ab - d a cf>d a ct> -V(<f>)- 9 -6F ab F ab , (5.1.4) 

107T 4 

where may represent an axion or a general pseudo-Goldstone boson. If the scalar field 
performs coherent oscillations, the evolution equation of the transverse circularly polarized 
photons, A±k, is again given by a Mathieu-like equation 

5A^ k + (k 2 ± Airg(f/k) 5A± k = . (5.1.5) 



In this case the resonance occurs when k ~ 4irg<f)uj, with u being the oscillation frequency 
of (which is very small). For V = A0 4 /4 and 47n?/ = 1 (with / being the Peccei-Quinn 
symmetry scale), A±k grows linearly for k/ui <C kitgf [189]. 

All the previous description did not take into account dissipation due to the presence of 
other charged fields, i.e., plasma effects. When they are taken into account (basically in the 
form of electric conductivity), their effect on the magnetic field depends on the wavelengths 
considered. For wavelengths larger than the plasma collision length, the equations acquire 
a damping term proportional to the conductivity, 47iaq5A' k , with a = 0(77) being the scale 
factor of the universe and q the conductivity. If aq is constant (as is the case considered in the 
literature, where <; oc T oc 1/a [191]) and larger than the Floquet exponent, the parametric 
resonance could be completely suppressed. For wavelengths shorter than the plasma collision 
length, the plasma frequency changes as uj\ (77) — > io\ (77) + 4ne 2 n (77) /m, with n (77) being the 
number density and m the mass of the plasma particles. This term plays the role of an effective 
mass that decays as a -3 , thus allowing for resonance, especially for large coupling constants. 

Whether or not primordial magnetic fields are amplified by parametric resonance during pre- 
heating, depends on the existence of an oscillating scalar or axion field, which the e.m. field is 
coupled to. If, for example, that charged scalar field is the inflaton, the maximum amplifying 
factor obtained for the gauge field is ~ 10 12 , which is not enough to give the minimum seed 
fields that can sustain a dynamo action [189]. 

Exponential growth of large scale magnetic fields could also be achieved by considering the 
Lagrangian [35] 

pi a 

£ = -1 „ -p, — -.F a bF a — — RA a A a + Cm nation , (5.1.6) 

107rG 4 Z 

with R the scalar curvature and j3 a real constant. From the action S = J d A xC, it is obtained 
that the evolution equation for the Fourier components of the magnetic field is given by [193], 



SB' k ' 



k 2 + e(ri) 5B k = 0, (5.1.7) 



with 



V 



6(77) = 6/3—. (5.1.S 

a 
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When the inflaton enters the oscillatory regime, the scalar curvature, given by 



R ** 



M 2 Pl 



AV{4>) 



a 2 



(5.1.9) 



also oscillates and this fact can lead to efficient amplification of the magnetic fluctuations. Con- 
sidering V ((f)) = (1/2) m 2 <j) 2 , we have that during reheating, the inflaton condensate evolves 

as 

Mpi , x , 

sin (mi) , (5.1.10) 



_'?>nmt 
and therefore the scalar curvature as 



i?~ — -[l-3cos(2mt)l ■ (5.1.11) 

3t 2 



Defining B k = a l l 2 B k , eq. (5.1.7) recasts into [193] 

d 2 5B, 



, _, + [C k -2qcos(2mt)]5B k = 0, (5.1.12) 



with 



when ft > 0, and 



3 ± m 2 a 2 ' ' m 2 t' 2 



C k = ~q + ^, g = 12^-, (5.1.13) 



2 k 2 \B\ 
C* = -t?+-Vj, 9 = 12^, (5.1.14) 

3 m z a z m z t z 

for (3 < 0. Parametric resonance occurs when q < 1 initially [193], but to get relevant growth 
of £>*; one needs \f3\ > 1. However, for \f3\ ^> 1 the growth of the magnetic fluctuations is sup- 
pressed [194]. In the former case, as the super-horizon 5B k modes are exponentially suppressed 
during inflation, we do not expected to obtain high intensities from parametric amplification. 
On the other hand, for sub-Hubble scales, the suppression is weaker and consequently mag- 
netic fields can be amplified during preheating. For the latter case, namely when ft < 0, the 
amplification is mainly due to inflation rather than parametric resonance. When finite elec- 
tric conductivity is taken into account, a term of the form — q e a5B' k is added to the r.h.s of 
eq. (5.1.7), which counteracts the parametric amplification of the fields. In summary, despite 
the exponential growth of the magnetic fluctuations due to parametric resonance, the main 
amplification occurs during inflation [193]. 

The possibility of further amplification by parametric resonance during reheating of a seed 
hypermagnetic field generated during inflation, was investigated by Dimopoulos et al [195]. 
However the authors also concluded that the fields do not grow substantially during preheating. 



5.1.2 Magnetogenesis by stochastic electric currents 

Another possibility to induce magnetic fields during reheating is tied to the fact that abrupt 
changes in the metric at that stage may result in the abundant creation of charged particles. 
This could generate stochastic currents, which would eventually decay into the Maxwell field 
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[190,191]. As the inflation is a gauge singlet, it will not decay directly into charged species, so 
this mechanism assumes the existence of another field, a charged one, that is in its vacuum 
state during inflation. It becomes a particle state by the gravitational field, due to the changes 
in the equation of state of the inflaton [196]. Spin 1/2 particles, such as the electrons, would 
be conformally invariant at the high energies prevailing during inflation, and consequently 
are not created in large numbers. Therefore, we must seek for a minimally coupled charged 
scalar field, of which none is included in the standard model but only in its supersymmetric 
extensions [197]. 

The scalar field can be decomposed into its real and imaginary parts as $ = (0i + 02) /\/2, 
and the associated current as 

J a = J? + J%, (5.1.15) 

where 

J\k = ie OA02 - </>2<9 fc 0i) (5.1.16) 

and 

hk = ~e 2 A k U\ + <j>l) . (5.1.17) 



For a crude estimate of the field we can neglect Jik- Assuming that Ohm's law holds, J\ - 
<; (77) E % , where <; (77) is the electric conductivity. Then, the equation of the magnetic field is 



8 2 + a (77) <r (77) d v - V 2 j B a = (V x J ± J . (5.1.18) 

Because the current is stochastic, the induced field must be evaluated through its two-point 
correlation function, 



(B t (X a ) B V ( X ' a )) = H^ f ^l D ret (^ ^ jjret ^ yla) ^ {y a ? y , a) ? (fi ± 10) 

J (2tt) 
where the "noise kernel" Nn> is the Hadamard two-point function 

No, (y a ,y ,a ) = ({(V x J 1 ). (x a ) , (V x Jj) v (x' a )}) , (5.1.20) 

and D ret (x a , y a ) is the retarded propagator of eq. (5.1.18). We are interested in fields coherent 
over a scale A, so the spatial integral in eq. (5.1.19) must be weighed by a window function, 
W (A), that filters out frequencies higher than the one associated to the field's scale of coher- 
ence, A -1 . After weighing, the magnetic energy density stored today in a region of size k~ x 
can be directly inferred, giving 



1 tod c^ m l 



2 ^V^( 1 + ri/2 ) +e # + 2 ^( 1 + ^ 



(5.1.21) 



where tuq is the bare mass of the scalar field, Trh the reheating temperature of the universe 
(i.e., the temperature at the beginning of radiation dominance), % = e 2 Tjm the electric 
conductivity at the beginning of radiation dominance, T1/2 the mean lifetime of the current 
and H the Hubble parameter, treated as time independent during inflation. Observe that the 
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field intensity depends very weakly on T1/2. Assuming instantaneous reheating, Trh = ^/HMpi 
an estimate of (B\) tod on a comoving galactic scale X ga i ~ 1 Mpc is 

(B x ) tod ~ e 3 — ^ ~ IO- 51 Gauss (5.1.22) 

which is about 15 orders of magnitude weaker than the minimum required to feed the galactic 
dynamo. 

Calzetta et al have considered the effect of the "London current", eq. (5.1.17) [198]. In this 
case the evolution equation of the magnetic two-point function shows two kernels, a local 
and a non-local one. Of these, the local (non-dissipative) one dominates over the non-local 
(dissipative) one by several orders of magnitude throughout reheating, which means that 
dissipation in this system is not due to ordinary electric conductivity. 

The equations for the magnetic field can be recast in the form a Langevin equation, which 
due to the local kernel looks like the London equation for a superconducting medium: 



£ + C^[,(T) 



B s (r) = F 2 [z (t)] . (5.1.23) 



Here F (z) = z 1 ^ ' J 3/27 (z), C 2 = e 2 z~ 4/ \ with e 2 = e^/^T 2 [(2 7 + 3) /2 7 ] In (A/T), 7 be- 
ing a parameter that determines the temperature evolution during reheating, A = g 1 ' 2 TnH / H , 
(with < g < 1 being a coupling constant of the mass to the thermal bath and Trh being the 
reheating temperature), and T is the dimensionless wavenumber corresponding to the origi- 
nal inflationary patch. Finally, J^/2-y (z) is a spherical Bessel function and T 2 (■ ■ ■) a Gamma 
function. 

Due to this current, the heavily amplified long- wavelength modes of the scalar field act as 
a Landau-Ginzburg order parameter in a superconductor, and as in the Meissner effect, the 
photon acquires a time- dependent mass. This allows for an exponential growth of the Maxwell 
field during reheating. The obtained intensities, however, are too weak (~ 10~ 53 G) to seed 
the galactic dynamo. Besides, in this model the amplification is very sensitive to the details of 
the reheating scenario, so it is not possible to obtain generally valid estimates for the resulting 
magnetic intensity. 



5.1.3 Primordial magnetic fields from metric perturbations 

Amplification of electromagnetic vacuum fluctuations can also be achieved by scalar pertur- 
bations in the metric during the transition inflation- reheating, i.e., by breaking the conformal 
flatness of the background geometry, rather than the conformal invariance of the Maxwell field 
equations [193,199]. The main effect is due to super-horizon scalar perturbations, specially by 
those modes that reenter the horizon right after the end of inflation. These fluctuations create 
an inhomogeneous background, in which the magnetic field evolves in a non conformally in- 
variant way: the mode-mode coupling between electromagnetic and metric fluctuations mixes 
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positive and negative frequency modes of the former field, thus breaking its conformal invari- 
ance. 

The line element of a flat FLRW model with scalar metric perturbations, in the longitudinal 
gauge, reads [200,199] 

dS 2 = a 2 (n) [- (1 + 2$) drf + (1 - 2$) dxrfx*] , (5.1.24) 

with $ (77, x l ) representing the gauge invariant gravitational potential. To first order in the 
cosmological and e.m. perturbations, 5Ai, the evolution equation of the Fourier transform of 
the latter is [199] 

j^SAi (k, rj) + k 2 5A % (k, n) - J t (k, rj)=0, (5.1.25) 

with Jj (k,r)) being a source term that depends only on the Fourier transform of the metric 
perturbations and on their time derivatives [199] . The resulting field strength depends on the 
power spectrum for super-Hubble metric perturbations, which is given by [199] 

(7 \ ra-1 
—J , (5.1.26) 

where As — 5.10" 5 sets the normalization at the COBE scale (Ac — 3000 Mpc). At decoupling, 
the field strength on a coherence scale corresponding to a galaxy, kc, is 

where k max is a cut-off that must be introduced in the case of a blue spectrum (n > 1) to 
avoid excessive primordial black hole production, and that for negative tilt, is related to the 
minimum size of the horizon (k max < aiHi, I denoting the end of inflation). Observe that 
eq. (5.1.27) is a function of k max , i.e., of the mechanism that generated the perturbations. 
The resulting magnetic field spectrum is thermal (B^ ~ /c 3//2 ) in the low-momentum tail. The 
relation between the energy densities in magnetic field and photons, for a suitable wavenumber 
kc, at decoupling turns out to be 

^ (k G ) ~ 1.4 x 10 36 ( -^-Y . (5.1.28) 

Pb \ k„ 



h max , 



The obtained intensities are upper limits, as dissipative effects were not taken into account 
when deriving expression (5.1.27). 

Scalar metric perturbations can grow exponentially during preheating [201,202], thus inducing 
strong enhancement of magnetic fields. Let us consider the Lagrangian [193] 

L = -^ - - A F ab F ab - ^dy - V (0) - \d aX d a X ~ ^ W , (5.1.29) 

i07rG 4 2 11 
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with V ((f)) = (1/4) A0 4 , being the inflaton and x the scalar field it is coupled to. In this 
case metric perturbations are expected to grow due to enhancement of the scalar field pertur- 
bations, and in turn the former would stimulate the growth of magnetic field perturbations 
through gravitational scattering. 

Assuming that on super-Hubble scales $ depends only on time, and adopting the Coulomb 
gauge (A = 0, d l Ai = 0), the Fourier component Ai (k) of the gauge potential reads 

A!! (k, rj) + k 2 A, t (k, r,) = 2&A[ (k, rj) . (5.1.30) 

Defining A t (k,rj) = (1 — <&) Ai (k,rj), to eliminate the term in A\ (k,rj), eq. (5.1.30) recasts 
into 

A!! (k, n) + k 2 A t (k, r]) = &% (k, ri) , (5.1.31) 

with its solution given in integral form by 



Ai (k, ri) = Af ) (k, v) + j T &'A (k, rf) sin [k (rj - rj')] drj' . (5.1.32) 

k Jna 



Decomposing the scalar fields as tp — > (p + Stp (where tp denotes either <fi or x), the evolution 
equation for the Fourier transformed metric perturbations reads 

$ (k, n) + H$ (k, ri) = 4vrG [<p (k, rj) 5<p (k, rj) + x (k, v) $X (k, r})] . (5.1.33) 

When the fluctuations Sx (k,rj), with low k, are excited during preheating, the corresponding 
metric and inflaton perturbations, $ (k, rj) and 5(f) (k, rj) respectively, grow on large scales and 
thus enhance the magnetic fluctuations [193]. However, with the increase of g 2 /X, the long 
wavelength modes Sx (k, rj) are suppressed during inflation. Sub-Hubble fluctuations, on the 
contrary, do not suffer from suppression and exhibit parametric amplification during reheating 
[203]. Therefore, the mode-mode coupling between small-scale metric perturbations and large- 
scale magnetic fields in eq. (5.1.30) can enhance the latter. This model, however, has many 
uncertainties and complexities, which require further research because they make it difficult 
to obtain reliable estimates for the final magnetic intensities. 



5.2 Magnetogenesis in phase transitions 



The actual state of the particles in our universe is the result of phase transitions that occured 
in the early phases of the expansion. At least two phase transitions are believed to have 
taken place in that epoch: the electroweak (EW - at Tew ~ 100 GeV) and the quantum 
chromodynamical (QCD - at Tqqd ~ 200 MeV). In the former case, the transition is from a 
symmetric, high temperature phase with massless gauge bosons to the Higgs phase, in which 
the SU (2) x U (1) gauge symmetry is spontaneously broken and all the masses of the model 
are generated. 

In the QCD case, the transition is from a quark-gluon plasma to a confinement phase with 
no free quarks and gluons. At the same energy scale, it is expected that the global chiral 
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symmetry of QCD with massless fermions is spontaneously broken by the formation of a 
quark pair condensate. 

First order phase transitions occur via bubble nucleation. Domains of new phase of broken 
symmetry form, whose sizes are at most of the order of the horizon at that time. As the 
horizons grow, different domains come into causal contact and bubble walls collide with each 
other. Magnetogenesis occurs through violent processes that take place during these colli- 
sions: reconnection of magnetic field lines carried by the walls, MHD dynamos induced by 
the turbulent flows produced by the collisions, etc. In every case, the question is whether the 
generated fields can explain the intensity and morphology of the observed magnetic fields, 
or to seed further amplification mechanisms, such as turbulent dynamos that could operate 
within galaxies. In general, the magnetic fields that are produced during phase transitions can 
be very strong, but typically have very small coherence lengths. 

Second order transitions occur in a smooth and regular way, with approximate local thermal 
equilibrium being maintained throughout the process. In spite of this, magnetogenesis can 
also be possible as shown below. 

It was recently proved [36] that the QCD transition in the hot universe was an analytic cross- 
over rather than a phase transition. In this sense, the results on magnetogenesis obtained 
by considering the QCD transition as first order are invalid and new research needs to be 
done. We shall therefore review in this section only the EW phase transition, which also 
seems to provide a very suitable scenario for magnetogensis, since it facilitates the separation 
between electric and magnetic fields as classical fields. Besides, while the Standard Model 
predicts a smooth cross-over for this transition, its extensions can give a strong first-order 
phase transition, which is a fundamental ingredient for electroweak baryogenesis and the 



generation of primordial magnetic fields] 19 1 Supersymmetric extensions of the Standard Model 
have been the most intensively studied [208,209,210,211,212,213], but it is also possible to get 
a strong transition from more generic two-Higgs doublet models [214,215], from technicolor 
theories [216], etc. [8-17] [217,218,219,220]. 

Phase transitions in the early universe lead to another class of mechanism for generating 
primordial magnetic fields, based on the Kibble mechanism [221], i.e., on the generation of 
cosmic strings. If the vacuum manifold Ai of the broken gauge theory that exhibits a phase 
transition has a nontrivial first homotopy group [222,223], then a cosmic string network will 
form generically [221]. This network has a characteristic length scale £ (£), which expands with 
the expansion of the universe. Infinitely long strings and loops are formed, the smallest of the 
latter decaying away via gravitational radiation. The result is a scaling solution, in which 
the string properties such as £ (£) are proportional to the time passed [224,225,226,227,228]. 
This means that, if cosmic strings can produce randomly oriented magnetic fields, these could 



19 In the Standard Model of particle physics, the electroweak phase transition is of first order if the 
Higgs boson mass has mg < 72 GeV [204,205,206]. Also, to suppress sphaleron processes in the 
broken phase, would actually require that mn < 35 GeV. However, the current experimental limit 
of rriH is well above these values, at m# > 114 GeV [207], thus turning the Standard Model an 
inadequate theory for baryogenesis. 
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be coherent by the Vachaspati mechanism [229] over galactic scales at the time of galaxy 
formation, as required by the dynamo paradigm. In the last subsection we review some works 
done on this mechanism. 

Recently, a new mechanism of early magnetogenesis was proposed by Dolgov et al [230], 
whereby ferromagnetic domains of condensed W bosons would form in the broken phase of 
the standard electroweak theory. These domains could create large-scale magnetic fields that 
would survive after the decay of the domains due to flux-freezing. Although the authors do not 
give estimates for the produced fields, their work points towards a new direction that should 
be explored further. 



5.2.1 Magnetogenesis in the electroweak phase transition 

In his seminal work of 1983 [231], C. Hogan was the first to propose a mechanism of magne- 
togenesis based on a small-scale dynamo induced by a first order phase transition in the early 
universe. His aim, however, was not to explain the fields observed in galaxies, but to study 
the effect of the induced fields on structure formation. The dynamo that Hogan proposed 
would be induced in the wall of the bubbles by the ordered release of free energy during the 
transition. Each bubble would be an independent dynamo, producing fields correlated only 
on the scale of the bubble radius. The result is a pattern of randomly oriented field lines 
that, properly averaged, would produce a large scale field spanning over regions that are not 
causally connected, i.e. coherent on larger scales, whose spectrum is of the form B\ oc Z~ 3 / 2 
i.e., a dipole field, but of weaker intensity. 

Since then, several mechanisms for magnetic field generation by first order phase transitions 
have been proposed in the literature. In 1995 Baym, Bodeker and McLerran [232] proposed 
an also dynamo-based mechanism, whereby seed fields (produced by thermal fluctuations) are 
amplified by a turbulent dynamo induced by the collision of supersonic shock waves created 
by the expansion of the walls of the broken symmetry bubbles. Their work was framed within 
the standard model, since at the time it was still believed that a first order phase transition 
was possible in it. In this sense, the resulting magnetic intensities would be incorrect. 

Concretely, when the expansion of the universe supercooled the cosmos below a critical tem- 
perature, of the order T cr ~ 100 GeV, then (at random locations) the Higgs fields tunnels 
from the unbroken SU (2) x U (l) y phase to the broken U (l) em phase, forming bubbles that 
expand and convert the false vacuum energy into kinetic energy. As the shock fronts collide, 
turbulence is developed in the cones associated to the bubble intersection, with Reynolds 
number 

i?e~^, (5.2.1) 

A 

where v w is the wall velocity, Rb is the size of a bubble at collision time and A is the scattering 

length of fluctuations in the electroweak plasma. Baym et al found that for scales smaller 

than R b , Re ~ 10 12 . Assuming v w ~ f/^id ~ lO^ 1 , that the typical bubble radius after the 

completion of the phase transition is R b ~ f^H^ly (with f b ~ 10 -2 — 10~ 3 being the fractional 

size and H"j^ ~ Mpig* T~ 2 ~ 10 cm that the event horizon at the electroweak scale - Mpi 
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is the Planck mass and g* ~ 10 2 is the number of massless degrees of freedom in the matter) 
and that A ~ Tg EW a 2 | log en | (with a the fine structure constant and qew ~ 9* the number 
of degrees of freedom that scatter by EW processes), Baym et al obtained 

Re ~ 10~ 3 — "a 2 |log a\ ~ 10 12 . (5.2.2) 

J- c 

Such a large Reynolds number means that turbulence is fully developed on scales smaller 
than Rb. Assuming that the electric conductivity is large at that epoch [49], strong magnetic 
turbulence should exist and in that situation kinetic and magnetic energies are in equipartition, 
allowing us to estimate that 

B 2 (R b ) ~ g,Ty flmd . (5.2.3) 

To obtain the intensity of the large-scale field, Baym et al assumed that the small-scale field 
formed a pattern of continuously distributed dipoles, with distribution being a Gaussian. 
Therefore the correlation function of the dipole density is 



y(0)i/>(r)) = «<W(r) , (5.2.4) 

while the one of the magnetic for r ^> f^H^r reads 

iT(0)B*(r))^bg(^). (5.2.5) 

From eqs. (5.2.3) and (5.2.5), Baym et al. obtained that 

where (• • • ) R means averaging on a scale R. The present-time estimate for this magnetic field 
on a galactic scale, l gal ~ 10 9 AU, is B (l gal ) ~ 1(T 17 - 1(T 20 . 

In 1991 T. Vachaspati [229] proposed a mechanism of magnetogenesis based on second order 
cosmological phase transition. These transitions would produce domains of different vacuum 
expectation values for the Higgs field, with these differences amounting to gradients in the 
field. The latter would ultimately lead to electromagnetic fields after the completion of the 
transition. This mechanism can produce fields associated to other (unbroken) symmetries (like 
SU(3)) as well. When applied to the electroweak transition, and assuming that the initial 
correlation scale, \% ~ 2 {gTi)~ , is of the order of the inverse mass of the W boson, with 
Ti ~ 10 2 GeV, an initial intensity of B ~ gT 2 /2 ~ 10 23 G is obtained for that correlation 
length. The initial energy density of the field is comparable to that of the universe, Qb (ti)- 
For a region of size £i = Nxi, with N ^> 1, the Higgs field is randomly oriented. Consequently 
the initial magnetic intensity on that scale would be B^ ~ gT 2 /4N, which at the electroweak 
scale gives B N (t EW ) ~ lO 23 ^" 1 G. At the QCD scale, we finds B N (t QCD ) ~ lO 18 ]^ 1 G and 
today Bn (t to day) ~ lO^iV -1 G (with iV > 10 13 in all cases). For a scale of 100 Kpc today, 
N = 10 24 and thus B ~ lO" 30 G. 
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The work of Vachaspati was questioned by Davidson in Ref . [233] . She computed the electric 
current due to the dynamics of the Higgs field and showed that it vanishes during the EW 
phase transition. Her conclusion was that no large-scale magnetic fields are generated by the 
classical rolling of the Higgs vacuum expectation value during the electroweak phase transition. 
Later, in 1998, Grasso and Riotto [234] reanalysed the generation of magnetic fields during 
the EW phase transition and found that the Vachaspati mechanism was plausible. Grasso and 
Riotto analyzed the two possibilities for the phase transition: first order and second order. 
They showed that the magnetic induction is connected to some semiclassical configurations 
of the gauge fields, such as electroweak Z-strings and H^-condensates. The initial Higgs field 
configuration is 

$ m (X) = -L exp Li -^ln a A ° , (5.2.7) 

V2 V 2 / \p(X)e 1 ^ 2 ) 

with r a representing the Pauli matrices, n a a unit vector in the SU (2) isospace, 9 (X) the 
U(l) Higgs field phase and p (X) the modulus of the Higgs field. The equation of motion for 
the SU (2) gauge field in the adjoint representation is 

D C F% = -g \p\ 2 d h 9 (X) (n a - n c c a ) , (5.2.8) 

where it is assumed that the initial gauge fields W® and their derivatives are zero at t — 0. 
Also, = $t r a$/$t$ = cos #0 O + s i n ft x O + 2 sin 2 (9/2) (n • O ) n, with 4>l = - (0, 0, 1). 
As n does not depend on the space coordinates, it is always possible to assume that it is 
perpendicular to 0o- I n other words, can be always obtained by rotating 0o at an angle 9 in 
the (0o, ft-) -plane. Then, eq. (5.2.8) becomes 

D C F% = -g \p\ 2 d b 9 (X) n a . (5.2.9) 

This clearly shows that only the gauge field component along n, namely A a = n b W a i, is 
created by a nonvanishing gradient of the phase between the two domains. When the full 
SU (2) x U (1) Y group is considered, it is no longer possible to choose n arbitrarily, because 
the different orientations of n, with respect to O , correspond to different physical situations. 
Setting n parallel to O and assuming that the charged gauge field does not evolve significantly, 
Grasso and Riotto found the following complete set of evolution equations, which is valid for 
a finite (though short) time after the bubbles first contact: 



d a F* b = ^-P 2 (X) (dw + tt-^-Z,) , (5.2.10) 

2costV V 2 cos 9 W J 



and 



d a d a (p (X) e^ 2 ) + 2A (> (X) - i? 2 ) p (X) e^' 2 = . (5.2.11) 

Here d a = d a + ^ 2 cos6> ^ a , with rj being the vacuum expectation value of $ and A the quartic 
coupling. Expressions (5.2.10) and (5.2.11) are the Nielsen-Olesen equations of motion [235]. 
Their solution describes a Z- vortex with p = at its core [236,237]. The geometry of the system 
implies that the vortex is closed, forming a ring whose axis coincides with the conjunction of 
the bubble centers. 
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To determine the magnetic field produced during the process described above, it is necessary 
to give a gauge-invariant definition of the electromagnetic field in the presence of a non-trivial 
Higgs background. Grasso and Riotto chose 

d a FX = 2t^e w d a [Z a d b \ogp{X)-Z b d a \ogp{X)} , (5.2.12) 

while Tornkvist used [238] 

jpem = _ sin w ^ F C b + cos wF Y + !^_ £cd£ ^ ( D jy (p^y ^ ( 5 _ 2 _ 13 ) 

obtaining no electric current and so no magnetic field. It must be stressed that both definitions 
(5.2.12) and (5.2.13) fulfill the requirement that they reproduce the standard definition in the 
broker phase with an uniform Higgs background. Grasso and Riotto noted that the presence 
of an inhomogeneous PU-condensate, produced by string decay, gives rise to electric currents 
that can sustain magnetic fields even after the Z-string has disappeared. 

An attempt to predict the strength of the magnetic field at the end of the EW phase transition 
was done by Ahonen and Enqvist [239] and by Enqvist [240], who analyzed the formation of 
ring-like magnetic fields in collisions of bubles of broken phase in an Abelian Higgs model. 
Under the assumption that magnetic fields are induced by a process similar to the Kibble 
and Vilenkin mechanism [241], it was concluded that a field of the order of B ~ 2 x 10 20 G, 
with a coherence length of about 10 2 GeV -1 , could be induced. In addition, assuming that the 
plasma was endowed with MHD turbulence, Ahonen and Enqvist found that the coherence 
scale could be enhanced by the inverse cascade of the magnetic helicity, and so a field of 
B rms ~ 10 -21 G on a comoving scale of 10 Mpc could be present today. As stated earlier, 
however, the problem with first order phase transitions in the standard model is that they are 
incompatible with the experimental lower limit for the Higgs mass. 

Grasso and Riotto also analyzed the creation of magnetic fields when the EW phase tran- 
sition is of second order. In this case domains where the Higgs field is physically correlated 
are formed near the critical temperature. The formation of topological and non-topological 
vortices, is a common phenomenon in second order phase transitions via the Kibble mecha- 
nism. It is also known that the non-topological vortices share many common features with 
the electroweak strings [242] . In this sense, Grasso and Riotto argued that electroweak strings 
are formed during the second order EW phase transition. To estimate the density of vortices 
(and consequently the mean magnetic field), it is necessary to know the Ginzburg tempera- 
ture, Tq- This sets the threshold at which the thermal fluctuations of the Higgs field, inside a 
given domain of broken symmetry, are no longer able to restore the symmetry. The Ginzburg 
temperature was computed by the authors of Ref. [234], after comparing the expansion rate 
of the Universe with the nucleation rate per unit volume of sub-critical bubbles of symmetric 
phase with size equal to the correlation length of the broken phase. The latter is given by 



, _1 

ub „a 



e- 5 a/ T , (5.2.14) 



where i b is the correlation length in the broken phase and S^ is the high temperature limit 
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of the Euclidean action [243]. For the EW phase transition, Tq — To, and the corresponding 
size of a broken phase domain is determined by the correlation length at T = Tq, i.e., 

4(T G )- 2 = y"((0(T G )),T G ), (5.2.15) 

where V (0, T) is the effective Higgs potential. Using the fact that 4 (Tq) depends weakly on 
Mil-, Grasso and Riotto estimated the magnetic field strength, on a scale £b (T G ) at the end 
of the EW phase transition, to be B^ ~ 4e _1 sin 2 Q-w^\ (^g) ~ 10 22 G. To obtain the intensity 
on cosmologically interesting scales, the authors of Ref. [234] followed the procedure of line 
averaging sugested by Enqvist and Olesen, i.e., {B) rms L = B^/y/N, where N is the number of 
domains crossed by line, obtaining that a field coherent on a scale of 1 Mpc today would have 
an intensity of Bq (1 Mpc) ~ 1CT 21 G. It must be pointed out, however, that all these studies 
do not take into account the dissipative effects of the primordial plasma. Consequently, the 
corresponding numerical results should be treated as upper limits. 

The mechanism proposed by Vachaspati [229] and later analyzed by Grasso and Riotto [234] 
(see also Cornwall [32]) was recently numerically confirmed and improved by Diaz-Gil et al. 
[244,245]. The authors considered the full SU (2) ® U (1) model in the framework of hybrid 
inflation. After a short period of hybrid inflation that ends at the EW scale, where non- 
linearities in the Higgs and gauge fields can be neglected, tachyonic preheating develops and 
non-linearities in the fields cannot be neglected anymore. During this period the Vachaspati 
mechanism operates, and magnetic string-like configurations appear due to the gradients in 
the orientation of the Higgs field. The important feature of the induced magnetic fields is that 
they are helical, i.e., they posses a non null r.m.s. magnetic helicity. During the subsequent 
phase of (first order) EW symmetry breaking, the magnetic fields are squeezed in string-like 
structures localized in the regions between bubbles, where the gradients of Higgs fields are still 
large. The evolution of the coherence scale of these fields can be tracked for a short period of 
time after the end of the phase transition. At that time it is important to track the evolution 
of the low momentum part of the spectrum, which is the one that can seed the fields for 
galaxies and clusters of galaxies. It is seen that it carries a fraction of ~ 1CT 2 of the total 
energy density, which would be enough to explain the magnetic fields observed in clusters. 
The correlation length grows as fast as the particle horizon (i.e., linearly in time) and this 
behaviour is interpreted as an indication that an inverse cascade of magnetic helicity is in 
operation. However, it is not possible to extrapolate this behaviour to later times, due to our 
limited knowledge on the primordial plasma features. 

Stevens and Johnson [246,247] analyzed the possibility of magnetogenesis by a first order 
EW phase transition, possible for some choices of parameters in the minimal supersymmetric 
Standard Model (see also [248]). They considered the Lagrangian 

C = L\ + £2 + ( leptonic, quark and supersymmetric partner interactions) , (5.2.16) 

with 

A = -\w l ah W iab - \B ab B ab (5.2.17) 
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and 



£ 



2 



id a - 9 -r-W a - 9 -BA $ 



2 



V($,T). (5.2.18) 



Also, T represents the temperature, while W l and -B^ are given by 

W l ab = d a Wt - d h Wi - ge) k WlW b k (5.2.19) 

and 

B ab = d a B b - d b B a , (5.2.20) 

respectively. In the previous equations, W l (with i = +, — ) are the W + , W~ fields, $ is 
the Higgs field and r l is the SU (2) generator (fermions are not considered in this model). 
In the framework of the MSSM the bubbles that consist of a region of space filled by the 
Higgs field with a cloud of the other constituents of the MSSM in the broken phase. From this 
Lagrangian, one obtains the linearized equations of motion with O (3) symmetry, which are 
suitable to study collisions where the Higgs field is relatively unperturbed from its mean value 
within the collision volume. Stevens and Johnson [247] found that the coherent evolution of the 
charged W fields within the bubbles is the main source of the electric current that generates 
the magnetic field. In their model, fermions are taken into account as a background that 
provides dissipation through electric conductivity. They numerically integrated the equations 
of motion of the model, paying special attention to the role of the surface thickness of the 
bubbles, finding that the main sensitivity is due to the steepness of the bubble surface: the 
steeper the transition, the more enhanced the seed field becomes. Despite this, the authors 
of Ref. [247] did not attempt to give the present-day value of the generated magnetic field, 
because of uncertainties in the properties of the host plasma. 

5.2.2 Magnetogenesis from cosmic strings 

The interaction between cosmic strings and magnetic fields was first discussed in 1986 by 
Ostriker et al [249] , while their connection with primordial magnetogenesis was first suggested 
by Vachaspati in 1991 [229]. Later the mechanism was further developed by Brandenberger et 
al [250], also for superconducting cosmic strings, who showed that these models are severely 
constrained by cosmological arguments: the only stable confirgurations for those strings are 
springs and vortons, which produce matter overdensities in the same manner primordial mag- 
netic monopoles do. So, these models had to be ruled out. 

In 1999 Brandenberger and Zhang [251] studied magnetogenesis by anomalous global strings 
and discussed for the first time the importance of the coherence length in these models. The 
authors proposed a mechanism based on the realization that anomalous global strings couple to 
electromagnetism [252] through an induced F^F^ term in the low-energy effective Lagrangian 
and therefore magnetic fields can be generated. The major advantage of the mechanism is 
that the coherence scale of the induced magnetic fields is basically the curvature radius of 
the inducing string. The mechanism is realized within QCD, namely there exists a class of 
stringlike classical solutions of the linear sigma model, that describes strong interactions below 
the confinement scale, called pion strings [253]. At low temperatures those strings are not 
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topologically stable, decaying at a temperature Td ~ 1 MeV, but within the plasma they can 
stabilize because the plasma interactions break the degeneracy among the three pions. 

Since a pion string is made of a and n° fields, it is neutral under the U em (1) symmetry. 
However, the 7r° couples to photons via the Adler-Bell-Jackiw anomaly. In the linear sigma 
model, the effective coupling of 7T° to photons is obtained from the contribution of the quark 
triangle diagram [254]. At low energies only pions and photons are important, hence the 
effective Lagrangian to leading order reads 

L low = f -fTr (d a Y)d a Y) - - A F ab F ab ~^j e abcd F ab F cd , (5.2.21) 

where N c = 3, £ = exp (if ■ 7t/f n ), f are the Pauli matrices and a is the electromagnetic 
fine structure constant. From this Lagrangian one also obtains the classical equation of the 
electromagnetic field 

d a F ab = --d a (^-)p ab . (5.2.22) 



The key effect is due to the anomaly term in eq. (5.2.22). Charged zero modes on the string 

will induce a magnetic field circling the string that falls off less rapidly, as a function of the 

distance from the string, than it is classically expected. Zero mode currents are automatically 

set up by the analog of the Kibble mechanism [221] at the time of the phase transition and 

therefore magnetic fields coherent in a scale of the string size are automatically generated. 

The coherent magnetic field, as a function of the distance r from the string, can be expressed 

as 

en / r \ a l 1T 1 
B(r) = N c e ^(^) ± (5.2.23) 

2n \r J r 

with n being the number density of charge carriers on the string, r giving the width of the 

string and a < 1. By dimensional analysis, Brandenberger and Zhang obtained that at the 

time t c = t c , when the strings form, r ~ T c _1 and n ~ T c . The initial correlation length of the 

string network, £ (T c ) increases rapidly, approaching a scaling solution of the form £ (t) ~ t. 

During this evolution the charge density is diluted as the strings stretch, while at the same 

time the merger of small strings into larger ones leads to an increase of charge. Assuming 

that the initial separation of the strings is microscopic, and that they decay during radiation 

dominance, Brandenberger and Zhang obtained 

n(t d )~(^Jn(t c ), (5.2.24) 

with p = 5/4 or 3/2 [251]. Also, the corresponding magnetic field at t d is 

B {td) ~ 1Q5 K^V 7 ™ 1 (I)" {rTcr/7T GaUSS • (5 ' 2 - 25) 

Brandenberger and Zhang assumed that between td and the present time, to, the field prop- 
agates through a perfectly conducting plasma and found that today the magnetic intensity 
should be 

■ 14 T c „-i ( T d\ P (Tq\ . . a /ir 



B (to) ~ 10 !GeV r ^ \¥J \f d ) {rTc) GaUSS • (5 ' 2 - 26) 
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Note that To is the present time temperature, tk pc is the present distance from the original 
comoving location of the string expressed in Kpc and r the physical distance at T d . Considering 
t as an estimate of the string separation, T c ~ 1 GeV, T^ ~ 1 MeV and p — 1, Brandenberger 
and Zhang obtained 

B (t ) ~ KT 26 (rT c ) a/n Gauss , (5.2.27) 

arguing that, if r is of cosmo logical order and T^ 1 < 1, we can have (rT c ) Q/7r ^> 1. 

To analyze the coherence scale of the fields, the authors assumed that after T d the field lines 
are frozen in comoving coordinates, obtaining 

f (td) c = PUz (t d ) = /31(T 2 Kpc T d [MeV]" 1 , (5.2.28) 

where z (td) is the redshift at t d and /3 ~ 1 for scaling strings. On scales larger than £ (t d ) c the 
fields have random orientation, yielding an average of 

B = -j=B(to), (5.2.29) 



with 

2 

N 



ci[kpc]^- 2 10 4 T,[MeV] 2 , (5.2.30) 

where d is the scale the coherente field is calculated. The authors found that B ~ 10 _2 -B (to)- 
Obviously, if the resistivity of the host plasma is accounted for, the suppression will be larger. 

Recently, Gwyn et al [255] extended the mechanism to heterotic cosmic strings arising in M 
theory. Those strings, being stable, would produce even stronger fields. This work is reviewed 
in the following section. Another possible way that cosmic strings could produce primordial 
magnetic fields was proposed by Dimopoulos [256], by Dimopoulos and Davies [257] and by 
Battefeld et al [258]. In these scenarios, the magnetic fields are be induced by vortices produced 
by cosmic strings via the Harrison- Rees [152] effect. This mechanism, however, was recently 
strongly criticized by Hollenstein et al [259], who showed that the Harrison- Rees effect is quite 
inefficient in producing cosmologically intersting magnetic fields. 
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6 Magnetogenesis beyond the standard model 



In this section we are going to review different types of magnetic field generation mechanisms 
involving theories beyond the standard model. Primordial magnetic fields will be generated 
during inflation. As it was shown in § 4, on a non flat background perturbations in the elec- 
tromagnetic field can be efficiently amplified during inflation within the standard model, that 
is within the standard linear theory of electrodynamics. On the contrary, on a flat background 
the amplification during inflation is not sufficient in order to be cosmologically relevant. Fol- 
lowing [35] (see also § 4.3 earlier) the ratio r of the energy density in the magnetic field p-g 
over the energy density p 7 in the background radiation is introduced, thus 

r=— , (6.0.31) 

p 7 

where pe = f^- In the case of linear electrodynamics on a flat background the magnetic energy 
density decays as a -4 , where a is the scale factor. Hence the ratio r is a constant as the universe 
evolves. This is also true in the radiation dominated era when the universe is dominated by 
a highly conducting plasma. The interstellar magnetic field in our Galaxy is of the order of 
a few pG. Assuming that a galactic dynamo, contributing an exponential factor in time, is 
operating since the time of the formation of the Galaxy requires an initial seed magnetic field 
at the time of galaxy formation of at least B s ~ ICT 20 G [260] which corresponds to a minimum 
magnetic to photon energy density ratio r given by r ~ I0~ 37 . There is some controversy about 
the efficiency of such a galactic dynamo, thus working under the hypothesis that there is no 
efficient amplification of the initial magnetic seed field due to a galactic dynamo but only the 
amplification due to the collapsing protogalactic cloud requires r to be at least of the order 
of r = 10~ 8 [35]. These lower bounds on r were derived assuming no cosmological constant. 
In a flat universe with a large positive cosmological constant assuming a galactic dynamo 
operating these bounds can be lowered significantly. In particular for reasonable cosmological 
parameters an initial seed magnetic field of at least B s ~ I0~ 30 G is enough to explain the 
present day galactic magnetic field strength [30]. This corresponds to r = 10~ 57 . In typical 
inflationary scenarios on a galactic scale of 1 Mpc r ~ I0~ 104 at the beginning of the radiation 
dominated era (cf. equation (4.3.9)) which is much below the required minimal value even 
in the presence of a cosmological constant. Therefore, in the case of a flat background, it is 
necessary to go beyond the standard model. There are different possibilities of modifying the 
standard four dimensional electromagnetic Lagrangian. 

Currently models of modified gravity enjoy an intense activity due to the fact that they can 
be used to describe the late time evolution of the universe at a global scale as well as, say, the 
observed rotation curves of galaxies. Thus these models combine the effects of dark energy, 
which is used to model the accelerated expansion of the present universe, and dark matter, 
which is postulated to exist in the form of halos around most galaxies. 

The gravitational sector in theories of modified gravity is usually described by a Lagrangian 
of the form (e.g., [261]) 



58 




S = — d*xy/=jf(R), (6.0.32) 



where f(R) is a function of the Ricci scalar R, most often chosen to be of the form f(R) — 
R + aR n , where a and n are constants. It should also be noted that one of the original 
realizations of inflation is given by f(R) = R + R 2 which can be shown to be equivalent to a 
conformally coupled scalar field [262]. Usually it is argued that modified gravity theories are 
some kind of effective description resulting from taking into account quantum corrections to 
the classical Einstein-Hilbert action. 

In order to study the generation of primordial magnetic fields in this type of theories the 
electromagnetic field has to be included in the Lagrangian [263]. Considering flat space the 
conformal invariance of the Maxwell Lagrangian in four dimensions has to be broken in order 
to generate magnetic fields strong enough to seed the galactic magnetic field. In the following 
a survey of models will be given which are used in this respect. 



6. 1 Gravitational coupling of the gauge potential 



Models involving the gravitational coupling of the gauge potential are described by La- 
grangians of the form RA m A m and R mn A m A n where R mn is the Ricci tensor. Gauge invariance 
of the Lagrangian is broken explicitly and for that matter it does not seem very appealing. The 
term RA m A m describes a massive photon with its mass given by m 7 ~ R l l 2 . Electro magnetism 
is then described by Proca theory and A m is the Proca field. 

The strongest bound on the photon mass within our galaxy is obtained by assuming a Proca 
regime on all scales. The Proca field contributes to the magnetic pressure of the intergalactic 
medium which has to be counterbalanced by the thermal pressure of the plasma. Observations 
assuming standard electrodynamics indicate that within our galaxy the interstellar medium 
is approximately in equipartition. From this it can be concluded that the magnetic pressure 
due to the Proca field has to be subleading with respect to the standard magnetic pressure 
[264]. This implies the bound m 7 < 10~ 26 eV [265]. Furthermore, in [265] it is pointed out 
that these limits depend on the mechanism on how the photon aquires mass. If it is via the 
Higgs mechanism then it is possible that large scale magnetic fields are effectively described 
by Maxwell's theory. In this case the strongest bound comes from the validity of Coulomb's 
law and is given by m 1 < 10~ 14 eV [266,265]. Using this type of Lagrangian on cosmologically 
scales, the typical scale is given by the value of the Hubble parameter today, Hq, which 
results in the estimate, m 1 ~ R? and R^ ~ H then at present the photon mass is given by 
m 7 ~ Hq ~ lCT 33 eV which is well below the above mentioned present limits on the photon 
mass [35]. 

The idea is that the initial magnetic seed field is created from the amplification of perturbations 
in the electromagnetic field during inflation. In [35] the resulting magnetic field at the end 
of inflation in this type of theories was calculated using the assumption that at the time of 
horizon crossing during inflation the energy density in this mode is determined by the Gibbons- 
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Hawking temperature. This was critically reconsidered in [267] were it was found calculating 
the spectral energy density from first principle quantizing the corresponding canonical field 
that this assumption actually is an over-estimation of the actual energy density. However, 
here we follow the original calculation of [35] . The equations of motion are derived from the 
Lagrangian [35 . 
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1 be 

*- = — T-Tmn-T — ntiA — —R mn A A 



(6.1.1) 



where A 2 = A m A m . Furthermore b and c are constants. The equations of motion lead together 
with the parametrization of the Maxwell tensor in terms of the electric and magnetic field, 
E a and B a in the "lab" frame, respectively, 



— E~ — E,. — E-. 



E x B z 



E y —B z 



-By 



E, B„ —B T 



(6.1.2) 



to the equations 



1 d ( 2 p \ id n A * n 
^■t— [a E a ) -cvj\B a ^^- = 

n A rin V / Tt Or 



a 2 drj 



(6.1.3) 
curlE a = 0, (6.1.4) 

where rj is conformal time. Thus the line element is given by ds 2 = a 2 {—drj 2 + dx 2 ) and 



cr or] 



n = 77 66 — h c 
\ a 



-+ - 
a \ay 



(6.1.5) 



where a dot indicates a derivative with respect to conformal time. Taking the curl of equation 
(6.1.3) and using it in equation (6.1.4) results in a wave-type equation for the magnetic field, 



a 2 drj 2 



a 2 B r 



n 



D z B a + —B, 



(6.1.6) 



Expanding in Fourier modes using F a (k,r)) = a 2 J d?xe B a (x,r]), yields 



20 Recall the notation convention: Latin indices take values between and 3. Greek indices take 
values between 1 and 3. 
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F a (k) + k 2 F a (k) + —F Q (k) = 0. (6.1.7) 

So that the additional terms in the Lagrangian can potentially act as pump terms, that is 
amplifying the spectral energy density of the magnetic field. The averaged magnetic field 
energy density is given by p mag (^) = ((B a B a )(x,ri))/(8ii). Using the correlation function 

(F a (k,v)F* (lv)) = I VKM) (<W - ^!r) , (6-1.8) 

results in 

PmagW = ^ J 'd 3 k\F,F»\. (6.1.9) 

Thus, with pft(k,ri) = k p ^ ag , the spectral energy density is given by 

MM) = -^ 3 |F„F"|. (6.1.10) 

Therefore, solving equation (6.1.7) for different types of scale factor, the magnetic energy 
density is estimated after the end of inflation using pe oc \F^F^\/a 4 [35]. Note that for standard 
electrodynamics, n — 0, and thus the magnetic energy density simply scales at the usual a~ A - 
rate of a frozen-in magnetic field. 

As can be seen from equation (6.1.7) modes well inside the horizon, that is those modes 
with comoving wave number \ki]\ ^> 1, simply oscillate, since the last term can be neglected 
and equation (6.1.7) reduces to the equation of a harmonic oscillator. In the opposite case, 
for modes well outside the horizon, satisfying \ki]\ <C 1 the second term in equation (6.1.7) 
becomes subleading and the resulting equation can easily be solved giving the solutions 



F^l (x n m± , (6.1.11) 



where m± = 1 1 1 ± a/1 — An ) . Calculating the behaviour of pp during the different stages of 
the universe: de Sitter (inflation), radiation and matter dominated epoch, the fastest growing 
modes are determined by p = m_ = 1 1 1 — a/1 — 486 — 12c I which is calculated during the de 

Sitter stage and q = m+ — |( 1 + \/l — 486 — 24c ) calculated in the radiation dominated era 

[35]. Thus, at the time of galaxy formation, the ratio of magnetic over photon energy density, 
r, is found to be [35] 

/ M \ 4(<?-p)/3 /TLtt\ 2(2<?-p)/3 / T \ - 8( ?/ 3 / A \ -2(p+2) 

r^(7xl0 25 )- 2 ^ +2 )f— V fi™V (-±) (—) , (6.1.12) 

V ; \MpJ \MpJ \M Pl ) VIMpcy ' V ; 
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where T* is the temperature at which plasma effects become important during reheating. 

1/2 / v l/3-i 

In [35] this is estimated to be of the order of T* ~ min 



T RH M 



T£hMpi 



) 



. Taking 

typical values for the physical parameters, there is a wide range for the exponents p and q 
such that r is bigger than the minimal required value in order to seed the galactic magnetic 
field with a galactic dynamo operating, r > 10~ 57 , or without, r > 10~ 8 . In figure 1 logr is 
shown for typical values of the parameters. 



M=10 17 GeV, r*=10 123 GeV, r RH =10 9 GeV 




-2.5 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 

q 

Fig. 1. The logarithm of the ratio r, magnetic field energy density over photon energy density, is 
shown for different values of p and q in a model including RA 2 and R^A 11 A v terms [35]. The numbers 
in the graph refer to the values of log 10 r along the closest contour line. 

6.2 Quantum corrections in QED in a curved background 



The QED one-loop vacuum polarization of the photon in a general curved background gives 
rise to terms coupling the Maxwell tensor to the curvature [268]. Vacuum polarization de- 
scribes the effect of virtual electron positron pair creation thus giving the photon a size of 
the order of the electron Compton wave length which interacts with curvature. This leads to 
an interesting space-time structure including the phenomenon of gravitational birefringence, 
where the photon propagation depends on its polarization and it can be faster than the speed 
of light [268,269,270,271]. In general the Lagrangian has the form [268] 



r — p p mn 



Am 2 



mn rplk 



bRF mn F mn + cR mn F mk F\ + dR mnlk F mn F lk + f(V m F mn )(V a F a n ) , (6.2.1 



where b, c , d and / are constants. The last term can be neglected with respect to the other 
terms since it leads to higher order derivative terms in the equations of motion [272] , 
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bRF mn + - ( R m Fi n — R n Fi r 



+ -J-2 (V a V a V b F 6n + R a n V b F ba 



0. 



dR 



ik 



iFik 



(6.2.2) 



In order to proceed here a different approach from the one used in the previous section will be 
followed. This will allow to actually determine the spectrum of the resulting magnetic field. 
Instead of assuming that the energy density of the magnetic field at the time of horizon cross- 
ing corresponds to the one calculated using the Gibbons-Hawking temperature, the spectrum 
of the resulting primordial field will be calculated by determining explicitly the Bogoliubov 
coefficients. This gives the particle production and thus the spectral energy density of the 
primordial magnetic field. The background model will be described by two stages, an infla- 
tionary stage during which the correction terms coupled to curvature are important and a 
radiation dominated stage determined by standard Maxwell electrodynamics. Since the re- 
sulting field strength will be compared with observational values at the order of galactic scale 
corresponding to 1 Mpc today, it is not necessary to include the evolution during the matter 
dominated period. Galactic length scales re-enter the horizon during the radiation dominated 
stage, which can seen from the temperature T 



t ~n 



A 



Mpc 



eV, 



(6.2.3) 



at which a scale for A < \ eq ~ 



U^h~ 2 



Mpc crosses back into the horizon. Thus a galactic 
scale of order 1 Mpc enters the horizon at a time when the universe was at about 78 eV and thus 
inside the radiation dominated era, since radiation-matter equality occurs at T eq 



5.6Q m h 2 eV. 



The background cosmology is chosen to be such that 



0(77) 



01 (£) V<Vi 
ai (=EgO V > Vi- 



(6.2.4) 



In the following a± = 1. The matching between the inflationary phase and the radiation 
dominated era takes place at rj = r]\. De Sitter inflation corresponds to fi = —1 and for 
/3 < — 1 power-law inflation is taking place. 

The Maxwell tensor is written in terms of the gauge potential A m , that is F mn = d m A n — d n A m . 
Furthermore, the radiation gauge, A = 0, d\A\ = will be used. Then the gauge potential 
is given in terms of the expansion in Fourier modes by 



Apfax) 



(2tt 



Ct rv ^ — ^ 



(A) 



M) 



A k (rj)e 



lk-X 



,(A)t 



Al(v)e 



-ik-x 



(6.2.5) 



A=l 
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where the sum is over the two physical polarization states and e^ are the polarization vec- 



,(A) 



,W 



,(A)t 



tors satisfying, e- • • k — 0. a- and a- are annihilation and creation operators, respec- 



tively, satisfying the standard commutation relations, [cb ,ctp 



: = [4 A)t ,ag' )+ ] and 
[ab , dp ] = 5\\i5^p. However, since the Lagrangian has additional terms coupling the elec- 
tromagnetic field to the curvature, the commutation relations between the operator Aj and 
its canonical momentum n a = m» C a \ are not the canonical ones. However, as will be done 
explicitly below, it is possible to define a canonical field which does satisfy the standard 
commutation relation. It is that field that is used to quantize the theory and calculate the 
production of particles. 



The mode functions satisfy 



F 1 (r ] )A k + F 2 ( V )A k + F 3 ( V )k 2 A k = 0, 
where a dot indicates 4- and 

dr/ 



(6.2.6) 



F 1 (ri)=l 
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^2 / V 



T]\m 2 \r]i 



F 3 (V) = 1 



-2(/3+l) 



m 2 e rjl \i]i 



-2/3-3 



/4$ / V 



-2(/3+l) 



rj{m 2 e \i]i 



1*1 = 



6b(/3 - 1) + c(/3 -2) -2d 



fx 2 = -2((3 + l)fi 1 



f*3 = 13 



6b((3 - 1) + c(2/3 - 1) + 2d(3 



(6.2.7) 



In the case where the additional terms in the lagrangian are absent, that isb = c = d = 0, 
the mode equation for A k reduces to a simple harmonic oscillator equation. In this case A k 
itself can be used to implement the standard quantization scheme. In general, however, it is 
necessary to use the canonical field ^/^(i],x) and its Fourier amplitude ty(r),k) defined by, 
respectively, 



#„(77,x) = F[ 2 (r])A^r],x) tf = F?(rj)A k . 

With this the mode equation for ^f(r], k) is given by 



(6.2.8) 



■$" + P^j = 0, 



(6.2.9) 



where a new dimensionless variable z = —krj has been defined and = ■$-. Moreover 



P 



K\Z 



-4/3-6 



1 + K 2 Z~ 2 ^ +1 ) 



1 K 3 Z- 2 ^ 1 + /t 4 ^~ 2(/3+1) 

2 1 + K 2 Z- 2 ^+V + 1 + K 2 Z-W+ l ) ' 



(6.2.10) 
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and 

«! = ^2«o K>2 = HlKQ k 3 = (2/3 + 3)/i 2 K K 4 = /i 3 K 

m„\-V^ 2< " +1) 



where Ko ^_,j |_j . (6.2.U) 

Here the maximally amplified (comoving) wavenumber k\ has been defined by k\ = -p-r. Fur- 
thermore if i is the value of the Hubble paramter at the beginning of the radiation dominated 
stage at r\\. It is related to k\ by k\ ~ H\. Thus the canonical field satisfies the equation of 
a harmonic oscillator. This is also the case of a free scalar field in flat space-time. Therefore 
the canonical quantization procedure will be applied to the canonical field \l/, which will be 
written as 



flu ,W 



S x )r , „wu* 



%(r,, x) = J d*h e^ [affr + o^7|j . (6.2.12) 

Hence it will be required that the mode functions f%(x) = tye lkx /(2tt)* and /| form an 
orthonormal set, that is satisfying [273], 

(/*/*) = 6i3) $- k'), Uhfg) = S^(k-k'), (4,/*,) = 0. (6.2.13) 
Furthermore, the scalar product is defined by 

(/* U) = ~i L hi*) dm f*Md^ m , (6.2.14) 



where <i£ m = n m c?S and n m is a future-directed unit vector orthogonal to the space-like 
hypersurface S which is taken to be a Cauchy surface. Moreover, c?S is the volume element 
in E. Also the notation f%(x) dm fg = fk^mfc, ~ {dmf%)fc, was used. Since \I> is quantized 
in flat space-time, dYT 1 = 5™d 3 x and the normalization condition on the mode functions f% 
reduces to 



V(ri, k)d^(r], k')* - [d„tf (77, k)] ^(77, k')* = i, (6.2.15) 

which is the Wronskian of the solutions of the differential equation (6.2.9). The field equation 
in real space for the Fourier transform ty(r),x), assuming it to be real, can be derived from 
the Lagrangian 

£* = \ (V^ A - m 2 c ^ 2 ) . (6.2.16) 

The effective mass can be determined by going back to /c-space. Using equations (6.2.9) and 
(6.2.10) 
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m 2 cS = k 2 (l-P). (6.2.17) 

It can be verified that ^^{t], x) and its canonical momentum satisfy the standard commutation 
relations. The time-dependent effective mass m e s(r]) reflects the dynamics of the cosmolog- 
ical background. It also indicates that there is no unique vacuum state. Having one set of 
orthonormal functions ft another orthonormal set of mode functions ft can be found. Then 
the canonical field ^ has the expansion in terms of the annihilation and creation operators 

4 A) ' 4 A)t and h 



*p(ri, x) = j d 3 k efl [a™ ft + 4 A)t /|] ( 6 - 2 - 18 ) 

and this defines a new vacuum state |0), 

4 A) |0) = 0, (6.2.19) 

for all k and A and a new Fock space. Since both sets of mode functions are complete, they 
are related by the Bogolubov transformation [273], 

where «£- and /?£- are the Bogolubov coefficients satisfying Z^C^Vj^-r — Pqkfl*h) = ^qr- 
Moreover, it is found that, suppressing the index A, [273], 



EK-^+%4)- ( 6 - 2 - 21 ) 



Thus, equation (6.2.21) implies that the vacuum state |0) is in general not annihilated by a%, 
but rather gives 

%H0)=E%IV)^0. (6.2.22) 

q 

This means that the expectation value of the number operator JVg = at a^ of /^-mode particles 
in the state |0) is given by 

<0|JV £ |0> = £ |%| 2 . (6.2.23) 

q 

In order to determine the particle production due to the time-dependent cosmological back- 
ground the mode functions are matched at the transition time r\ — r\\. Furthermore, on 
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subhorizon scales corresponding to z >> 1, the mode equation (6.2.9) reduces to the equation 
for a free harmonic oscillators and therefore does not give any important contribution. Only 
modes on superhorizon scales are relevant, since in that case, for z <C 1, the mode equation 
can be approximated by 



where 



(6.2.24) 
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6b(j3 - 1) + c(2/3 - 1) + 2d(3 
Qb(/3 - 1) + c(/3 - 2) - 2d ' 



(6.2.25) 



The particular choice (3 — — 1 describes de Sitter inflation, and in this case £i = and £2 = 1 
leading to a plane wave solution which was also noted in [268,35]. Furthermore, (5 = —2 implies 
£1 = 0, but £ 2 = igb+4c+2d • Equation (6.2.24) solved during power law inflation, /3 < — 1 and 
/3 7^ —2, results in the following solution in terms of the Hankel function of the second kind, 
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which gives the correctly normalized incoming wave function for rj — > —00 for £2 > . This 
means that the incoming vacuum solution at infinity is a plane wave solution and moreover 
approaches the positive frequency solution in Minkowski space-time. It is assumed that elec- 
trodynamics becomes standard Maxwell electrodynamics at the beginning of the radiation 
dominated stage at rj = rji. Thus the terms due to the interaction between curvature and the 
electromagnetic field in the mode equation (6.2.24) can be neglected which leads to a free 
harmonic oscillator equation which is solved by the superposition of plane waves, 
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(6.2.27) 



where Z\ = k\r)i\ and c± are the Bogoliubov coefficients, corresponding to «£- = c + 5^-, and 
pt = c_<5^^. Therefore to determine the magnetic field energy spectrum, the Bogolubov coef- 
ficients are calculated by matching the solutions of the gauge potential and its first derivative 
at 77 = 771 on superhorizon scales. Using the small argument limit of the Hankel functions 
[274], this leads to f3 ^ -§ [272] 
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(6.2.28) 



where it was used that in the approximation used here, F^rji) ~ /ii f ^ J .In the case f3 — — | 
the limiting behavior of the mode function on superhorizon scales leads to a a divergent factor 
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In 2 (y/^2T-) i n l c -| 2 - Thus, we will not pursue this case any further. Including both polarization 
states the total spectral energy density of the photons is given by (cf., e.g., [275]) 



Since the electric field decays rapidly due to the high conductivity of the radiation dominated 
universe, the spectral energy density (6.2.29) gives a measure of the magnetic field energy 

density, pe- Using the density parameter of radiation, f2 7 = f^j (^H, the ratio of magnetic 
over background radiation energy density r is given for j3 7^ —2, — |, — 1, by [272] 

where Mpi is the Planck mass. The magnetic field energy density can also be calculated 
using the two point function of the magnetic field, {B^(k)B*(k')). This leads to an expression 
similar to (6.2.30). Furthermore, the form of the magnetic field spectrum (6.2.30) imposes the 
constraint v < | . This implies the range for (3 given by —3 < f3 < — 1 taking into account the 
constraint from power law inflation. Using the constraint which was used to derive the mode 
equation (6.2.24) 



Mir) >1; (6 - 2 - 31) 



the maximal value of r which can be achieved within this model can be estimated. It is found 

3 
2' 

111 ( H x 



to be, for (3 ^ —2, — §, —1, at uq = 10 14 Hz corresponding to a galactic scale of 1 Mpc, and 



using the maximal amplified frequency evaluated today, wi(??o) = 6 x 10 11 \j^)' 2 Hz, [272] 



a-79+52*/ |W,,M 2 I ' ^ l&\ ( Hl 



r max ^G) = lO- 79+ ^[T(u)Y^--uj \^j \j£.J . (6.2.32) 

In figure (2) log 10 r max is shown for the case that the parameters determining the contributions 
due to the quantum corrections are all of the same order, b ~ c ~ d. In this case the constraint 
on Hi leads to a lower bound on the parameter b, given by 

10" 45 f H 1 Y 2 . 

h - " m - io) \wj ■ (6 - 2 - 33) 

As can be appreciated from figure (2) there is a region in parameter space where the resulting 
magnetic field strength is larger than B s > 10 -20 G corresponding to r > 10~ 37 , namely, 
(3 < —2.4, -^- > 10~ 18 bounded by the corresponding contour line. In the case of De Sitter 
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Fig. 2. In the case b = c = d the contour lines are shown for the maximum value of the logarithm of 
the ratio of magnetic to background radiation energy density. The values of log ( jj 2 - J correspond 

to reheat temperatures between 10 7 GeV and 10 19 GeV. The numbers within the graph refer to the 
value of log 10 r miB along the closest contour line. 

inflation, corresponding to /3 = — 1, there is no significant magnetic field generation since the 
mode functions during inflation as well as during the radiation dominated stage are plane 
waves. Furthermore, it can be checked that the resulting maximum magnetic field strength 
satisfies the bound due to gravitational wave production [107]. It was shown in [107] that 
for magnetic fields created before nucleosynthesis conversion of magnetic field energy into 
gravitational wave energy takes place. This leads to a maximal value of tqw given by, [272] 
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(6.2.34) 



at the galactic scale used here, A = 1 Mpc. Thus, the requirement r max < raw leads to an 
upper limit on jp- , that is, 
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(6.2.35) 



Hence the allowed range is given by (jjfM < (ifj ■ This * s & l wa y s satisfied since Hi < 
M Pl . 

It is also important to check that the fluctuation in the electromagnetic field during inflation 
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are within the perturbative regime and thus there is no strong backreaction on the dynamics of 
inflation. This effect can be estimated by calculating the energy density in the electromagnetic 
field and comparing it with the total energy density during inflation given by, 



3 /#!\ 2 ^r 2(/m) 



M Pl 8tt \MpiJ {rnj 
The average value of the electromagnetic field energy density is found to be [272] 



(6.2.36) 



4s(*r*erte)"*~ 

where fc* is the wave number corresponding to the scale which becomes superhorizon at the 
time T) during inflation. Thus with k* ~ — rj -1 the ratio (p(cm))(v)/P ~ P/Mp which is always 
smaller than one in the classical domain. Therefore, no backreaction effects have to be taken 
into account. In other models of magnetic field generation during inflation backreaction does 
play a role [276,277]. 



6.3 Trace anomaly 



Linear electrodynamics is scale invariant at classical level. Taking into account quantum cor- 
rections it is known that this classical symmetry is broken. Defining the energy momentum 
tensor by (see for example, [278]) 

Q ab = 2—^- f d 4 xC m , (6.3.1) 

then if the classical theory is scale invariant there will be a conserved current C a = Q ab Xb such 
that 



d a C a = e a a . (6.3.2) 

A scale transformation is equivalent to a conformal transformation of the metric, such as 

g mn {x) ->■ e 2a g mn (x). (6.3.3) 

Since the trace vanishes of the electromagnetic field in linear electrodynamics, the current 
C m is conserved at the classical level. However, when quantum corrections are included scale 
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transformations are no longer a symmetry, since the renormalized coupling constant depend 
on the scale. Namely, the renormalized coupling constant changes under the conformal trans- 
formation (6.3.3) as [278] 

g^g + a(3(g), (6.3.4) 

where (3(g) is the beta function. The Lagrangian changes as a(3(g)-^-C Thus the current 

satisfies, 

d m C m = QZ = /%)!"£■ (6-3.5) 

dg 

In massless QED the trace of the energy-momentum tensor can be found explicitly as [278] 

6™ = ^F ab F ab . (6.3.6) 

A similar expression for the trace of the energy-momentum is also found in QCD and other 
gauge theories. 

The trace anomaly was used in [279,280] (see also [281]) to study the generation of primordial 
magnetic fields during inflation. It induces a new term in Maxwell's equations, namely, in a 
flat Friedmann- Robertson- Walker background with scale factor a, [279,280], 

d m F n m + K^F n m = 0. (6.3.7) 

a 

The constant k depends on the theory which is used. For example, for the SU(iV) gauge theory 
with Nf generations of fermions in the fundamental representation [279,280], 

a fllN 2N f 



where a is the fine structure constant taken at the time of horizon crossing of the scale k~ x dur- 
ing inflation. Quantizing the gauge potential and finding the spectrum of the electromagnetic 
field in de Sitter inflation it is found that [280] 



\A k \~(^~y. (6.3.9) 

Thus for large values of k, corresponding to a large number of light fermions during inflation the 
trace anomaly could provide an efficient mechanism to generate large magnetic fields during 
inflation to serve as seed magnetic fields for a subsequent amplification by a galactic dynamo. 
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6.4 Coupling to other fields and varying couplings 



In [35] the coupling of a pseudoscalar axion to electrodynamics was proposed which for energy 
scales below the Peccei-Quinn symmetry breaking scale f a can be described by the effective 
Lagrangian, 

C = - l -d m 6d m 6 - - A F rnn F mn + g a 6F mn *F mn , (6.4.1) 

where g a is a coupling constant and the vacuum angle = <f) a / fa, where a is the axion field. In 
[282] a similar model has been considered in detail, namely the coupling of a pseudo Goldstone 
boson to electrodynamics. It is interesting to note that in these models the created magnetic 
fields have non zero helicity [282,283,284,285]. In [282] the Lagrangian is assumed to be of the 
form 

£ = ~\ (F mn F mn + g$F mn * F™) , (6.4.2) 

where g = a/(2irf) and / is the coupling constant and a the fine structure constant. The field 
equations in Fourier space are found to be [282] 

J2rp 

' ' "' ' '"" 0, (6.4.3) 




where F± = a 2 (B y ± iB z ) are the two circular polarization modes. The electric field satisfies 
the equation, [282] 
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d 2 G± f l2 d<j) 7 \^ d\ 

where G± = a 2 (E y ± iE z ). Using a potential of the form V((f>) = A 4 [l — cos(0//)] for the 
scalar field, the resulting amplification of the magnetic field during inflation is too weak in 
order to provide a seed field for the galactic dynamo [282]. In [286] a model with N pseudo 
Goldstein bosons has been investigated. It has been found that even in the case of one pseudo 
Goldstein boson due to the helical nature of the generated magnetic field the process of the 
inverse cascade will result in a strong enough seed field at the time of galaxy formation. 

A time-dependent coupling of the electromagnetic field provides another possibility of ampli- 
fication of perturbations in the electromagnetic field during inflation. This was first studied 
in [287] (see also, [288,289,290,291,292]) considering a lagrangian of the form [287] 



C ~ e a *F mn F mn , (6.4.5) 
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where a is a constant. Perturbations in the electromagnetic field are amplified during de 
Sitter inflation. The resulting magnetic field for the choice of a = 20 is found to be as large 
as 6.5 x KT 10 G today. 

In [293] (see also [294]) the generation of cosmologically relevant magnetic fields and their 
subsequent signature in the CMB has been discussed in a model where instead of the usual 
U(l)em gauge field, the photon, a hypercharge field Y m [295], which is associated with the U(l)y 
hypercharge group before the electroweak phase transition when the SU(2)xU(l)y symmetry 
is still unbroken, is coupled to a spectator field during inflation. After the electroweak phase 
transition the photon field is determined by the hypercharge field by A, m = Y m cos #w which 
gives rise to a primordial magnetic field in the postinflationary universe. 



6.5 Magnetogenesis in string theory 



A natural candidate for a scalar field coupled to the electromagnetic field is provided within the 
low energy limit of string theory. In the low energy limit string theory leads to Einstein gravity 
coupled to additional fields, such as the dilaton, which is a scalar field, and the antisymmetric 
tensor field which in four dimensions can be related to a pseudo-scalar field, the axion. To 
lowest order in the inverse string tension ol and in the loop expansion controlled by the string 
coupling g s the action in the so-called string frame is given by 



S = - ft n \ d D x^/-g D e 



Rd + d„ 



+ _ p T? mn 
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(6.5.1) 



where Gd is Newton's constant in D dimensions, Rd the Ricci scalar in D dimensions and 
is the dilaton. Indices take values between and D — 1. The string coupling is given by 
g s = e^. Superstring theory can be consistently quantized only in D = 10 and M-theory 
predicts 11 space-time dimensions. In order to reduce the resulting model to the four observed 
space-time dimensions, the extra space-time dimensions can either be treated as compactified 
to small extra dimensions, following the paradigm of Kaluza-Klein compactification, or one 
could model the observable universe as a four dimensional hypersurface embedded in a higher 
dimensional background space-time, which is the procedure followed in the models of brane 
cosmology [296]. Here the simplest model is used where the extra dimensions are compactified 
on static tori with small, constant radii. Thus the action (6.5.1) is used in D = 4 dimensions 
[297,298,299]. 

It is difficult to implement the standard slow roll inflation paradigm in string cosmology, 
derived from the low energy limit of superstring theory. The reason for that is that the dilaton 
does not have an appropriate potential. The potential resulting from supersymmetry breaking 
is far too steep to allow for a slow roll phase in the evolution of the dilaton. Inflation driven 
by the kinetic energy of the dilaton, however, can be realized. This is the pre-big-bang model 
[297,298,299] where inflation takes place for negative times (pre-big-bang phase) and is matched 
to the standard radiation dominated stage for positive times (post-big-bang phase). Since in the 
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low energy limit of superstring theory to lowest order Einstein gravity is recovered at cosmic 
time t — there is a space-time singularity, which follows from the theorems of Penrose 
and Hawking. Thus higher order corrections have to be included in order to regularize the 
transition between the pre- and post-big-bang era. In general, when calculating perturbations 
in pre-big-bang inflation it is assumed that the background evolves from an asymptotically 
flat initial state at t — > — oo to a high curvature phase at around t — which, however, 
never reaches a singular state. Only a few explicit, non-singular solutions are known and it 
seems difficult to determine the generic behaviour of pre-big-bang cosmologies [300]. During 
pre-big-bang inflation in four dimensions the scale factor behaves in the string frame in which 
the universe is expanding and accelerating, 




(6.5.2) 

assuming the end of inflation at r]i. The dilaton behaves in the low energy phase, r) < r] s , 

as [301] 

= — v3 In \r)\ + const. (6.5.3) 

After some time f] s higher order corrections in the inverse string length a' become important 
and the universe enters into a string phase which lasts until the end of inflation at r]\. During 
the string phase the dilaton evolves as [301] 



(f) = -2/3 In y+ const. (3 = - , , (6.5.4) 

2mz s 

where z s = a\/a s and a s and oi are the scale factors at the beginning of the string phase and 
at the end of inflation, respectively. 

Maxwell's equations derived from the action (6.5.1) imply in the radiation gauge in Fourier 
space the mode equation for the gauge potential, [301,302] 

Al + [k 2 - V( V )]A k = 0, V( V ) = g^g; 1 )" , (6.5.5) 

where a prime denotes the derivative with respect to conformal time rj. Matching the stage 
of pre-big-bang inflation to the radiation dominated era, quantizing the gauge potential and 
calculating the appropriate Bogoliubov coefficient results in the following values of the ratio 
of magnetic energy density and background radiation energy density [301], 

2 , ^ x 4-2/3 

r(w)~— 1- ( — ) u s < u < wi (6.5.6) 

1671^ \u>i/ 

i-J-Y'^z^e-^ u<u s , (6.5.7) 
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where g\ is the value of the string coupling at the beginning of the radiation era, A<p s = <p s — <fii 
and oj s = Ui/z. In figure 3 the ratio r(u) is plotted at galactic scale IMpc which corresponds 



(gi/4/r) 2 = 10 -3 , r>10 -57 
(gi/4w) 2 = 10 _1 , r>10~ 57 
(gi/4zr) 2 = 10 _1 , r>10~ 8 
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Fig. 3. The logarithm of the ratio r, magnetic field energy density over photon energy density, is 
shown at galactic scale, corresponding to IMpc, for different values of the string coupling at the 
beginning of the radiation dominated era, <?i, and the minimal value r m ; n necessary to seed the 
galactic magnetic fields after amplification by a galactic dynamo (r(wc) > ^min = 10~ 57 ) or without 



(r{u G ) > 



10 ). H\ ~ Mpi is assumed. The shaded areas indicate the region in parameter 



space for which the two conditions r{uja) > r m j n and r(u>) < 1 are satisfied. The dark grey area 
corresponds to the domain in the z s — g s parameter space for which the galactic magnetic field could 
be seeded directly without the assumption of a galactic dynamo. 

to ujg = 10~ 14 Mpc. Imposing that r{uj) < 1 for all frequencies, leads to the condition [301,302] 
z~ 2 < g s /gi- This implies a lower bound on the value of the coupling at the beginning of the 
string phase, g s . From figure 3 it can be appreciated that for a duration of the string phase 
determined by z s > 10 20 and a string coupling g s less than 10~ 42 the resulting magnetic field 
is strong enough to seed the galactic magnetic field directly. Furthermore, figure 3 shows that 
even for a very short string phase the resulting magnetic fields can be as strong as 10 _30 G which 
is the limiting value in case of action of a galactic dynamo in a universe with non vanishing 
cosmological constant. Without taking neither the string phase nor the cosmological constant 
into account, corresponding to a minimal required value of 10 -37 it was concluded in [303] 
that it is not possible to generate cosmologically relevant magnetic fields during pre-big-bang 
inflation. 

Even though in this section we focus on mechanisms which rely on the amplification of elec- 
tromagnetic perturbations during inflation, we will briefly comment on a different model of 
magnetogenesis in string theory. In [255] the generation of primordial magnetic fields from 
heterotic cosmic strings is studied. Heterotic fundamental cosmic strings were ruled out by 
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Witten for stability reasons [304]. However, as was shown in [305] the presence of branes 
offers a solution to the stability problem. In [306] heterotic cosmic strings are constructed 
by wrapping M5 branes around the 4-cycles of the Calabi-Yau manifold present in heterotic 
string theory. In [255] it was found that in a generalisation of the model of [306] the resulting 
heterotic strings are superconducting and as such can generate strong magnetic seed fields 
(see also section 5). 



6. 6 Magnetogenesis from extra dimensions 

Extra dimensions played a role in gravity ever since the proposal by Kaluza [307] to ex- 
plain gauge fields geometrically. Postulating a fifth dimension the components of the metric 
involving the fifth coordinate can be interpreted as the components of the gauge potential 
A m of electrodynamics and a scalar field 0, whose effective coupling to electrodynamics in 
four dimensions is that of a dilaton (see, e.g., [308]). Einstein's equations in vacuum in five 
dimensions imply Einstein's equations in four dimensions as well as Maxwell's equation for 
the gauge potential A m and the massless Klein-Gordon equation for the scalar field 0, if the 
dependence on the fifth coordinate is suppressed. However, despite its successful unification 
of gravity and electrodynamics there is still something missing in this picture. The point is 
how to explain that we have not observed the fifth dimension and why there is no dependence 
on the extra dimension. These problems were solved by Klein [309] assuming that the extra 
dimension is a circle of such a small radius that it is beyond observational limits. Using a 
Fourier expansion in the extra coordinate at each point in the four dimensional space-time 
there is an infinite number of four- dimensional fields. The zero mode results in the original 
theory of Kaluza where the fields have no dependence on the extra coordinate. The remaining 
part of the spectrum corresponds to massive modes. 

Extra dimensions appear naturally in models of string/M-theory which also admits solutions 
with large extra dimensions [310]. Contrary to the Kaluza-Klein picture in the case of large ex- 
tra dimensions inspired by string theory, our observable four dimensional universe is described 
by a four-dimensional hyper surface (brane) embedded in a higher dimensional background 
space-time. The cosmological solutions on the brane are influenced by the curvature of the 
higher dimensional space-time projected onto the brane. This leads for example to additional 
terms in the Friedmann equation [296]. 

In models derived from higher dimensional gravity the four dimensional Planck mass, which 
in this section will be denoted by M 4 , is no longer a fundamental parameter, but the IP- 
dimensional Planck mass Mp. Assuming for simplicity that all extra dimensions are of the 
same characteristic size R, the four-dimensional and the D-dimensional Planck masses are 
related by, 

M 4 2 = R n M™ +2 , (6.6.1) 

which follows directly using Gauss' law [311,312,313]. Newton gravity is observed to be valid 
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downto length scales of the order of 1mm [314]. This can be used to put a bound on Mo/M^. 

Dynamical extra dimensions can generate a primordial magnetic field in four dimensions. 
This rests on the fact that the conformal invariance of Maxwell's equations is broken in four 
dimensions in the presence of higher dimensions of time-dependent size. Assuming the metric 
to be of the form [315,316] 



ds 2 = -a 2 (ri) [drf - 5 a ^dx a dx 13 ] + b 2 (r])5 AB dy A dy B , (6.6.2) 

where a,/3 = 1, .., 3 and A, B — 4, .., 3 + n, n > 1. a(rj) and b(rj) are the scale factor of the 
external, 3-dimensional space and the internal, n-dimensional space, respectively. Assuming 
that before a time 77 = —771 inflation takes place in the external dimensions while the extra 
dimensions are collapsing. After this time the universe enters the standard radiation dominated 
era with the extra dimensions frozen to a small size. The first stage is described by a generalized 
vacuum Kasner solution. Thus the behaviour of the scale factors is determined by, 



airf) = d , b(rj) = bi , for 77 < —% (6.6.3) 

V ViJ \ ViJ 

, 6(77) = 61, for 77 > —771 . (6.6.4) 

m J 

In the following we set a± — 1 — bi. The Kasner exponents o and A are given in terms of the 
number of extra dimensions by [315], 



° = --[\ /— "I I A = J , 3 v (6.6.5) 

2\\jn + 2 J 1 V n(n + 2) v ; 

In D dimensions Maxwell's equations are given by V A F AB = with F^ B = Vr^A™, A,B = 
0, .., n + 3. Assuming that A^ = A^(x a , y B , 77) , A B = and using the radiation gauge A = 0, 
D^A^ = 0, Maxwell's equations imply 



- 1ft [b n d A,] + ]T d u d u A, + (^V ]T d B d B A, = 0, (6.6.6) 

v=l ^ tj ' / B=A 

where do = ^-, d^ = -£^ and d B = t^b-. Defining the canonical field \P M = 6^A^ the following 
expansion is used 



^(77,^,7/) = /^E^l) Ka^^e^ + al^^r^e-^] , (6.6.7) 
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where l m is a (3 + n)— vector with components / M = k^ 1 , l A = g^ 4 . Moreover, 1-X = kx + qy. 
a runs over the polarizations. During inflation for rj < —rji the mode equation is given by 
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(6.6.8) 



where ' = 4- and iV 

Of] 



j (nA — 1) — |. Moreover, /3 = <j — A. /3 < since only solutions with 

contracting extra dimensions are of interest here. — 1 < (3 < —1/(1 + \/3), where the lower 
boundary corresponds to n — 1 and the upper bound gives the value for large n. For n — 1 
and n = 6 there are known solutions in closed form of equation (6.6.8). Hence for n = 1 the 
exact solution is used to find the spectrum of the primordial magnetic field and for n > 1 
approximate solutions are employed. 

The solution for one extra dimension, n — 1, is given by [316] 
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(6.6.9) 



where i^ 2 •* (z) is the Hankel function of the second kind. 



The approximate solution for n > 1 is found by solving the mode equation (6.6.8) in two 
regimes determined by whether the term due to the modes q in the extra dimensions, (—f]/rji) 2 ^q 2 
is larger or smaller than k 2 , where k are the comoving wave numbers in the observable three 
dimensional space. When the contribution due to the wave numbers q in the extra dimensions 

is subdominant ( — — J q 2 < k 2 , or u q < Uk in terms of the physical frequencies Uk = k/a(rj) 
and u q = q/b(r)), the canonical field is approximately given by [316], 
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where H ^ is the Hankel function of the second kind and /i 2 = l; + iV=^/i=l (nX — 1 ) . In 
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the other case, that is for [ — JL ) q > k , or u q > Uk, it is found that [316] 
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(6.6.11) 



where k = ^- and /i = \{n\ — 1). 

During the radiation dominated era rj > —rji the mode equation is given by 
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(6.6.12) 



which is solved in terms of parabolic cylinder functions E(a, z) [315], 
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(6.6.13) 



where z = f — J 2 (77 + 2771) and a = 
satisfying the normalization |c + | 2 — 
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. Furthermore, |c±| are the Bogoliubov coefficients 



The total magnetic energy density is given by [317] 
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(6.6.14) 



where, assuming that the volume consists of two spheres, dV = -^j-^^-k dk A T^wrq n dq. 
At i] = —i]i the comoving wavenumbers k and q are equal to the physical momenta, since 
a\ = 1 = b\. The spectral energy density p{oJk) = dp/dlogUk is then given by 
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(6.6.15) 



where Y = - 1 , and c^ = -, w g = |. To calculate the ratio of energy density in the magnetic 
field over the background radiation energy density it will be assumed that U\ — — and 
k\ ~ Hi is the maximal wave number, leaving the horizon at the end of inflation and thus at 
end of the dynamical higher dimensional phase rj = —rji- Equally it is required that there is a 
maximal wave number q max in g-space corresponding to the modes in the extra dimensions. 
These assumptions are justified by the sudden transition approximation used here. At the 
time of transition i] = —r]i the metric is continuous but not its first derivative. For modes 
with periods much larger than the duration of the transition the transition can be treated as 
instantaneous. However, in order to avoid an ultraviolet divergency an upper cut-off has to be 
imposed [318,319,320]. The following expressions for r(uk) are obtained [316]: 



(1) For q — 0, n — 1 
r(uj k ) 



3tt 2 VM 4 
(2) For q > and n = 1 



2 /^V^V ln ^y 

OJiJ 



UJlJ 



(6.6.16) 
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where o; (ri) = ^^- and it was assumed that u a > Uk- 
(3) For q = and n > 1 



2nA-2 



<-(*¥)*-«r®'(*r- — » 



Furthermore, n\ = J^r^- Since nA < 4 the resulting spectrum for r(uk) is increasing in 
frequency. 
(4) For q > and n > 1 

r W ~V„»*~ (|L)" +2 (^ )- 2 <" (^ (6 ,, 9) 

where 

A r 16 8tt 7t 1+ § 2 2 ^" 3 _ 2 . . ^2^/ l^ 2 

where subleading terms have been omitted and uj qmax > oou was assumed. The resulting 
spectrum is growing in frequency. 

The resulting spectrum of the primordial magnetic field is characterized by the Hubble pa- 
rameter at the beginning of the radiation dominated era Hi, the Z)-dimensional Planck mass 
Md and the number of extra dimensions, n. In addition, in the case where the modes lying 
in the extra dimensions are taken into account, there is the maximal physical frequency u qmax 
which is estimated assuming q ma x ~ fci- The spectrum is constrained by r(u) < 1 for all 
frequencies. The ratio of the D-dimensional over the four dimensional Planck mass is limited 
by the observation that Newtonian gravity is valid at least down to scales of the order of 
1 mm [314]. This leads to the lower bound -jf- > (1.616 x 10~ 32 )"+2. Furthermore, with 2~i 
the temperature at the beginning of the radiation epoch, big bang nucleosynthesis requires 

that T\ > 10 MeV. This imposes a bound on Hi by using |j- = 1.66g* (Ti) (j^-) , where for 
T\ > 300 GeV the number of effective degrees of freedom is given by g*(Ti) = 106.75 (see, 
e.g., [163]), namely, log g > -40.94. 

Imposing the various constraints leads to upper limits on r(uk) calculated at the galactic 
scale corresponding to 1 Mpc, that is, ug = 10~ 14 Hz and taking the maximally amplified 

frequency evaluated today to be uj\ ~ 6 x 10 11 Hz (jjfM 2 [316]. Not taking into account modes 
in the extra dimensions leads for one extra dimension to magnetic field strengths B s < 10~ 39 
G. However, taking into account these modes substantially increases the upper value of the 
magnetic field to upto 10~ 8 G. Imposing the constraint T\ ~ M 5 leads to magnetic seed fields 
B s < 10~ 23 G. In models with more than one extra dimension, n > 1, strong magnetic seed 
fields can be created if the internal momenta are taken into account. In particular, without 
the assumption that the temperature at the beginning of the radiation epoch is of the order of 
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the .D-dimensional Planck scale allows for the creation of seed magnetic fields with strengths 
of upto 1CT 10 G. For more than three extra dimensions, this also holds assuming 7\ ~ Mjy. 
With this assumption for two and three extra dimensions results in weaker magnetic seed 
fields, with maximal field strengths, B s < 1CT 18 G for two extra dimensions and B s < 1CT 13 G 
for three extra dimensions. 



6. 7 Magnetogenesis in theories with broken Lorentz symmetry 



The spontaneous breaking of Lorentz invariance is present in certain solutions of string field 
theory which leads to a non-vanishing photon mass described by the Lagrangian [321] 



C = - l -F mn F mn + M 2 L a- 2e A m A m , (6.7.1) 

2 

where M£ = -tjw, ttil is a light mass scale in comparison with the typical string energy scale 
m pi 

and 2£ is a positive integer. In [321] cosmo logically interesting magnetic field strengths are 
found for a diverse choice of parameters of the model. In [322] an extension of the standard 
model is presented in which due to new physics at the Planck scale Lorentz symmetry is broken 
spontaneously. In the pure photon sector of the extended QED the Lagrangian is given by, 

[322] 
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£ = —jF mn F mn — -{k F ) p i mn F p F mn + -(k AF ) p e p i mn A F mn , (6.7.2) 

where the coupling (kF) P imn is real and dimensionless and (/caf) p is real and has dimensions 
of mass. In the context of the generation of primordial magnetic fields during inflation the 
Lagrangian (6.7.2) has been investigated in [323,324]. Analyzing the model resulting by taking 
into account only the first two terms in (6.7.2) it has been shown in [323] that magnetic fields 
of nanogauss field strength on a megaparsec scale at present can be generated for a wide 
range of parameters. In [324] (see also, [325]) primordial magnetogenesis during inflation has 
been discussed in a model resulting from considering only the first and third term in (6.7.2). 
In this case the generated magnetic field is found to be maximally helical at the end of de 
Sitter inflation. The subsequent inverse cascade of the magnetic field spectrum taking place 
in the turbulent plasma during the radiation dominated era results in a magnetic field with 
an interesting field strength and correlation length at the time of the protogalactic collapse. 

Noncommutativity in space provides a different possibility of breaking Lorentz invariance 
which in this case is explicitly broken so that all the amplitudes are frame dependent. In the 
context of generation of primordial magnetic fields this was first discussed in [326] . Noncom- 
mutative spaces occur in string theory in the Seiberg-Witten limit [327] and are described by 
the commutation relation for the coordinate operators x m , 
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[x m ,x n ] =i9 mn , (6.7.3) 

where 9 mn is a constant of dimension length which is conveniently parametrized in terms of 
the noncommutativity scale Anc, defined by 9 mn = 4r— where c mn is an antisymmetric tensor 
with components of order unity [326]. Moreover, in order to avoid problems with unitarity 
and causality, 9 0fl = is chosen such that only space is noncommutative. In [328] it was 
shown that the magnetic dipole moment of a charged massive particle, such as the electron, 
receives quantum corrections at one loop which are spin independent and proportional to 
9^ = e^ UK 9 UK . This leads to a non vanishing magnetic field proportional to 9^ when summing 
over all possible states. However, choosing the noncommutativity scale A^c — 10 3 GeV the 
authors find the resulting magnetic field to be too weak in order to successfully seed the 
galactic dynamo. In [329] noncommutative quantum field theory was used for the U(l) gauge 
field leading to a modified Lagrangian describing the photon which is of the form of the 
Lagrangian (6.7.2) including the first and the third term. Moreover, in this case (&UF)m is 
nonzero only for the spatial components and given by the noncommutativity parameter 9^. 
Using the approach of [330] to implement the stringy spacetime uncertainty relation which 
leads to an effective noncommutative space-time [331] investigate primordial magnetogenesis 
in dilaton electromagnetism. In [332] the generation of primordial magnetic fields in inflation 
with a cut-off is investigated. The effect of the cut-off is to add extra terms to the action which 
in the model under consideration describes a photon with a mass term during inflation. The 
free parameter of the model can be chosen such that cosmologically relevant magnetic fields 
are obtained. 



6.8 Magnetogenesis and nonlinear electrodynamics 



So far the models in this section proposed to generate primordial magnetic fields are all situ- 
ated within linear electrodynamics. In order to amplify perturbations in the electromagnetic 
field during inflation the electromagnetic field is coupled to a scalar field or curvature terms, 
quantum corrections resulting in the trace anomaly, symmetries are broken or dynamical extra 
dimensions are taken into account. All of these leading to the breaking of conformal invariance 
of Maxwell's equations in four dimensions. 

Nonlinear electrodynamics provides yet another possibility of breaking conformal invariance 
of the electromagnetic field. On large scales present day observations confirm the linearity in 
the electric and magnetic fields of Maxwell's equations in vacuum. However, as smaller and 
smaller scales are approached one might expect deviations from linearity due to the fact that 
charges become more localized (see e.g. [127]) and hence increases the energy density. This 
led to the hypothesis that there is some upper bound on the field strengths avoiding thus an 
infinite self-energy of a charged particle. A first example of a classical singularity-free theory 
of the electron was proposed by Born and later by Born and Infeld [333,334,335,336]. The 
modified field equations can be derived from the Lagrangian of the form [334] 



82 



C = b 2 (l-^/l-(E 2 -B 2 )/b 2 ), (6.8.1) 

where b is a maximal field strength. In this section vector notation will be used to make 
expressions easier to read. The electromagnetic field is modified at short distances and its 
energy density is finite. One of the problems with this type of theory is its quantization 
[127]. Nonlinear electromagnetism had been considered before by Mie [337]. However, it was 
discarded since it depended on the absolute values of the gauge potential [334]. 

Another place where nonlinear electrodynamics arises is in quantum electrodynamics. Virtual 
electron pair creation induces a self-coupling of the electromagnetic field. Heisenberg and Euler 
calculated the self-interaction energy for slowly varying, but arbitrarily strong electromagnetic 
fields [338,339,340]. It is described by the lagrangian [338] 



C = -X + — dss- 3 e- m ^ s 
8ir 2 jo 



f x 2v -Re cosh esM 2 2 

im cosh esM 3 



(6.8.2) 



where 



X = - A F mn F mn Y = - A F mn *F mn (6.8.3) 

and M 2 = 2X - %Y . Expanding the lagrangian leads to [338,339,340] 

L = X + k X 2 + K t Y 2 . (6.8.4) 

This describes the Heisenberg- Euler theory for the choice n = J^j and /ti = ^j , where a is 
the fine structure constant and m e the electron mass. Furthermore, the propagation of a pho- 
ton in an external electromagnetic field can be described effectively by the Heisenberg-Euler 
langrangian. Moreover, the transition amplitude for photon splitting in quantum electrody- 
namics is nonvanishing in this case. Photon splitting is a process in which an electron-positron 
pair is created and one of the particles emits a photon before annihilating with the other par- 
ticle to generate the second photon. Thus the initial one photon state transforms into a two 
photon final state. In principle, this might lead to observational effects, e.g., on the electro- 
magnetic radiation coming from neutron stars which are known to have strong magnetic fields 
[340,341,342]. In particular, certain features in the spectra of pulsars can be explained by 
photon splitting [343,344]. 

Finally, Born-Infeld type actions also appear as a low energy effective action of open strings 
[345,346,347,348]. As was shown in [349] the low energy dynamics of D-branes is described by 
the Dirac-Born-Infeld action. 

To test whether nonlinear electrodynamics can lead to the generation of cosmologically rele- 
vant primordial magnetic fields the following model will be considered. A stage of de Sitter 
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inflation followed by reheating is matched to a standard radiation dominated era. During in- 
flation quantum fluctuations are excited within the horizon. Upon leaving the causal domain 
they become classical perturbations. It is assumed that electrodynamics is nonlinear during 
inflation and becomes linear once the universe enters reheating and subsequently the radia- 
tion dominated stage. This latter assumption ensures that the evolution during the radiation 
dominated era and subsequent stages of the universe are described by the standard model of 
cosmology. To study nonlinear electrodynamics in this setting was put forward in [350] and 
independently in [351]. 



6.8.1 Field equations 

The Born-Infeld or Heisenberg-Euler lagrangians are particular examples of theories of non- 
linear electrodynamics. In general the action of nonlinear electrodynamics coupled minimally 
to gravity can be written as, see e.g. [347,348,352] 



S = — j^- J d A x^R + -^ J ' d 4 x^gL(X,Y), (6.8.5) 

where X and Y are defined by equation (6.8.3). Maxwell's electrodynamics corresponds to the 
choice L = —X. The equations of motion are given by 



V m P mn = , (6.8.6) 

where P mn = — (LxF mn + Ly *F mn ), where the dual bi- vector *F mn is given by 

*pmn _ _L e m,nabp ab anc } e mnab ^ e Levi-Civita tensor with e i23 = +1- Furthermore La 

denotes La = dL/dA, and 



V m *F mn = 0, (6.8.7) 

which implies that F mn = d m A n — d n A m . Assuming the background metric to be of the form, 

ds 2 = a 2 {r]) [-drf + dx 2 ] . (6.8.8) 

Furthermore writing the Maxwell tensor in terms of electric and magnetic fields in the "lab" 
frame (cf. equation (6.1.2)) equations (6.8.6) and (6.8.7) imply [350], 
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(6.8.9) 

(6.8.10) 

(6.8.11) 
(6.8.12) 



From these equations two wave type equations can be derived which, however, contrary to 
the case of linear electrodynamics do not decouple the electric and magnetic field. Taking the 
curl of equation (6.8.10) and using equations (6.8.11) and (6.8.12) a wave type equation for 

the magnetic field B can be found [350]. 
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(6.8.13) 



Similarly, taking the time derivative of of equation (6.8.10) and using the remaining equations 
results in a wave type equation for the electric field E [350], 
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(6.8.14) 



In the long wavelength approximation spatial gradients can be neglected [353] . Thus neglecting 
spatial derivatives equation (6.8.13) implies, 



1j X J^x J^x 
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(6.8.15) 
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where Bk = a 2 Bk, £k = a 2 Ek and a prime denotes the derivative with respect to conformal 
time 77, that is ' = j-. Assuming that the lagrangian depends only on X, that is Ly = 0, 
equation (6.8.15) implies 



k — T > 

L*X 



(6.8.16) 



where Kk is a constant vector and Lx 7^ 0. For Kk = linear electrodynamics is recovered, 
for which Bk = const. In the long wave length limit, the wave like equation for the electric 
field, equation (6.8.14), yields to 



d_ 
drj 



L x 
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Ly 
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(6.8.17) 



Thus integrating equation (6.8.17) results in 
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(6.8.1* 



where Pk is a constant vector. The equations determining the magnetic and electric field, 
(6.8.15) and (6.8.17), are coupled nontrivially for Ly 7^ 0. Thus in order to find solutions, 
the Lagrangian will be considered to be only a function of X, L — L(X). Furthermore, since 

X = \{B 2 — E 2 ) it is useful to find equations for £ k 2 and Bj 2 which are given by, for Pj* > 0, 
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(6.8.19) 
(6.8.20) 



Assuming that the constant vector in equation (6.8.18) vanishes, Pk = 0, leads to a significant 
simplification. In this case, equation (6.8.18) for L = L(X) can be solved immediately, giving 
for the electric field 



z M k 

Z-k — 7 — , 
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(6.8.21) 



where Mk is a constant vector. Thus for Pk = equation (6.8.20) leads to an equation only 
involving X and Lx, namely, 
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6.8.2 A particular model 

In order to find explicit solutions of equation (6.8.22) a particular lagrangian has to be chosen. 
For simplicity the lagrangian is chosen to be of the form 




(6.8.23) 

where 5 is a dimensionless parameter and A a dimensional constant. This is the abelian 
Pagels-Tomboulis model [354,355]. An effective model of low energy QCD is provided by its 
nonabelian version [356]. Clearly, linear electrodynamics is recovered for the choice 5 = 1. 
The lagrangian (6.8.23) is chosen since it leads to a significant simplification of the equations, 
but still allows to study the effects of a strongly nonlinear theory of electrodynamics on the 
generation of primordial magnetic fields. In general, the energy-momentum tensor derived 
from a lagrangian L(X) is given by 



T =!. 

4"7T 



Lx9 a F ma F bn + g mn L . (6.8.24) 

Furthermore, for the lagrangian (6.8.23) the trace of the energy-momentum tensor is given by 



T = - — -L, (6.8.25) 

IT 

which vanishes only in the case 5 = 1 that is for linear electrodynamics. The energy-momentum 
tensor is calculated explicitly to check whether there are any constraints on the parameter 5. 
Decomposing the Maxwell tensor with respect to a fundamental observer with 4-velocity u m 
into an electric field E and a magnetic field B, implies [357,114,115,119], 



where rj mn k s = \/—g^mnks and u m u m = — 1. Thus the electric and magnetic field are given, 
respectively, by E m = F nm u n and B m = \r\ mn uu n F kl . The lab frame is defined by the proper 
lab coordinates (t, r) determined by dt = aalrj, air = aolx. Applying a coordinate transformation 
then gives the relation between the fields measured by a fundamental observer and the lab 
frame. Using the four velocity of the fluid u m = (a -1 , 0, 0, 0) results in the relation [358] 



E„ = aEp, B„ = aB^. (6.8.27) 

As shown in [357,114,115,119] the energy-momentum tensor of an electromagnetic field can be 
cast into the form of an imperfect fluid. The energy-momentum tensor of an imperfect fluid 
is of the form (see for example, [357,114,115,119]), 
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where p is the energy density, p the pressure, g m the heat flux vector and 7r mn an anisotropic 
pressure contribution of the fluid. h mn = g mn + u m u n is the metric on the space-like hyper- 
surfaces orthogonal to u m . With q m u m = and n mn u m = 0, 



n = T v m v n a = -T v m h n 

l J - L mn u ' •* <ia ± mn Llj lb a 

Q ab = T mn h™h n b Q ab = ph ab + n ab . (6.8.29) 

Therefore using equations (6.8.24) and (6.8.26) the energy density and the heat flux vector 
for the Pagels-Tomboulis model (6.8.23) are found to be 

P = ~ k? [(2^ - 1) E a E a + B a B a ] (6.8.30) 

qa = ^VabmnU b E m B n . (6.8.31) 

Imposing the condition that ir a h is trace-free then the pressure and it a b are given by 

15-1 , 

p = p - L 6.8.32 
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- (E m E m + B m B m ) h ab — (E a E b + B a B b ) 



(6.8.33) 



Thus considering p (cf. equation (6.8.30)) in general there is a constraint on 5 which is 5 > | 
required by the positivity of the energy density p. 



6. 8. 3 Estimating the primordial magnetic field strength 

During de Sitter inflation electrodynamic is nonlinear and described by the Pagels-Tomboulis 
lagrangian (6.8.23). Thus in the very early universe electrodynamics is highly nonlinear and 
very different from standard Maxwell electrodynamics. At the end of inflation electrodynamics 
is assumed to become linear so that the description of reheating and the subsequent radiation 
dominated stage are unaltered. 

Recalling that the scale factor during de Sitter is given by a{rf) — a± ( - 2 -J , where 77 < 771 < 0. 
The end of inflation is assumed to be at rj = r]\. The equations determining the electric 
and magnetic field in the long wave length limit, (6.8.19) and (6.8.20), are coupled, since X 

depends on E 2 and B 2 , in particular the invariant X reads, 2a 4 X ~ Bj 2 — ££. Therefore to 
find approximate solutions three different regimes will be considered. 

(1) Bj ~ 0{£ 2 ). 

(2) Bl ^> £ fc 2 . This implies the approximation 2a 4 X ~ Bj*. 



(3) ££ 3> B%. This implies the approximation 2a' 



A X ~ -El 



Following [35] it is assumed that quantum fluctuations in the electromagnetic field lead to 
initial electric and magnetic fields. The energy density at the time of first horizon crossing 
during inflation, say at a time 772, is estimated to be of the order of the energy density of 
a thermal bath at the Gibbons-Hawking temperature of de Sitter space. Furthermore it is 
useful to recall the energy density in the magnetic field at the time of first horizon crossing 

corresponding to a value of the scale factor a 2 given by pb(&2) — H 4 ~ ( -^- ) , where M 4 is 

the constant energy density during inflation. Assuming that initially the magnetic and electric 
energy densities are of the same order, there is an equivalent expression for the electric energy 
density at the first horizon crossing. While after the end of inflation, during the radiation 
dominated epoque, the electric field rapidly decays due to plasma effects, the magnetic field 
remains frozen-in. 



;i) El ~ 0{£i): 

In this case equation (6.8.22) can be approximately solved and leads to the magnetic 
field strength at the end of inflation corresponding to the value of the scale factor 01 
determined by [350], 

B k( a i) „ 4NW cos^\rn{xi + (5- l)Ci)] 

Blip*) ~ cosh 2 [m(x 2 + (0 - l)d)] ' l ' ' ' 

where N(X) is the number of e-folds before the end of inflation at which A left the horizon, 
that is, e N ^ = ai/oo- Moreover, m = : J ,7-. . and x = tj=t- 

11 A ' M P i\M k \ M p l 

Furthermore, the constant C\ is chosen such that (5 — l)Ci = — x 2 . Using that during 
de Sitter inflation, a = 01(771/77) and the number of e-folds, results in the magnetic energy 
density ps at the end of inflation, 

Pb(oi) ^ p B (a 2 )e' 4Nix) cosh^-mx^^ - 1)], (6.8.35) 

where pp. = §-■ Together with the expression for the number of e-folds remaining until 
the end of inflation after the comoving scale A has crossed the horizon during inflation 
(cf. equation (4.3.1)), this results in the ratio of magnetic to background radiation energy 
density r at the end of inflation [350], 



r(a±) ~ 10 



-104 



A \ 4 /M\t 




MpcJ \T RH ) 




9.2xl0 25 f A ) 
\MpcJ 


/My /T RH y 

\M Pl ) \M Pl ) 



x cosh^ -9.2 x 10^ — — \{—^) \^ L ) mil 



(6.8.36) 



There are bounds on the parameter m coming from the requirements that r should be 
larger than a lower value in order to be strong enough to seed the galactic magnetic 
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field and smaller than an upper bound coming from the fact that r has be less than one. 
Therefore assuming A = 1 Mpc and the inflationary energy scale M = 10 17 GeV, the 
following values for —mx\ are found [350]. For a model with reheat temperature Trh — 
10 9 GeV [35] the parameter —mx\ has to be in the range 2.7 x 1CT 20 < —mx\ < 5 x 1CT 20 
in order to achieve a magnetic seed field with a field strength to be at least B s ~ 10~ 20 G, 
corresponding to r$ = 10~ 37 . For a higher reheat temperature Trh = 10 17 GeV [35], for 
the same magnetic seed field strength — mx\ has to be in the range, 9.5 x 10~ 23 < —mx\ < 
1.5 x 1CT 22 . And similarly, for the less conservative bound r = 10~ 57 , for Trh = 10 9 
GeV, —mx\ has to be in the range 1.4 x 1CT 20 < —mxi < 5 x 1CT 20 and for Trh = 10 17 
GeV it is found that 6.7 x 1(T 23 < -mx 1 < 1.5 x 1(T 22 . 

(2) Bl » Si 

For 5 7^ \ and 5 ^ | the r is found to be [350] 
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For P% > 0, solutions consistent with the assumption used in the approximation yield 
cosmologically interesting magnetic fields for 5 > 1.9 for Trh = 10 17 GeV and for 5 > 3.0 
for Trh = 10 9 GeV. In these cases the ratio of the energy density of the magnetic field over 
the energy density of the background radiation r is found to be r > 10 -37 corresponding 
to a primordial magnetic field of at least B s = 10 _20 G. 

In the cases 5 = \ and S = | the solutions show a behaviour different from the power- 
law solutions for X used in the former case (cf. equation (6.8.37)). On the one hand, for 
S = | an implicit solution depending on Euler's (3 function can be found which however 
makes it difficult to estimate explicitly the magnetic field strength at the end of inflation. 
On the other hand, for S — | a solution depending on simple functions can be found, 
namely the ratio r at the end of inflation is given by, [350] 
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Here the constant C\ has been chosen as Ci = — (ff) 2;r i ( — ) • Furthermore, a.2 — 

4i? 2 

M -2 .4 4 . In this case the solutions found for the electric and magnetic field are consistent 

m pi a a i 

with the approximation for P% > and P 2 = 0. Moreover, the resulting magnetic field is 
strong enough to seed the galactic dynamo [350]. 

(3) El » Bl 

In this case the approximation implies that at the end of inflation the electric energy 
density is much larger than the magnetic energy density. 
The ratio r at the end of inflation is given by [350] 
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For P% > the resulting magnetic fields very weak magnetic fields and r(ai) <C 10 -37 for 
typical values of the cosmlogical parameters. However, in the case P% = 0, for 5 > 19.5 
and a reheat temperature Trh = 10 9 GeV primordial magnetic fields result which could 
successfully act as seed fields for the galactic dynamo. 

In summary, using the Pagels-Tomboulis model an example of a theory of nonlinear elec- 
trodynamics has been provided which can lead during inflation to sufficient amplification of 
perturbations in the electromagnetic field in order to seed the galactic magnetic field. 

In [351] the resulting magnetic field generated during inflation is estimated for Lagrangians 
of the form C = C(X). The magnetic field strength is obtained by neglecting the magnetic 
field contribution to X on superhorizon scales since equation (6.8.12) implies B^ ~ krjE^ on 
these scales. Moreover, neglecting the spatial gradient terms in equation (6.8.10), [351] find 
the scaling 



(£ X ) 2 X oca~ 4 . (6.8.40) 

Using this relation for Lagrangians of the form C = —X + J27=2 c i-^- % an d £ — —X exp (— cX), 
where q and c are constants, and the Born-Infeld Lagrangian (6.8.1), the resulting magnetic 
field strength at present is estimated. There is a range of parameters for which magnetic fields 
strong enough to directly seed the galactic magnetic field can be generated. 

In [359] the generation and evolution of primordial magnetic fields has been discussed during 
de Sitter inflation, reheating and the radiation dominated era in theories of nonlinear electro- 
dynamics described by Lagrangian densities C ~ X + 7X 5 and C ~ X + fi 8 /X, where 7, 5 
and \x are constants. It was found that not only primordial magnetic fields of interesting field 
strengths can be generated but also the baryon asymmetry by gravitational coupling between 
the baryon current and the curvature of the background. In [360] primordial magnetic field 
generation has been discussed in DBI inflation. 
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7 Summary and outlook 



Although the origin of cosmic magnetism is still the subject of debate, the ubiquitous presence 
of large-scale 5-fields with similar (of /xG order) strengths in galaxies, galaxy clusters and high- 
redshift protogalactic structures seems to suggest a common, primordial origin for them. Very 
recent reports indicating the presence of extragalactic magnetic fields close to 10~ 15 G in low 
density regions of the universe also point towards the same direction. 

The possibility of a cosmological (pre-recombination) origin for all the large-scale magnetic 
fields is a relatively old suggestion and there have been numerous studies looking at the gen- 
eration, the evolution and the potential implications of such primeval fields. There are still 
serious difficulties to overcome, however, especially when trying to produce the initial 5-fields 
that will seed the galactic dynamo. Primordial magnetogenesis is still not a problem-free exer- 
cise, which probably explains the plethora of mechanisms proposed in the literature. Roughly 
speaking, magnetic seeds produced between inflation and recombination are too small in size, 
while those generated during inflation are generally too weak in strength. In either case, the 
galactic dynamo will not be able to operate successfully. The former of the aforementioned 
two problems is essentially due to causality, which severely constrains the coherence length of 
almost every 5-field produced during the radiation era. The latter problem is attributed to 
the dramatic depletion suffered by typical inflationary magnetic fields. Primordial turbulence 
and magnetic-helicity conservation can in principle increase the initial coherence scale of mag- 
netic seeds, especially of those generated during phase-transitions in the early radiation era. 
Considerable effort has also been invested in the search for viable physical mechanisms that 
could amplify weak inflationary magnetic fields. Solutions to the magnetic strength problem 
are typically sought outside the realm of classical electromagnetism and/or that of standard 
cosmology, although conventional amplification mechanisms can also be found in the litera- 
ture. The aim of this review is to provide an up-to-date and as inclusive as possible discussion 
on the current state of primordial magnetogenesis. 

Deciding whether the large-scale magnetic fields that we observe in the universe today are of 
cosmological origin, or not, would be a step of major importance for cosmology. If confirmed, 
such primordial fields could have affected in a variety of ways a number of physical processes 
that took place during the early, as well as the subsequent, evolution of the cosmos. Although 
the argument in favour of cosmological 5-fields may not settle unless an unequivocal magnetic 
signature is detected in the CMB, their case gets stronger as more reports of magnetic fields 
at high-redshifts and in empty intergalactic space appear in the literature. Upcoming obser- 
vations may also help in this respect. A new generation of radio telescopes, like the Expanded 
Very Large Array (EVLA), the Low Frequency Array (LOFAR), the Long Wavelength Array 
(LWA) and the Square Kilometer Array (SKA) have large-scale magnetic fields in their lists 
of main targets. If nothing else, the expected influx of new data should put extra constrains 
that may allow us to distinguish between the various scenarios of magnetic generation and 
evolution. Information of different type, but of analogous importance, may also come from 
CMB observations, like those associated with the ESA PLANCK satellite. At the same time, 
structure-formation simulations are becoming more sophisticated by the day and a number of 
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research groups have started systematically incorporating magnetic fields into them. This in 
turn should help us understand and interpret better the non-thermal regime of galaxy forma- 
tion. So, hopefully, we will soon have cosmological and structure formation models with fewer 
free parameters and more physics. 
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